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MECHANICS OF WASHOUT OF AN ERODIBLE FUSE PLUG 


By E. Roy Tinney,! M. ASCE, and H. Y. Hsu,? F, ASCE 


SYNOPSIS 


aboratory and field experiments are described that were undertaken to 
onstrate the feasibility of afuse plug inthe spillway of a major dam. From 
jJata and from tractive force theory, similitude relationships are derived 
permit the application of model test results to rates of washout. 


INTRODUCTION 


‘he Idaho spillway of the Oxbow project, on the Snake River between Idaho 
Oregon, will have a fuse-plug control instead of conventiona) gates. The 
otype fuse plug as designed and constructed for controlof this spillway will 
ide an economical, safe, and effective control for the maximum designed 
of 150,000 cfs. The combined design capacity of the fuse-plug spillway and 
sonventional spillway is 300,000 cfs. 

he fuse plug is simply a rockfill dam with special arrangement of the im- 
ious core and strict control of the rockfill and filter materials. It hasa 
or of safety equal to that for ordinary rockfill or earth dams when it func- 
3as adam; however, it will be washed out by the reservoir water auto- 
cally at a suitable rate when the spillway capacity is required. 
he general arrangement of the fuse plug controlled Idaho spillway is shown 
ig. 1. The fuse plug, located across the upper end of the approach channel, 
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a total length of 440 ft and varies inheight from 27 ft to 29 ft above a level 
crete sill. 
\t the upstream end of the embankment, a pilot channel, 10 ft wide with an 
rt 3 ft lower than the crest, is provided to start the washout of the fuse 
-at a predesignated location. The location of this pilot channel is chosen 
hat the flow of water, during washout, will be smoothly channeled into the 
Iway chute. 
Vhen it is required to pass an extraordinary flood exceeding the capacity of 
Oregon spillway, the fuse plug would be washed out to provide additional 
city. To initiate the fuse-plug breach, the Oregon spillway gates would be 
ed slightly raising the reservoir water surface to the required predeter- 
ed level, one foot above the invert of the pilot channel. After breaching, 
fuse plug would wash out progressively from the upstream to the down- 
am end in an orderly manner. During the washout, the Oxbow reservoir 
ld be kept more or less at constant level by closing the Oregon spillway 
Ss, thus shifting the discharge from the Oregon spillway to the Idaho spill- 
‘with no appreciable variation of outflow. 
‘use-plug controlled spillways have previously been designed and construc- 
The Aluminum Company of America has seven such spillways in service 
orth Carolina which were constructed under the directionof James P. Grow- 
F. ASCE, the originator and pioneer inthe use of the fuse plug for the con- 
of spillways in the United States. Fuse-plug sections have been used on 
es for passing extraordinary floods on major rivers in this country and 
mad. In case of unusual floods, the fuse plug is designed to wash out either 
matically or artifically to minimize flood damages. 
‘iS most existing fuse-plug spillways have not been in service long enough 
ass their design floods, thereby demonstrating their washout performance, 
‘ature describing the mechanics of washout of a fuse plug is insufficient to 
ride practical data substantiating design assumptions. Although scattered 
rmation is known from observations of washout of cofferdams and levees, 
guld be difficult to arrive at any reliable statistical deductions as to the 
nout characteristics. A project as large asthe Oxbow project on the Snake 
r required more thorough investigation. 
efore the fuse-plug spillway could be approved for construction, several 
amental engineering aspects had to be thoroughly investigated, sothat there 
d be no doubt that the fuse plug would function as designed. 


. How does the pilot channel breach? What is the water depth over the 
- channel before it breaches ? 

. Whatis the nature of the washout process after the initial breach? What 
ie speed of washout? 

. How complete is the washout? Will any portion of the plug remain in- 
? 


he fuse plug was designed according to theoretical computations based 
ily upon bed load transportation principles and modern concepts of soil 
hanics. However, the assumptions made by the design engineers required 
fication by a large scale model. Furthermore, some basic information, 
h could only be found through a large scale model test, was needed to fi- 
se the design. After considerable planning, it was decided to conduct a 
scale model test at the project site where the same materials could be used 
oth model and prototype, the construction equipment and control personnel 
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would be readily available, and alarge quantity of water and constant head cow 
be provided in the Oxbow reservoir upstream from the diversion cofferdan 

Because of the large expenditure involved in the construction and testing + 
the 1/2-scale model, as well as interference with upstream power generatic 
at Brownlee Dam during the test, it was impractical to repeat the 1/2-sca 
model test. In other words, the information which was desired from the lars 
scale model had to be obtained in one test. Therefore, preparatory tests | 
small scale models in the laboratory, and of large scale flume models in tl 
field were deemed necessary. 

From experience with river models which usually are of small scale al 
deal with small particle sizes, it is known that model results are, at best, on 
indicative. However, for large particle sizes and large scale models, mod 
test results are surprisingly accurate, as demonstrated at the Fort Randé 
closure.? 

It has also been established experimentally that the bed-load transport 
non-cohesive material hasa definite relationship to thetractive force. In ord 
to avoid the troublesome problem of simulating the particle size, shape, a 
compaction, it was decided touse the same materialand the same constructi 
control for both the large scale field model and the prototype and to establi 
the similitude relationship by using the tractive force, the only predomina 
variable, as the basis. It was also considered important to have a scale rat 
as largeas possible, to minimize not only the scale effect, but also the influen 
of critical tractive force. 


LABORATORY INVESTIGATIONS 


Laboratory investigations, including experimental tests and theoretical de 
ivations, were performed to determine (a) the mechanics of the washout pr 
cess, and (b) the similitude relationship for the rate of washout. : 

For convenience, the model tests were made in the downstream portion 
an existing model. A concrete pad and concrete blocks were placed to forn 
non-erodible base and a vertical wall was constructed on-the right (looki 
downstream) bank. 

The fuse plug models were constructed with axes angled 64° from the fa 
of the right wall. The typical cross sections of the models are shown in F 
2. These are for 1:20 and 1:40 scale models of the proposed Oxbow fuse pli 
The mechanical analysis of the materials usedin each model are shown in F 
3. 

No attempt was made to duplicate the density or other characteristics of t 
materials which would be used in the prototype. However, considerable effi 
was made in constructing the different 1:20 and 1:40 scale models to achik 
identical placement and compaction of like materials. 

In all, ten models were built and washed out as shown in Table 1. The ne 
mal core and thin core, referred to previously, represent the two types 
cores used, In the thincore section, the rockfill material was increased to 
maximum possible areafor the purpose of testing its influence on washout tir 
The dimensions of these cores are given in Fig. 2. 


3 “Dredge Fill Closure of Missouri River at Fort Randall? by Lorenz G. Straub, E 
ceedings, Minnesota Internatl. Hydr. Convention, September, 1953. 
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In the 1:20 model, the invert of the pilot channel was 3 in. lower than the 
st of the fuse plug, and the width (parallel tothe axis of the fuse plug model) 
5 6 in. In the 1:40 model these dimensions were 2 in. and 3 in., respectively. 
The operation of the tests was as follows: 


1. The forebay was filled to the level of the invert of the pilot channel, 

2. The forebay water level was then raised rapidly the prescribed amount 
proximately 0.4 in.) to achieve initial breach action, 

3. The forebay level was held more or less constant at this elevation dur- 
the washout process by continually increasing the inflow. 

4. A record was kept of the time at which the top and bottom of rows of 
ds spaced at 1-ft intervals along the model were washed away. 

5. Photographs were taken and the time recorded for each photograph, 

6. A movie was taken of Test Nos. 3 and 5. 


Embankment 
Material 


Rockfill 


Dimensions of Model in Inches 


FIG. 2,—DIMENSIONS OF LABORATORY MODELS 


Test Results. 

Yescription of Washout Process.—The washout process as observed in these 
S consists essentially of two phases, (1) the breaching of the pilot chan- 
and (2) the eroding of embankment face. The breaching action starts with 
nitial slip at the downstream end of the pilot channel (Fig. 4). This slip- 
= rapidly increases until the clay core is exposed. The flow then cascades 
r the clay core and begins to erode the material downstream nearly to the 
2 of the model. At approximately this time the clay core, which is left 
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PARTIGLE SIZE, IN MILLIMETERS 


FIG. 3.—MATERIAL USED IN LABORATORY MODELS 


FIG, 4.—INITIAL SLIP IN PILOT CHANNEL 
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anding aS a cantilever by the removal of the supporting downstream material, 
ils under the water loadand the breaching action is complete. The flow then 
gins to attack the base of the embankment, carrying away the material as it 
ides into the flow and producing a lateral advance of the eroding face along 
e length of the fuse plug. When a sufficiently large portion of the core is ex- 
sed, it breaks off. This process proceeds in a most regular manner grad- 
lly increasing the opening and washing out the fuse plug. 

The shape of the eroding face depends on the nature of the materials. For 
iformly graded material without fines, the slopes on the eroding face are 
adual and the material rolls continuously into the flow (Fig. 5), but for ma- 
rials that have some apparent cohesion or mechanical interlocking, the face 
ands nearly vertical (Fig. 6). With the latter materials, the face recedes by 
idden large slips of the material and sudden failure of portions of the clay 
re. 

Experimental Data.—The dataon the rate of washout for eachof the tentests 
e plotted as graphs in Figs. 7, 8, and 9, which show the eroded length of the 


TABLE 1 


Code Number of Test 


1:20 Seale 1:40 Scale 


Composition of Model 


arse sand - normal core 
a gravel - normal core 
ushed rock - normal core 
lter sand - normal core 
ushed rock - thin core 

a gravel - thin core 
ushed rock - normal core 
(Repeat of Test 3) 


(yu) tesco 


sdelas afunctionof time. In Fig. 7, Ris the rate of recession of the eroding 
se, in feet per hour; and Ry denotes the ratio of the rate of recession for the 
0 scale to the 1:40 scale model. The average values of the top and bottom 
shout times were used. Rows affected by breaching action were not included. 
In the interest of determining laboratory test accuracy, the 1:20 crushed rock 
del test was repeated. The results of the two tests were remarkably close. 
e differences are well within the expected accuracy of measurement of this 
ye of test (Fig. 10). 

Analysis of Experimental Results.—Analysis of the laboratory washout tests 
sults indicates the following: 

1, The washout rate is afunctionof the grain size of the material; the larger 
» size, the slower the rate of washout. However, the effect isnot pronounced. 
2. The washout rate for a fuse plug is between 1.05 and 1.20 times faster 


in the corresponding half-scale model. 
3. Increasing the volume of rockfill (reducing the thickness of the clay 


re) decreases the washout rate slightly. 

Theoretical Analyses.—After the pilot channel has been breached, the main 
ibankment washes out ina regular manner, as described previously. The 
lowing treatment of the mechanics of washout has been made to establish 
. similitude relationship. 
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FIG, 5,—SLOPING FACE WITH UNIFORMLY GRADED MATERIAL 


FIG, 6.—-VERTICAL FACE WITH MATERIALS HAVING APPARENT COHESION 
OR MECHANICAL INTERLOCKING (TEST 3) " 
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Referring to the definition sketch (Fig. 11) the rate of recession of the 
sroding face of the main embankment will be 


P|.Q 
-_— 
— 
— 


n which Ris the rate of recession of the eroding face (L/T), q denotes the rate 
wr erosionof material (L3/ T), and Arefers to the cross-sectional area of fuse 
lug (L2), 


EROSION OF FUSE PLUG, IN MODEL FEET 
EROSION OF FUSE PLUG,IN MODEL FEET 


TIME IN MINUTES TIME IN MINUTES 
COARSE SAND PEA GRAVEL 
LEGEND 


—O—O- =1:20 SCALE 
—X—X— =1:40 SCALE 


EROSION OF FUSE PLUG, IN MODEL FEET 


0 | 2 3 “4 
TIME IN MINUTES 
CRUSHED ROCK 


FIG, 7,—RELATIONSHIP BETWEEN 1:20 AND 1:40 MODELS 
OF IDENTICAL MATERIAL 
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The rate of erosion of material can be obtained from Du Boys equation 


in which G denotes the rate of erosion of material in weight per unit time/unit 
width of eroding face, 7 isthe sediment transportation characteristic (Fig. 12) 


Normal Core 7 
Test No. 3 


Thin Core, 
Test No. 8 


ieee core washout 


Steeper slope of normal 
core washout rates 

| indicates R is inversely 
proportional to the rock- 
fill area, 


EROSION OF FUSE PLUG IN MODEL-FEET 


=a 
PEA GRAVEL CRUSHED ROCK 


7 ee 


O I 2 3 4 5 
TIME IN MINUTES 


FIG. 8.—COMPARISON OF WASHOUT RATES FOR DIFFERENT 
ROCKFILL ZONE AREAS 


for riverbed material, t describes the shear or tractive force on the eroding 
material, 7; refers to the critical tractive force required to move material, 
(Fig. 12), and y is the unit weight of water. Now 
PhS SE 
== T(T-Te Sa dees ss RE es as te 
Bin at genet) 
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in which yg is the dry unit weight of material, and s denotes the width of the 
eroding face. 

t The shear or tractive force, Tt, for uniform flow in open channels is given 
el 


: a ce Tee ae ae 


Bottom washed out first, then 
top portion fell in slices. 


peaaae 


ce 
| 
z 
ma Test No. 7 
iS Fine Filter Sand 


O | 2 3 4 5 
TIME IN MINUTES 


EROSION OF FUSE PLUG, IN MODEL- FEET 


FIG, 9.—-WASHOUT CHARACTERISTICS FOR MATERIAL 
WITH APPARENT COHESION 


1 which y is the depth of flow, and S denotes the slope of energy grade line (or 


ater surface for uniform flow). 
For flow around the eroding face of the fuse plug, however, a more general 


Juation is required: 


T. EC 
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in which C¢ denotes the coefficient of drag, and p is the mass density of fluid. 
For uniform flow it can be shown that 


EROSION OF FUSE PLUG, IN MODEL-FEET 


| 
0 | ve 3 4 5 
TIME IN MINUTES 


fae e TEST NO.1I0, 1:20 SCALE, CRUSHED ROCK 
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FIG. 10,-COMPARISON OF WASHOUT TESTS FOR IDENTICAL MODELS 


maximum velocity available on the face, V=V2gh, the maximum ae 
force is 


j 
Ti fm Hes Ss ee (7) 
q 
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Jombining Eqs. 1, 3, and 7 
vy rr b[Cryn-ro]s 


R=, 7 (8) 


Ys y2 K H2 Mee ear Re RES! bee! 6 


n which K is the constant of proportionality inthe formula for the cross-sec- 
ional area of the fuse plug 


Bie eH eta makes oe See ae SESS (9) 
enerally speaking, s is proportional to h, that is 
Se ae Gee Ra (10) 


1 which a is a constant. 


Now, Cf is proportional to f, and f, in turn, is proportional to the 1/3 power 
f the relative roughness of the eroding face. This last relationship follows 


FLOW 


WATER SURFACE 


PATH OF ERODED MATERIAL 


PLAN DURING WASHOUT 


W.S. sn, Mae 


S. DOWNSTREAM 


SECTION A-A 


FIG, 11.—DEFINITION SKETCH FOR THEORETICAL CONSIDERATIONS 


m the fact that Manning’s “n” is proportional to the 1/6 power of the grain 
meter for flows over unriffled sediment beds and from the fact that f is pro- 
tional to n2. It follows then that 


d 
G = ¥(8 Peceanaere swig nniys Senversia (Ld) 
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in which d denotes the mean diameter of the material. Inserting Eqs. 10 an 
11 into Eq. 8 yields 


ah wy Eo ae (g) 7 : re ose (12 


Finally, for conditions in which 7, is much smaller than 7 (that is, small non. 
cohesive material in large fuse plugs) 


2 
_abk .2/3,1/3/h 

R= SEES g2/Sell3|o\ oye Soa eee 
¥g K H 


Eq. 13 provides a basis for determining the similitude relationship for rate o 
recession between model and prototype. Let the subscripts r, p, and m, refe: 
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to the words “ratio,” “prototype,” and “model,” respectively. Then, for ge 
metrically similar models using identical materials (which must be relative 


small in both model and prototype) it follows that : 
1/3 
y= ge =( 2] =, 1/3 
Rm \4m 


as all other ratios are unity. (Ly is the scale ratio of the model. ) 
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For two fuse plugs of identical materials with scale ratio of 1:2, the larger 
1€ will therefore erode 3V2 or 1.26 times faster. This relationshipis roughly 
srified by the laboratory tests described previously for which the 1:20 model 
roded 1.05 to 1.20 times faster than the corresponding 1:40 model. Although 
ie experiments indicate a rate of erosion for the larger model that is some- 
hat slower than the theory indicates, it should be remembered that neither 
e theory nor the experiments are precise and, at least froma practical stand- 
dint, this agreement between theory and experiment is acceptable. 

Summary.—From the laboratory tests and theoretical analyses on the me- 
ianics of washout of fuse plugs, the following has been found: 


1. The mechanics of washout can be explained on the basis of the laws for 
>diment transport. 

2. The rate of recession of the eroding face varies approximately with the 
‘3 power of the scale ratio for fuse plugs constructed withidentical materials 
1on-cohesive) and methods. 


LARGE SCALE FIELD MODEL TESTS 


General Considerations .—The tests of small scale models inthe laboratory 
d supplied the preparatory information necessary to correlate the results 
the large scale model test to prototype performance. In order to finalize 
e design of the Oxbow Idaho spillway, the construction of the 1/2-scale field 
odel was commenced in June 1959 with the test scheduled for late July. 
‘Purposes of Large Scale Field Model Tests.—The main purposes of the large 
ale model test were to verify the design assumptions, and to seek additional 
formation for the final design of the fuse plug. 

The assumptions to be verified were as follows: 


1. Mechanism of the breaching action at the pilot channel; 
2. Depth of water required to effect a successful breach; 
3. Completion of breach before lateral erosion; and 

4, Washout process for lateral erosion. 


The following additional information was needed: 


1. Time function of the breach process; 

2. Time function of the lateral washout process; 

3. Effect of the core on rate of washout; 

4. Completeness of washout; and 

5. Anything unforeseen which could cause malfunctioning of the fuse piug. 


Model Scale and Similarity.—The selection of a 1:2 ratio for the large scale 
ld model was prompted by a number of considerations. The primary factor 
s the existence at the project of a readily adaptable site where a 1/2-scale 
del could be constructed with a suitable length-to-height ratio. 

The predominant material in a fuse plug is the rockfill that occupies more 
n 75% of the volume (in the design shown), has relatively much larger size 
rticles than the filter, and, therefore, controls the rate of washout. As de- 
‘mined by the laboratory investigations, the similitude relationships between 
prototype and the 1:2 scale model withidentical materials and construction 
thods would be the same as between the 1:20 and 1:40 models. 
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Because the impervious core material was to have the same inplace density 
and moisture content for both the model and the prototype and, thereby, the 
same unit strength, the similitude relationship of the impervious core had to 
be determined according to its failure action. In this design, the core is ex- 
pected to break by cantilever load, both in breaching and lateral washout. 
Therefore, the core thickness was proportioned according to 


in which t, is the thickness ratio, and Ly denotes the scale ratio. In order not 
to disturb the similarity of the rockfill, the reduced thickness of the core (more 
than the 1:2 ratio) was made up by the fine filter material, which has very little 
effect on washout rate. 

For testing the initial breach action, the prototype head and dimensions of 
the pilot channel were used because the materials and slopes are identical in 
both model and prototype. 

Flume Model of Clay Core.—After having studied the general breaching and 
washout actionon small scale models in the hydraulic laboratory, it was more 
or less certain that the initial breach would start at the downstream edge of 
the pilot channel with a very small water head. However, because the imper- 
vious core material, the fine filter, and even the rockfill material may behave 
differently after being compacted to maximum density, it was proposed to test 
the breaching action in a wood flume at the site before proceeding with the 
large scale washout model test. The specific purpose of this model was as fol- 
lows: 


1, To observe closely the erosion phenomena of the downstream rockfill 
and filter material; j 
2. To verify the break-up action of the impervious core. 


To avoid all complications, the model in the flume was to have the same scale 
as that of the large scale field model, with the same materials and compaction. 
A wooden flume was constructed at the project site which was able to ac- 
commodate a 6 ft high crest section of the pilot channel portion of the 1/2- 
scale fuse plug. 
Two flume tests were performed, from which sufficient information was 
obtained for the construction of the large scale model. 


Test No. 1.—The flume model had a section as shown in Fig. 13. The im- 
pervious core was supported on processed fine filter material on the down- 
stream side, and covered by the same material on the upstream side. Flume 
floor and walls are of 1 in.-by-6 in, tongue and groove with polyethylene sheet 
lining. The mechanical analyses are shown in Fig. 14. : 

The model was constructed so asto obtain the same compaction in the cor 
and filter zones as in the prototype. After turning on the water and removin 
the downstream bulkhead, the downstream processed fine filter material ero- 
ded only to the extent of its contact with the lower nappe falling over the im- 
pervious core. Due to the apparent cohesion of the compacted, well-gradec¢ 
material, the processed fine filter material stood almost vertical for the en- 
tire height of the model. 3 

After mechanically removing the processed fine filter material, the im: 
pervious core was tested to determine the point of failure under cantileve: 
action. (The bond between the core and the walls of the flume was carefully 
broken at the top portion of the model.) The first break of the core occurre¢ 
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out 18 in. from the top with less than a foot of water over the top. The sec- 
1 break occurred about 4.5 ft from the top. 


Test No, 2.—The second flume model had a section as shown on Fig. 13. 
© impervious core was identical to that of Test No. 1, but the processed fine 
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SECTION FLUME MODEL NO.2 
FIG, 13.—LARGE SCALE FLUME MODEL TEST 


r material was replaced with a downstream rockfill zone of concrete ag- 
ate and a filter zone composed of a thin layer of concrete sand over pit- 
‘ine gravel. Thesubstitutes were of moreor less uniform materials. Their 
ations are shown on Fig. 14, | . 
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[he materials were placed in the same manner as for Test No. 1, but with 
visions to prevent bond and wedge action between the impervious core and 
walls of the flume. 

After the end bulkhead was removedand the test started, the material down- 
2am of the core was gradually washed out from the top. As the water hit 
downstream rockfill material and washed it out of the flume, the rockfill 
filter material beneath the nappe gradually receded by surface sliding. The 
t break of the unsupported core occurred about 18 in. from the top. The 


hout cycle was then repeated. The second break of the core occurred at a 
ht about 2.5 ft below the first. 


In the basis of the flume model tests, it was concluded that the impervious 
> of the large scale model was Satisfactory with respect to thickness and 
ungement, 

‘hese tests also indicated that the compacted, well-graded riverbed sand 
gravel, because of its apparent cohesion, would not be Satisfactory for use 
ie pilot channel section. It was demonstrated that clean, uniform material 


necessary for the downstream filter and rockfill to effect a rapid and suc- 
ful breach. 


lalf Scale Field Model Tests. 

escription of Half Scale Model.—The location chosen for the construction 
€ model was a natural rock saddle at the left abutment of the Oxbow cof- 
am. Here a large water supply was available from the natural river flow 
the Oxbow reservoir above the cofferdam could be used as a regulating 
to keep the water surface almost constant throughout the test. The pool 
een the cofferdam and the main dam could store all the water passing the 
21 without any flow over the Oregon spillway. 

he first step in preparing the model site was to remove part of the exist- 
offerdam, and expose the rock foundation. Because the rock surface was 
- irregular, a concrete sill was poured for the cutoff under the impervious 
- The rock onthe right side of the concrete sill sloped upto form anatural 
nent. However, on the left side, the rock was lower and very irregular, so 
‘tical steel sheet pile wall was constructed to form the abutment for the 
21, The piling was set in a trench excavated to rock and backfilled with 
rete. The top of the piling was securely anchored to logs buried well back 
e cofferdam. Hand-tamped clay was placed in back of the piling. 

he pilot channel in the model was adjacent to the sheet pile abutment. In 
rototype, the pilot channel will be adjacent to a concrete-faced rock abut~ 
, €xcavated as steeply as possible. No difference in breaching action is 
een due to the probable minor difference in abutment slopes. 

compacted clay foundation pad, 1 ft to 3 ft thick, was placed under the en- 
nodel to the level of the concrete sill. 

ie length of the base of the fuse-plug model, measured along the concrete 
was 51.6 ft. The crest length was 75 ft, and the height was 13.5 ft mini- 
except at the pilot channel. The crest was 1-ft higher in elevation at the 
abutment than it was adjacent to the pilot channel. The height at the pilot 
el was 11.5 ft. Typical sections of the model as constructed are shown 
y. 15, and a picture of the complete model is shown on Fig. 16. 

the lower 5 ft of the model, the rockfill materials were dumped directly 
ce by trucks, and spread by bulldozer. For the remainder of the model, 
aterials were stockpiled near the site, transferred to the embankment 
imshell as required, and spread by hand. 
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|Compaction of the riverbed sand and 
S accomplished by vibrating roller. 
ed tothese materials. The concrete 
iced. 


The clay core was spread in 3-in. lifts, sprinkled as required, and com- 
sted by pneumatic hand tampers. In order to get satisfactory bond between 
S it was necessary to roughen the contact surface and moisten it before 
cing the next lift. In general, compaction was between 85% and 90% of peak 
dified American Association of State Highway Officials (AASHO) density, 
ithe moisture content was within 1-1 /2% of optimum, 

Material Selection and Laboratory Tests.—The model was constructed of 
Same materials economically available for use in the prototype. The model 
3 designed as a composite section, using several rockfill materials, in order 
t washout rates for the various materials could be established and thereby 
ure selection of the most advantageous materials for the prototype. 


gravel, and the processed fine filter, 
Liberal quantities of water were ap- 
aggregates and fine gravel were heavily 


FIG, 16.—SANDBAGS JUST CLEAR OF PILOT CHANNEL 


b: 


he flume model tests had indicated that pit-run riverbed sand and gravel, 
9 slaced wet and compacted, had considerable apparent cohesion. Because 
roperty might slow the erosion process, it was decided to use this ma- 
1 where it would be the last to wash out; that is, at the end opposite the 
channel, a 
ithe pilct channel it was necessary to have an easily moved, cohesionless 
rial. Concrete aggregate passing a 3/4-in. screen and retained on.a No. 
2en was sure to meet these requirements and was therefore used. Be- 
there was an excess of the larger sizes of concrete aggregates available, 
3 decided to test a section using rock between the 3/4-in. and 1-1/2-in, 
_ This section, 20 ft long, was placed adjacent to the pilot channel sec- 
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The processed fine filter used in the main dam showed the same propert 
of apparent cohesion as the riverbed sand and gravel. Therefore, this filte: 
material was restricted to use only inthose sections using the sand and grav 
el. For the sections using concrete aggregates, a pit-run fine gravel with : 
minimum thickness of concrete sand onthe face in contact with the imperviou 
core, was used for the filter. The fine gravel showed no cohesion under am 
conditions. The impervious core was constructed of the same material as th 
main dam impervious core. The upstream and downstream slopes were face 
with the 3/4-in. to 1-1/2-in. concrete aggregate. 

For the results of standard laboratory tests performed on the various em 
bankment materials, see Table 2, and mechanical analysis curves, Fig. 14. 

In addition tothe standard laboratory tests referred to previously, a serie; 
of special flexural tests were performed on thecore material. The tests wer 


TABLE 2,—LABORATORY TESTS ON FUSE PLUG MATERIALS 


Shear Strength _ 


Specific] Maximum/Optimum|Coefficient 


Material 


(1) 


Clay core 
Fine Filter 


Riverbed Sand and Gravel 
Processed Concrete Sand 
Pit-run Fine Gravel 
Concrete Aggregate 

No, 4 to 3/4 in, 
Concrete Aggregate 
3/4 in, to 1-1/2 in, 


Se ee eae eee ed 


4 Maximum density determined by: (1) Modified AASHO method for clay, (2) fiel 
density test after compaction for fine filter, and (3) dry, rodded method for all othe 
materials, : 


required to quantitatively evaluate the maximum cantilever projection of th 
core that could occur in the pilot channel section during removal of the down 
stream supporting rockfill. This value, in turn, permitted a rough computa 
tion of the time required to start initial breakup of the core, and the effecti 
breach of the section. Flexural strength values were desirable for both a 
as-placed material and the saturated material to get acomparison between tk 


be required to wash out. 


The test specimens were prepared by compacting the soil to approxima te 
the same density, and at about the same moisture content, as the materi 
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aced in the model. The specimens were 2 in. square in section, and 10 in. 
ng. The soil was placed and compacted in about five layers. 

The test consisted of Supporting the specimen at its ends, and applying a 
iformly increasing load at the third points until rupture occurred. The load 
as applied normal to the directionof compaction, to reduce the effects of lam- 
ations in the test specimen, 

The test results showed a modulus of rupture between 3000 psf and 4000 psf 
r the material as placed. No conclusive results were obtained on the effects 

saturation. However, the tensile strength was apparently significantly re- 
ced due to prolonged submergence. 

Measuring Devices. 

Embedded Washout Rate Indicators:—To get informationon the rate at which 
> interior portions of the fuse plug were being eroded, a series of washout 
te indicators were installed during the course of construction. These were 
mply electrical switches which were held closed by the earth pressure with- 
the embankment. As the erosion progressed the earth pressure would be 
lieved, the circuit broken, and a‘light would go out at the observer’s station. 

The indicators were installed ina vertical plane along the center of the 
est on three vertical lines as shown on Fig. 17. In general, they were in- 
ulled at 2-ft increments of elevation. The wires from each indicator went 
rtically down to the clay foundation, upstream to the concrete sill and then 
- over the Idaho abutment. This routing prevented damage to the wiring be- 
‘€ operation of the respective indicators. 

Grid Lines on Downstream Face:—A 5-ft grid was painted on the downstream 
e of the fuse plug model, as shown in Fig. 16. Three independent observers 
-orded the elapsed time between the start of the test and the removal of the 
d intersections, as an aid in establishing the washout rates for the various 
terials. 

For the final determination of washout rates, a careful study of the motion 
ture, taken at the time of the large scale model washout, was made. The 
d lines were the dimensional control, and repeated stop watch measure- 
nts established the time. 
Description of Washout. 

Provisions for Starting Test:—A 10-ft long section of the crest was left 2 ft 
er in elevation than the remainder to form the pilot channel, through which 

water would flow and start the breaching action. Close control of the res- 
oir water surface elevation was required, but difficult to achieve. In order 

‘et the proper head for the test, without danger of premature washout, sand- 
S were employed in the pilot channel. 
four layers of sandbags were used to close the pilot channel to the model 
st at El. 1748.8. The sandbags were arranged ona wire rope sling at- 
ied toa truck crane. After the reservoir stabilized at the desired eleva- 
, the sandbags could be lifted quickly to start the test. 

‘inal Test:—In the morning of July 22, 1959, all preparations for the final 

had been completed. The reservoir water surface was brought up to the 
er level, one foot above the invert of the pilot channel. At 30 sec after 

1 a.m. the initial breach began with the lifting of the sandbags clear of the 

t channel (time consumed in clearing sand bags, 9 sec), as shown in Fig. 


iS Soonas the sandbags were lifted, the water gushed through the pilot chan- 
and, flowing down the slope, commenced the erosion of downstream rock- 
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fill starting at the top of the slope. Within 35 sec after the startof the test the 
erosion had progressed tosuch an extent that the upper portion of the core was 
acting as a weir, with the nappe clear of the material below, Fig. 18 (time 
11:12+). The first break of the core occurred 1 min and 7 sec after the start 
of the test. Fig. 19 (time 11:13+) shows the moment at which the second break 
of the core had just occurred. The time, as calibrated from the motion picture 
record was 1 min and 27 sec after the start of the test. The breach action was 
considered completed at this point as no further change occurred in surface 
profile of the water rushing through the gap or in the color of the water. 


Existing 
Cofferdam 
EMBEDDED WASH-OUT INDICATOR EL. |7524 
Line AN Line BN Line C Steel Sheet 


Piling Abutment 


16:50 14:10 


4164 EL 1745.35} 
#386 3:00 | 
; 

505 1__|350 Figo “t-'740.37 | 
+36 | 
630 J 
6:30) J 


BASE EL. 1735.3 


DOWNSTREAM ELEVATION OF MODEL 
OBSERVED WASHOUT TIME 


The numbers on the signs along the crest designate |O foot 
stations from the Sheet Pile Abutment. 

Time is in minutes and seconds from start of wash-out. 

yminzsec wash-out of grid intersections, 

qmin:ses wash-out of embedded wash-out indicators 


Ist break of core 1:07 
2nd break of core |:27 
(completion of breach) 


Approximate limit for 


concrete aggregate section 
NJ 
Approx. 


VS 


§ 
y 
: 
x 
4 
: 
; 
4 
x 


(Time |Om. 30s. avg) Time 2m. 3is, 


DOWNSTREAM ELEVATION OF MODEL 
FIG. 17,—RECORD OF WASHOUT RATE 


The lateral erosion of the model began during the initial breach period : 
the rockfill material along the side slidinto the water rushing through the g: 
and the core was left overhanging along the side. As the uniformly-gr 
rockfill and filter slid out from underneath, the core was left cantilevered a 
sections would break off at intervals. The surface of the uniformly-grad 


4 


_ FIG, 19, 


WASHOUT 


—AT MOMENT OF SECOND BREAK OF CORE. 
-BREACH ACTION COMPLETE 
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material sliding intothe water assumed a slope slightly steeper than the angle 
of repose of the material. This phase of the washout tookapproximately 6 min. 
Fig. 20 (time 11:17) shows the characteristic slope of the uniformly-graded 
rockfill and filter under full face lateral erosion. The change fromthe uniform 
to the well-graded material occurs at Station 35 and the vertical break can be 
seen at the downstream face above El. 1795, the upper white line. 

The washout of the well-graded material is characteristically different from 
that of the uniform material. Due to apparent cohesion and a small amount of 
fine material in the well-graded riverbed sand and gravel, the end of the par- 
tially washed out fuse plug model assumed anearly vertical attitude. The ero- 
sion of this portion of the fuse plug model proceeded as the churning water 
undermined the vertical face. Atintervals the undermined portions of the down- 
stream rockfill broke off and fell into the water and the material was carried 
away. As the core became exposed and was left unsupported, sections would 
break off. Fig. 21 (time 11:22) shows the vertical face of the well-graded ma- 
terial and the churning, undermining action of the water. 

As soon as the test was over, the spillway gates at Brownlee were closed, 
and the powerhouse release was orderedcut down as much as possible so that 
the foundation-could be examined as early as possible. 

At 12:30 p.m., the powerhouse release was reduced to 6000 cfs, and kept 
there until 2:30 p.m. The Oxbow reservoir receded rapidly. At 3:15 p.m. the 
water surface downstream of the fuse plug was at the same elevation as the 
Oxbow reservoir, and the water started to flow into the reservoir. At this time 
the reservoir water surface, as well as the tailwater pool, was at El. 1738.5, 
4 ft below the temporary crest of the Oregon spillway. 

Analysis of Test Results.—The large scale field model test has clearly 
demonstrated the details of the washout action. Information which had pre- 
viously been determined from the laboratory and flume model tests was con- 
firmed, and some new information was revealed which could not be found from 
the other model tests due to their limitation in scale and size. 

The continuous motion picture record’has been found most valuable, because 
it enabled the engineer to study the washout action closely.and repeatedly. 

After the Oxbow reservoir and the pond downstream of the cofferdam had 
receded, the model foundation was examined. 

All information obtained, including the river gage records at Oxbow and Hells ? 
Canyon, have been carefully analyzed, and the results are presentedas follows: 


Breach Action:—The breach action was entirely as expected. It showed 
close similarity tothe small scale model. The core broke up exactly as in the 
flume model, first, a small break, then a large break, and then repetition of the 
breach action all the way to the base. There was no significant seepage into 
the rockfill at the downstream edge of the pilot channel before the initial breach. 

Lateral Washout:--The uniform material washout process was exactly as’ 
estimated. The kinetic angle of repose was about 45° for the uniform material, 
For the material with apparent cohesion and a small amount of fines in it, the 
washout face was vertical, and the washout process was a three-step action; 
undermining the base of the vertical face, vertical slide of the rockfill and filter 
material, and the breakof the impervious core. The embedded indicator record 
also indicated vertical face erosion under water. 7 

Rate of Washout:—Basedupon direct observations (three observers) and the 
motion picture record, the washout time function is shown in Fig. 13, and the 


¥ 
i 


FIG, 21,—VERTICAL FACE OF WELL-GRADED MATERIAL 
(COMPARE WITH FIG. 6) 
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results summarized as follows: 


Breach (Complete) 1 min and 27 sec 
Rate of washout, uniform material 5.6 ft per min 
Rate of washout, well-graded material 1.4 ft per min 


Hold-up time of the core was insignificant. 

Order and Completeness:—The entire washout proceededin anorderly man- 
ner, Breach was completed before lateral erosion. There was no evidence of 
any portion of the plug remaining beneath the water during the washout and 
Fig. 22 shows the completeness of the washout after the foundation was ex- 
posed next morning. The reservoir water surface at the start of test was El. 
1747.80. The reservoir water surface at the end of test was El. 1747.98. The 
maximum water surface during test was El. 1747.98. 


FIG, 22,—LOOKING TOWARD LEFT ABUTMENT OF STEEL SHEET PILING 


4 
E 
Headwater and Tailwater:—From reservoir water gage records, it was foune 
that during the entire test, the headwater varied only 0.18 ft. The breach he: 
was 1 ft. The slope of the ground downstream of the model was approximatel 
10%. There was no settling of the material in the escape channel. The di 
downstream of the model was constructed to study its effect on the washou 
but no noticeable effect was found. 4 
Summary.—The large scale field model test has conclusively demonstrate 


the following: 


1. The structure showed the ability to withstand a water load with a resel 
voir surface at or above the pilot channel invert for 4.5 hr, and also a rap: 
drawdown actually more severe than expected for the prototype. 

2. The breaching action was initiated with a water head one foot above tl 
pilot channel, and completed to the base of the model in less than two minute 
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uniform, cohesionless material is essential for the downstream filter and 
/ckfill zones at the pilot channel to achieve a rapid and successful breach. 

3. The lateral washout process for thetwo different types of materials was 
early demonstrated. The rate of washout for the uniform, cohesionless ma- 
rial was found to be 5.6 ft per minute, and for the well-graded material with 
parent cohesion 1.4 ft per min. Through proper proportioning of the lengths 
volving the two kinds of materials, a suitable value for the total washout time 
the prototype fuse plug can be achieved. 

4. Theimpervious core design is satisfactory, because the “hold-up” action 
the core was found to be insignificant. 

5. The model fuse plug was completely washed out to the concrete founda- 
m sill in an orderly fashion, without unforeseen developments. 


CONCLUSIONS 


The small scale washout tests in the laboratory have established a definite 
tern of similarity between a 1:20 and 1:40 scale model, composed of the 
ne materials and constructedin the same manner. The 1/2-scale field model 
t has also shown similar washout mechanics. 

The ratio of washout rates from the small scale model tests, for identical 
terials and construction, and scale ratio of 1:2, checks very closely with 
theoretically derived value. Therefore, it can be logically concluded that 
actual washout rate of the prototype can be closely approximated by apply- 
this ratio of washout rates to the results fromthe 1/2-scale field model test. 
The design assumptions of the fuse plug have been verified, in all respects, 
the actual large scale field test. Based upon the performance of this test, 
prototype fuse plug design, including the selection of materials, will be en- 
ly satisfactory in effecting a successful breach, and a complete washout, in 
orderly manner, as planned. ; 

By using two types of rockfilland filter materials, it is considered definitely 
sible and practical to construct a fuse plug for suitable washout rates along 
erent portions and an appropriate total duration of the washout. 
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SEEPAGE THROUGH LAYERED ANISOTROPIC POROUS MEDIA2 


By David K. Todd,! M. ASCE, and Jacob Bear2 


SYNOPSIS 


The paper summarizes an investigation of seepage from leveed rivers into 
-lying adjoining agricultural lands. Only the flow rate and distribution of 
page as a function of the surface and subsurface boundary conditions were 
lied. Flow nets in idealized cross sections were determined by an electric 
logy model. Boundaries represented channel bed and bank, a levee base, a 
=r impermeable stratum, and the water table in adjoining drained agricul- 
U land. Variables studied included two subsurface layers of different per- 
bilities and anisotropies, layer arrangements and thicknesses, channel depth 
width, levee base width,. and water table slope. 


INTRODUCTION 


eepage from rivers and natural channels is a serious problem in many lo- 
ons. Significant contributions to the literature on this and other types of 
age have been made by several engineers (1), (2), (4), (5), (6), (11), 
, (17), (21), (26).3 When channel levels exceed the adjacent ground sur- 
elevation, water moves through and under the confining levees into adja- 
lands. If drainage facilities are inadequate, the soil becomes saturated, 


ote.—Discussion open until October 1, 1961. To extend the closing date one month, 
itten request must be filed with the Executive Secretary, ASCE. This paper is part 
e copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
‘Society of Civil Engineers, Vol. 87, No. HY 3, May, 1961. 

Presented at the March 1960 ASCE Convention in New Orleans, La. 

Assoc, Prof, of Civ, Engrg., Univ. of California, Berkeley, Calif, 

Formerly Graduate Research Engr., Univ. of California, Berkeley, Calif.; present- 
ccturer, Technion, Israel Inst. of Tech., Haifa, Israel. : 
Numerals in parentheses refer to corresponding items in the Appendix Bibliography, 
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water ponds on the surface, and crops are damaged. The Sacramento Valle 
of California is a typical example of an area which suffers extensive seepag 
damages. Losses in the period 1937-1953 (at 1953 prices) were reported a 
$1,834,000 per year, or $102 per acre affected (20). The problem of seepag 
in this area prompted the present investigation. 

The purpose of the investigation was to study the physical picture of flow ¢ 
water from stream channels through subsurface strata to the adjacent groun 
surface. Natural conditions governing these flows include geology and soils 
topography, ground water, rainfall, runoff, and channel stage. The range ¢ 
conditions for each of these variables, together with their possible combina 
tions, makes it apparent that not all situations could be studied. Instead, ke 
variables were selected for study so that results could be applied not only t 
the Sacramento Valley, but also to other areas with similar conditions. Th 
present work deals, therefore, with general aspects of seepage from rivers 

Most seepage studies in the Sacramento Valley have been concerned wit 
field conditions of the affected areas, periods and durations of flooding, an 
economic losses. Comprehensive surveys are Summarized in reports by th 
Joint Committee on Water Problems of the California Legislature (9) and k 
C. E. Plumb, M. ASCE, and others (20). Because the present study was a lak 
oratory investigation of seepage flow nets for steady flows, these reports an 
others (10), (16) furnished valuable information relating to construction ¢ 
idealized cross sections. Results from an electric analogy model were ana 
lyzed and presented in graphical form to facilitate interpretation and applica 
tion to various field conditions. In view of the number and complexity of cros 
sections to be investigated for quantities and distribution of seepage, mods 
investigations rather than analytic studies were necessary (application of ana 
lytic techniques to two simplified special cases is shown in Appendix I). Hy 
draulic models have often servedas a valuable tool to supplement analytic aj 
proaches (24, pp. 307-321). The electrical analogy was most suitable for tl 
present investigation. Analysis of thistype of modelandits application to se 
age studies have been presented by N. N. Pavlovsky (19), C. J. G. Vreende 
burgh and O. Stevens (22), (27), (28) and others (8), (12). A recent s 
mary of the theory of the analogy has been prepared by L. C. Malavard (14 

In most cases, subsurface strata through which seepage flows occur a1 
nonhomogeneous and anisotropic. Therefore, the main emphasis in the prese 
investigation was on seepage through layered anisotropic soils. : 


SEEPAGE CROSS SECTION 


Conditions along the Sacramento River servedas a basis for selecting cro; 
sections toinvestigate. The variety, irregularity, and wide diversity of layer 
alluvial deposits along the river precluded the possiblity of determining see} 
age through actual geologic sections without having extensive field data ava 
able. Instead, about fifty representative hypothetical cross sections were ca 
structed for model tests (25). q 

Almost all seepage occurs in the zone above the highest more or less ce 
tinuous clay layer, which acts as a lower impermeable boundary. In the pe 
meable zone various layers of deposits can be found. Based on types of ge 
logic deposits, layers could be reduced to those having relative permeabiliti 
of 50, 1, and 0. These layers are essentially horizontal, but may vary inthi 


ness and in their relative vertical positions. As most seepage occurs witl 
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3 than a mile of the river, the model was limited to a maximum distance of 
)0 ft inland fromthe levee. A maximum depth of 100 ft from the water table 
he impermeable clay layer was also selected. If Symmetry is assumed 
ut the center of the river, a section on only one side of the river need be 
sstigated. 
Other assumptions made for the model investigation were that the river 
mel is rectangular with a depth of 25 ft (a few cases of other depths were 
) studied) and a width of 600 ft. Levees were assumed to be Situated ad- 
ing the river channel, to be impermeable, and to have a base width of 100 
_ the elevation of the inland water table. An impermeable levee, of course, 
lies that all Seepage passes under the levee. Depending upon the mode of 
struction and materials employed, seepage through the levee may or may 
be a significant item. 
. final variable necessary in evaluating seepage is that of the anisotropy of 
neable zones. Most alluvial strata are deposited horizontally with a pro- 
iced tendency for agreater uniformity in grain size and depositional struc- 
horizontally than vertically. Commonly this results in a condition of ani- 
opic permeability where the horizontal permeability may greatly exceed 
of the vertical. Ratios of vertical to horizontal permeability of 10M and 
were selected to give a range of representative conditions. 
he position of the water table relative to river stage and to ground surface 
rmines relative seepage rates and the length of time before seepage will 
ar at ground surface. Seepage to a deep water table from a short duration 
[stage may not affect surface conditions, whereas at a location with a high 
r table it may have a marked surface effect. 
1 most irrigated agricultural land the water table is controlled to within 
ow limits by subsurface drainage facilities. The level is commonly main- 
d at 3 ft to 4 ft below ground surface. For steady flow conditions studied 
e electric analogy model, the water table was assumed to be stationary. 
might represent levels controlled either by subsurface drains or by pond- 
t ground surface. Both horizontal and sloping water tables were investi- 
1 


ield (9) and laboratory (23) investigations have confirmed that seepage 
flood flows is proportional toa duration-height factor. This study was 
ed to the steady case with emphasis on the distribution of seepage inland 
the river, Factors of river stage magnitude and duration are variables 
1 are necessary in applying these results to obtain quantitative seepage 
given situation and location. 

ised on these considerations the general cross section to be studied is 
n in Fig. 1. From Fig. 1 the channel depth, d, has a range of values of 
9 100 ft; the permeable zone depth, D, is a constant 100 ft; thickness of 
‘layer, u, ranges from 10 ft to 100 ft; thickness of lower layer, 1, has a 
of 0 ft to 90 ft; length of cross section from centerline of river, L, has 
ge of 1,000 ft to 3,300 ft; no range is assigned for the river stage above 
base, H; the range of channel half-width, B, is 100 ft to 300 ft; the levee 
width, b, is a constant 100 ft; distance landward from toe to levee, x, 
s from 600 ft to 2,900 ft; no ranges are assigned pa Siac appa of 
Ye r layer, Ky, or the lower layer, Ky; the ranges for both the anisotropy 
ee kik, and the lower layer, KyL/KxL, are 1/1 to 1/100; 
e ranges for the ratio of permeabilities to the two layers, Ky/Ky,. 

e channel depth d and the channel half-width B were constant for all but 
model tests. Significant variables are the thicknesses of the upper and 
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lower layers, u and 1, respectively. These are related by 
D2 41 100 ftyestevdie Setnnlasennelom (; 


The lengthof cross section L was varied as required to obtain a practical ve: 
tical scale to define the seepage flow pattern in each model test. Although @ 
river stage H determines the rate of subsurface flow, no values were assigni 
since it has no effect on the distribution of seepage. Similarly, the perme: 
bilities of the layers, Ky and Ky, were unassigned, but their relative valu 
expressed by their ratio Ky/Kz form a key variable of the investigation. T. 


IMPERMEABLE RIVER 


GROUND WATER TABLE LEVEE 


SURFACE 


UPPER 
LAYER 


FIG, 1.—DEFINITION SKETCH QF CROSS SECTION ASSUMED 
FOR MODEL ANALYSIS 


. 
anisotropies of the individual layers are expressed by the ratio of their : 
spective directional permeabilities, Kyu/Kxu and KyL/KxL- Based on - 
variables defined in Fig. 1 and on a study of field conditions in seepage ar 
of the Sacramento Valley, a series of model tests was outlined covering " 
ous combinations of the significant variables (25). 7 


THEORY OF THE ELECTRIC ANALOGY 


The electric analogy is a well-known device for studying the field of fl 
fluid through porous media (8), (18). The analogy is based on the simila 
between the differential equations which describe the flow of fluid through 
ous media, and those which govern the flow of electric current through co f 
ting materials (14). Methods for treating anisotropic layered media by tr: 
formation into equivalent isotropic fields by scale distortion have been 
cribed by M. Maasland (13). ~~ 
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Model scales and dimensions were based on 


2 
ev Syy 
= Roy saree Seco ic oo on (2a) 
ru 
X251, Kyz 
= rae Systema Mien sis Oemiemere Sacra tn 8 (2b) 
Y rL 


Xpu = Xpy = Xp 


ommon boundary between layers was made to coincide after transforma- 
An illustrative example of a typical cross section showing prototype and 
1 dimensions is shown in Fig. 2. 


DESCRIPTION OF THE MODEL 


1¢ electric analogy model for this study consisted of several components, 
ding the electrolytic tank, the electrical circuit, the conducting medium, 
1alogy boundaries, and the drawing table. A photograph of the assembly 
own in Fig. 3. 
} the conducting medium was an electrolyte, a large shallow lucite tank, 
dimensions 25 in. by 60 in. by 3 in., was constructed to hold the liquid. 
ottom of the tank was plane and was levelled to insure a uniform depth of 
‘olyte. In all tests the maximum length of the tank was utilized. It was 
that making a test of the largest possible size in the model minimized 
iracies and boundary errors. 
e electric circuit for the model is shown in Fig. 4. A 12.6v 60-cycle cur- 
jas connected to the electrolytic tank. The potential divider accurately 
d the voltage into 1,000 subdivisions so that readings to 0.1% were pos- 
Visual identification of the potential at the probe position in the tank was 
ed with a cathode ray oscilloscope connected to the probe and the poten- 
vider (Helipot). Whenthere was no potential difference between the probe 
e desired potential set on the potentiometer, no current flowed through 
cilloscope. 
2 resistance box provideda means for establishing a wide variety of fixed 
ials on electrodes in the tank. Most tests required only the limiting po- 
Sof 0% and 100%; but to represent sloping water tables, intermediate 
als with specified differences were necessary. The box contained forty 
ion resistors of 150 milliohms each, connected in series. In order to 
adjustments when intermediate potentials were being connected to the 
des in the tank, resistances in the resistance box were made relatively 
4s compared to the resistance of the electrolyte in the tank. In this way 
ry small amounts of current were consumed at the intermediate points 
main line. The deviation by variable current from a linear distribution 
ntials along the line was negligible. 


es 
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Various concentrations of a copper sulfate solution as well as tap wat 
were employed as the conducting medium for these tests. Liquid conducto 
are advantageous when a large number of tests with different boundaries a 
to be studied. Solutions could be poured into the model readily, and after ea 
test-could be drained into storage containers. Also, where apermeability rat 


L= 1000’ 


(a) PROTOTYPE 


(b) MODEL 


FIG, 2,—-ILLUSTRATION OF PROTOTYPE AND MODEL DIMENSIONS FOR 
ONE CROSS SECTION 7 


in two layers of 50:1 was under investigation, this difference was obtain 
selecting solutions witha corresponding conductivity ratio. Liquid depths 
the entire cross section were kept equal to about 2in. Withalternating c 

copper electrodes, and a copper salt in the tank, no problems of polariz 
were encountered at the electrode surfaces. ’ 
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FIG, 3,—VIEW OF ELECTRIC ANALOGY MODEL 
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FIG, 4.—ELECTRIC CIRCUIT FOR ELECTRIC ANALOGY MODEL 
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Three kinds of boundaries were required in the model: permeable and in 
permeable outer boundaries, and inner boundaries between layers of differe 
permeability. Permeable boundaries were made of copper strips 3 in. wit 
bent to fit the configuration of the boundary and carrying a given potential | 
means of a connection with the resistance box. For boundaries representil 
sloping water tables, a variable potential was required. Because the pressul 
is constant (atmospheric) on a water table, the potential is everywhere in prt 
portion to the elevation. Therefore, such a boundary was reproduced by a st 
ries of closely-spaced copper electrodes (Fig. 3), each connected to its co! 
responding potential in the resistance box. The potentials along the bounda1 
were expressed in percentages of the total head difference betweenthe chann 
water level and the lowest water table elevation. Impermeable boundaries we! 
nonconducting strips of lucite, 3 in. wide, or the lucite wall of the tank. 

The inner boundaries had to provide a waterproof seal between the two lit 
uids, and had to conduct electricity freely from one layer to the other, but ea 
point on a boundary had to maintain a unique potential value. The boundary wi 
constructed by modifying a procedure first suggested by Vreedenburgh al 
Stevens (22), (28). Lucite strips were notched and perforated at 3/16-in. ii 
tervals, and a short length of 1-mm diameter copper wire was inserted in eat 
hole and bent sothat the ends crossed in the notch above. All boundaries we: 
held in place by modeling clay at the bottom. 

Outer boundaries were interchanged during each test. Electrodes on pe: 
meable boundaries enabled equipotential lines of electricity and flow to | 
sketched, while placing electrodes on impermeable boundaries allowed strean 
lines to be located. Variable potential boundaries could be interchanged on 
after computing the redistribution of potentials along such boundaries. 

Movement of the probe was recorded on a drawing paper at a 1:1 scale | 
means of a pantograph connecting the probe and a pencil. By connecting poin 
of the same potential, each equipotential line was obtained. Sketching time f 
one set of lines of a cross section amounted to about 30 min. 

In general, with proper instrumentation, any degree of accuracy can be 
tained in an electric analogy model. As the present work was concerned w 
ground-water flow, basic data such as boundaries, permeabilities, and ani 
tropy ratios were only approximately known. A high degree of accuracy, ther 
fore, was not required in the idealized model cross sections. Results we 
estimated to be accurate to within 5%. 


THE FLOW NET 


The curves ¢(x,y) = constant (equipotential lines) and #(x,y) = cons’ 
(streamlines) form two orthogonal families of curves which together form 
flow net. Each two neighboring streamlines form a flow channel, in which 
amount of flow is equal to the difference between the two ~ -values of the 
responding streamlines. Usually a flow net is constructed withan integer m 
ber of flow channels, each transmitting the same amount of flow. Simila: 
the number of equipotential lines is chosen such that the drop of head betw 
adjacent ones will be the same throughout. If H is the total loss of head, 
denotes the number of equipotential drops AH, and Ng defines the num 
stream channels with a flow rate AQ per unit width in each, then by Dar 
law 

AH 


AQ = Bina «1 0.0068 4 
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which a is the streamline spacing and b is the equipotential spacing. The 
io a/b remains constant throughout any homogeneous isotropic region. 

For reasons of simplicity it was decided to draw the flow nets for the pres- 
work with Ne = Ng = 20. This gave a potential drop between two equipoten- 
lines equal to 5% of the total available head H, and a flow ineach flow chan- 
equal to 5% of the total flow from the river. The total flow per unit dis- 
ce along the river channel was given by 


an average of a/b from several measured “squares” was found sufficient 
nsure a reliable flow value. When two layers are present, the ratio a/b is 
erent in each zone, so that 


fo draw equipotential lines, the river bank and bed [ABCD in Fig. 2(b) ] 
e connected tothe 100% electrical potential, and the ground-water level and 
end of the model [ EFGH in Fig. 2(b)] were connected to the 0% electric 
ntial. The impervious base of the levee (ED), the impervious bottom (HI), 
the right edge of the model which is a line of Symmetry (IA) were made 
on-conductive lucite. From a theoretical standpoint the line FH should 
> been infinitely distant from the channel, but the error introduced by as- 
ing it to be about 1,000 ft to 2,000 ft from the channel is negligible. For 
e boundaries and a given Kyj/Ky, ratio, the equipotentials were drawn. 

‘or drawing streamlines, impermeable and fixed potential boundaries were 
rchanged. At the sametimethe conductivities of the two layers were inter- 
iged. 

‘he flow net drawn directly from the isotropic model is an orthogonal one. 
prototype itself is anisotropic; therefore, the actual flow net is not or- 
onal, When the drawing is reduced by distorted scales, that is, the hori- 
al reduction differs from the vertical, the resulting drawing is also a non- 
ogonal flow net (see Figs. 5 and 6). 

Imost all tests had an upper boundary formed by a fixed horizontal water 
}, the elevation of which was taken as the datum level (¢= 0). This ap- 
imated conditions existing in fields adjacent to a river where drains con- 
the water table. In this way the problem of an unknown phreatic surface 
avoided; a phreatic surface can be studied in an electric analogy only by 
and error. However, in some field cases of interest the controlled steady 
r table is not horizontal, but slopes away from the river with a gradient 
om about 1:500 to 1:1500. Two tests were performed with sloping water 
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tables to study these situations. The lowest point of the water table (at a dis- 
tance of 3,000 ft from the center of the river) was taken as the datum level 
(@=0), anda slope of 1:500 was assumed for the water table. An example 
giving prototype dimensions is shown in Fig. 7. 

The-equipotential lines for these tests were drawn by the electrical analogy, 
while the streamlines were drawn by constructing an orthogonal family of cur- 
ves to the equipotential lines according to the rules of graphical flow net con- 
struction (3). 


RESULTS 


Results of the model tests were summarized inthe form of plates and tables 
(25). Two plates, shown in Figs. 6 and 7, illustrate cross-sectional flow nets. 
Seepage quantities and distribution were summarized in tabular form. For 


ff VARIES FROM 0 TO 30% OF H 


948' 


Y+o% 


FIG, 7.—EXAMPLE OF CROSS SECTION WITH SLOPING WATER TABLE 


each cross section, a seepage factor was obtained from measured values of - 
in the upper layer, and from Eq. 6, ] 


om ts 


Ky H 


Therefore, the total seepage from one side of the channel for a given cro! 
section is obtained by multiplying the factor a/b by appropriate values of K 
and H. The fraction of total seepage emerging through the vertical side ¢ 
depth d was determined as wellas the distribution of Seepage at ground surfac 
as a function of distance inland from the levee. 

A dimensionless scaling factor a was computed for each cross section 
that the test data can be extended to other anisotropies of the upper and lowe 
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yers. The factor is related to the permeabilities and scales by 


’ = t Kxup GES) 
Y.y*a@x ru (ae) 52 he ee (9a) 
d 
t ’ Kx 1/2 
Me gh eeeX (ae Athen, Tie oes (9b) 


which X'yy,, X' yy, Y' rU, and Y'yy, denote the ratio of the figure length to 


> corresponding prototype distance and other terms are as previously de- 
ed. 


ANALYSIS OF RESULTS 


Seepage flow nets for the cross sections studied depend upon the variables: 
| channel depth, (2) channel width, (3) levee base width, (4) anisotropy, (5) 
‘angement and thickness of layers, and (6) water table slope. From the test 
sults an effort was made to evaluate the effect cf each factor on seepage, 


Effect of Channel Depth.—Discharges as ae are plotted against the depth 


io d/D in Fig. 8. In general, the discharge increases with channel depth, 
rate of rise being dependent upon the anisotropy ratio. For an increase in 
th ratio of 0.25 to 1.00, the seepage increases 8% for an anisotropy ratio of 
00; an increase of approximately 45% occurs for an increase in depth ratio 
m 0 to 0.25. For an anisotropy ratio of 1:10 the increase in discharge is 
for an increase in depth ratio from 0.25 to 1.00 and the same for an in- 
ase in depth ratio from 0 to 0.25. 
Analytic evaluations of channel depths d = 0 and d = D are discussed in Ap- 
dix I. 
ffect of Channel Width.—Most of the tests were performed with a half- 
mel width equal to 300 ft, so that a direct observation as to this effect for 
various anisotropy ratios and thicknesses of layers could not be made. 
ever, it was apparent that a relatively small portion of the total seepage 
irs through the channel bottom when the upper layer is a pervious one not 
nding appreciably below the channel bottom. 
n tests in which the bottom of the upper pervious layer extends to well be- 
the channel, the percentage of seepage through the bottom increases up to 
‘oximately 50%. The percentage depends upon the geometry of the cross 
ion and the anisotropy ratio. When the upper layer is semipervious, most 
€ seepage takes place through the bottom unless the thickness of this layer 
naller than the depth of the channel. 
he effect of channel width for a homogeneous anisotropic case was demon- 
led in tests for bottom widths of 300 ft, 200 ft, and 100 ft. The resulting 
Q 
Ky H 
1 differences were obtained inthe analytic treatment of a zero depth chan- 
Appendix I). 
fect of Levee Base Width.—Because all tests were performed with a levee 
width of 100 ft, no conclusions could be drawn as to the effect of this fac- 


age factors were = 1,08, 1.10, and 1.15, respectively. Similarly, 
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FIG, 8,—EFFECT OF CHANNEL DEPTH ON SEEPAGE DISCHARGE 
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tor. For the theoretical case of a channel with zero depth, described in Ap- 
pendix I, the effect of levee base width was determined (Fig. 14). For a given 
B/D ratio, the discharge decreases with increasing levee base width b/D. This 
trend is believed to be generally true for the investigated cross sections. 

Effectof Anisotropy.—Fig. 9 shows discharge as a function of the anisotropy 
ratio of the more permeable layer. In general, the seepage factor increases 
as the anisotropy ratio becomes smaller. This relation holds both for the ho- 
mogeneous cases (Curve I) and for two-layer cases (Curves II, IJ, and IV). 

It was apparent that flow in the semipervious layer, when it is the upper 
layer, is mainly vertical so that its anisotropy ratio is unimportant. Also, 
whenever the semipervious layer lies below the pervious one, the flow trans- 
mitted through it is very small sothat the effect of its anisotropy is again neg- 
ligible. 

The distribution of seepage through the side of the channel depends upon the 
anisotropy ratio. For both the homogeneous and the layered cases with the 
more permeable layer on top, the percentage of seepage through the banks in- 
creases as the anisotropy ratio decreases. When the less permeable layer is 
on top, the fraction of seepage through the banks increases with an increase in 
the anisotropy ratio. Trends are shown on Fig. 10; in some cases the effect 
is negligible. 

The distribution of seepage with landward distance from the levee is shown 
in Fig. 11 for three tests with homogeneous cross sections and for three tests 
with isotropic semipervious layers on top. Although data are limited by model 
dimensions required for transformations, these demonstrate that seepage is 
distributed further inland as the anisotropy ratio decreases. 

Effect of Arrangement and Thickness of Layers.—In comparing the effect of 
the position of the relatively permeable layer to that of the relatively imper- 
meable layer, it was noted that the discharge is greater when the more per- 
meable layer is on top. For the same situation, a greater amount of seepage 
appears near the levee. 

The effect of the relative thicknesses of the two layers is shown in Fig. 12 
for selected tests. These show, for different anisotropies, that the discharge 
decreases asthe thickness of an upper permeable layer decreases| Fig. 12(a)], 
but the converse is true if the upper layer is the relatively impermeable one 
Fig. 12(b)]. 

Effect of Water Table Slope.—The effect was studied in two tests. In both 
cases the water table slope away from the levee was 1:500, but the total drop 
(in feet) represented different percentages of the total head H. A direct com- 
parison between the two tests could not be made; however, individual compar- 
isons were possible between these tests and a test having a similar cross sec- 
tion with a horizontal water table. Compared toa discharge of 1.08 KH for 
the horizontal water table, discharges of 1.3 K H (for a water table drop of 0.3 
H) and 3.88 K H (for a water table drop of 0.6 H) were obtained. Slope of the 
water table also affects seepage distribution landward from the levee. Seepage 
is distributed farther inland with water tables sloping away from the levee than 
for horizontal water tables. The larger the portion of the total available head 
above the lowest point of the slope, the greater the inland extent of seepage. 


CONCLUSIONS 


| Model tests demonstrated that seepage from a leveed river into adjoining 
ow-lying lands is affected to only a minor extent by variations of channel depth 
| ; 

q 
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and width for the conditions studied. The greater the base width of an imper- 
meable levee, the less the seepage. It was found that witha smaller anisotropy 
ratio (greater horizontal permeability) the seepage factor was greater andex- 
tended farther inland. Only rough estimates of anisotropy are necessary to 
define a seepage flow pattern. In two-layered strata, seepage is greater if the 
more permeable stratum is the upper one and increases with the thickness of 
the more permeable stratum. With a water table sloping away from the river, 
the seepage quantity and inland movement increases with slope. 

The electric analogy model was found to provide a rapid and versatile means 
for evaluating seepage under a variety of boundary conditions. The model is 
general in application, being capable of representing an infinite number of ge- 
ologic conditions. Only lack of subsurface data precludes analysis of seepage 
from given cross sections. 

Results defined seepage for a variety of geologic conditions, particularly 
those where two strata govern the flow. By means of a dimensionless scaling 
factor, good estimates can be obtained of seepage quantity and distribution for 
many other cross sections which only approximate test conditions. 


APPENDIX I.—ANALYTIC EVALUATION OF EFFECT OF CHANNEL DEPTH 
ON SEEPAGE FOR d= 0 ANDd=D 


Analytic determination of seepage is practicable for only a very few sim- 
plified cross sections. Solutions are presented here for two isotropic cases, 
one with a channel depth equal to zero (d = 0) and a second with fully penetra- 
ing channel (d = D). These supplement the model results and their interpre- 
tation. 

Case 1 - Channel with Zero Depth (d = 0).—Fig, 13(a) shows the cross sec- 
tion for this case. With regard tothe resulting flow net, Figs. 13(b) and 13(c) 
are equivalent to Fig. 13(a), Fig. 13(c) is then used as a starting point for a 
set of Schwartz-Christofel transformations, The z-plane shown in Fig. 13(d) 
is transformed into the £1 -plane of Fig. 13(e) by means of the transformation 


fy = sin 72 and m and n in Fig. 13(e) are defined by 


2D 
™B -1B 
e2D_ 2D 
m= D) ee . . S Seee 2, eure (10a) 
and 
1(b+B) -7(b+B) 
2D 2D 
-e - . 
re 5 = OE Lathe, Parga ie ke 


The £5 -plane shown in Fig. 13(f) is equivalent to the ¢ -plane shown in Fig, : 
13(e). By this transformation the entire width (2B) of tha channel is consid- 
ered; accordingly, the resulting discharge will correspond to the whole chan- 
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nel. The £3 -plane (Fig. 13(g) is obtained from the £9 -plane by rotating 


The ¢4 -plane[ Fig. 13(h)] is obtained from the £3 -plane by reducing the scale 
by afactor m. Finally, the upper half of the £4 -plane is transformed into the 
inside of the rectangle in the £5 -plane | Fig. 13(i)] by the transformation 


¢ 
colic oie eae ea a 
0 y (1-t?) @-K2 £2) 
£ 
c= [1 t oe eee (12b) 
0 y (1-t2) (1-12 ¢2) 


and 


iE (12¢) 
i K' = SE a c 
1 (1-52) (1-k2 52) 


in which k = = and K and K' are shown in Fig. 13(i). 


The total flow from the channel (2Q) is equal to the flow through the rec- 
tangle of Fig. 13(i) which is 


2Q = (Ke) (2 K) (B) = Ken EE a We ey or (13) 


in which Ke is the hydraulic conductivity of the porous media and H is the head 
difference between the water in the channeland the water table. The quantities 
K and K' are shown in Fig. 13(i). 

A comparision can be made between the theoretical result and model tests 
for the same cross section. As the above considerations apply to isotropic 
media, any given cross section has first to be transformed into an isotropic 
one. The discharge is determined by the ratio & so that the scale can be re- 
duced arbitrarily. Consider the cross section where B=6 in., b = 2 in., and 
D= 20 in. Then 


Be Meee 
Sr OFS: 2D 0.471 
Bt+b _ m(B+b) 
D 0.4 2D 0.628 
0.471 -.-0.471 
k cs = = Se EA = 0.7 re 


00.628 _»-0.628 
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und from Eq. 9a 


sin a =k =0.727; a = 46° 40'; o=5 


' = 0.686; K= 1.880; K' = 1.830 
2q = eee Ke H = 2.06 Ke H; Q = 1.03 Ke H 


[he result of the model test was Q = 1.06 Ke H. For another cross section the 
nodel test gave Q = 0.9 Ke H while by computation Q = 0.935 Ke H. Fig. 14 


tives seepage discharges for various ratios of - and 5 : 


The distribution of seepage with distance from the levee can also be com- 
yuted. The complex potential w = $+ i for the rectangle [ Fig. 13(i)] is 
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wnits-i [ eee eee ee (14a) 
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and 


4et 1 denote the distance from the origin in the horizontal direction [ Fig. 13 
a)]. Inthe transformed z-plane this distance is represented by y along the 
maginary axis. Along this axis 
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The integral must be evaluated in three separate parts: 


1. O<y<B, or0<p<1. For this case 


e dt 
ere (17) 
[, (1-+2) (1-k? t2) 


When y = 0, p = 0, » = 0; and when y = B, p= 1, y= K (the total flow from hall 
a channel). 


FIG, 14,-SEEPAGE FROM CHANNEL OF ZERO DEPTH AS A 
FUNCTION OF b/p AND B/p 
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ee B<y<Btbor1<p<jor fs from K to K+ iK'. Along this line ~ 


remains constant and only the value of @ changes. This follows from the fact 
that the integral has a real part, which is the same as that of part (1), and an 
imaginary part which gives the change of potential along the base of the levee, 


In general, values for ¢ and canbe obtained from w; however, only the values 
of 2 are considered here. 


sD + bay <0, or t< p<, or ¢5 from K+ iK' to iK'. For this case 


sorta | eth a 
y=R if one hae (18a) 
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- fs Bee teee (hy cig eS. 
1/k J (1-t2) (1 -K2 ¢2) 
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where both integrands (and also the integrals) are real in their intervals. E- 
valuation by tables yields 


Die Rie! FC oO eer Seen ose eas ee (19) 
in which 
sin @=1/k pS! pes nla Sint (20a) 
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Thus, for y = 0, » = 0; for y= Bandy = B+b, #=k; and for y = ~, y= 0, 
The following example shows how the preceding analysis can be applied for 

letermining the distribution of seepage. 

_ Example.—At 100 ft from the levee in the same cross section described 

bove, y = 10 in. Width 2 m = 0.978, 


7 


a em /4 = e-77/4 


p Say es <4 1775 
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kc = 0.727; z= 1.372; K = 1.880 


= OOS 
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d 2(0.727) 2 
sin ¢ = 1.316 4a (1.775)*( )* = 0.762 
1 - (1.775) 2 


¢ = 49.6° 
From Eq. 9a 


sin a =k = 0.727; a = 46.7° 
F(0.727; 49.6°) = 0.92 
» = 1.88 - 0.92 = 0.96 


As a percentage the value of ~ equals aa x 100 = 51.2% as compared to ap- 
proximately 53% obtained by the model. 

Case 2 - Fully Penetrating Channel (d = D).— A sketchof this cross section 
appears in Fig. 15. After a series of transformations, similar to those des- 
cribed for Case 1, the following results can be written (25) 


: 
x= | Se, ES eee eee! © 
0 y(1-t2) @-K? t2) 
and 


1/k 
iK' = | Parc eae Re (21b) 


1 ya -s2) (1 -k2 s2) 


The total discharge (for one side) is 


The seepage distribution can be obtained in a similar way tothat described for 
Case 1. 


Fig. 15 shows the relationship between the discharge (expressed as Q 


b Kee 
and the width of the levee | expressed as 2). 


Cases 1 and 2 lead to the following conclusions: 


1. Increasing the widthof the channel increases the discharge very rapidl} 
at first, but has almost no effect for 7 Dk 


2. For a fixed ratio of 5 the discharge increases asthe ratio 5 decreases 


3. For a fully penetrating channel, its width has no effect on the discharge 
The discharge from the channel increases as the lengthof the base of the lever 


. decreases. 


Although these conclusions were obtained from the theoretical cases of : 
channel with no depth and a fully penetrating channel, it is believed that thes 
trends remain true also for the cross sections investigated in the model. Th 
analytic study was based on an aquifer of infinite extent, whereas the mode 
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- ‘a ae, presupposed an aquifer of limited extent. Comparison of re- 
Sults by thetwo methods reveal that differences were negligi justifyi 
ee ak gligible, thus justifying 
In general, the theoretical approach described previously is applicable only 
to a limited number of simplified cross sections ina homogeneous aquifer. 
Even for these cases, the determination of the complete flow net requires many 


2.0 


GHANNEL 


10) 1.0 2.0 3.0 4.0 


FIG. 15,-SEEPAGE FROM A FULLY PENETRATING CHANNEL AS 
A FUNCTION OF b/p 


hours of claculation for each case. From a practical point of view, therefore, 
the electric analogy model was essential to investigate the seepage through 
cross sections approximating field conditions. — 
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FISH PASSAGE THROUGH HYDRAULIC TURBINES2@ 


By Glenn H. Von Gunten,! F, ASCE 


SYNOPSIS 


One of the major problems in connection with water resource development 
in the Pacific northwest of the United States has been safe fish passage at 
dams. Particular concern has been felt in regardto fish safety at dams 100 ft 
or more high. This concept has been perhaps the greatest factor in establish- 
Ment of the basic plan for development of the Lower Columbia and Snake 
Rivers. Also, it may be that the apparent safety accorded to dams under 100 ft 
high has been an adverse factor in migratory fish protection. There is now 
promising evidence that successful passage can be achieved through hydraulic 
heads of several hundred feet, and that the hydraulic characteristics of the 
water passageways are far more significant than heights of dams. . 

This paper reviews basic concepts of fish passage through dams and their 
significance with particular reference given to passage through turbines. Test 
procedures and summarized results of fingerling survival in passing through 
high speed model turbines and a prototype turbine under a 450 ft hydraulic 
head are furnished. The numerous data presented demonstrating a correlation 
of fish survival with various turbine characteristics show promise of leading 
to the achievement of successful passage through high head turbines. 


INTRODUCTION 


Anadromous fish have had and will probably continue to have a great 
influence on water resources development in the Pacific northwest. Hydro- 


Note,—Discussion open until October 1, 1961, To extend the closing date one month, 
1 written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
siety of Civil Engineers, Vol. 87, No. HY 3, May, 1961. 

a Presented at the August 1960 ASCE Hydra, Div. Conf. in Seattle, Wash. 
hi Chf,, Planning and Reports Branch, U.S, Army Engr. Dist., Corps of Engrs,, Walla 
Walla, Wash, 
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electric power developments, primarily on coastal streams and in the middle 
Snake River area, have been extensively curtailed due to the concern of main- 
taining commercial and sport fisheries dependent on those streams as 
spawning grounds and migration passageways to and from the ocean. Curtail- 
ment of such development has been an exceedingly controversial issue, prima- 
rily because of the great economic potential for power development and the 
lack of evaluation of the effects such developments would have on the migratory 
fish population. 

The various species of salmon spawn in fresh-water streams. The young 
fingerlings move from their hatching grounds to the sea where they live until 
maturity. After maturing in the ocean, they return to the the streams to 
spawn and die, thus completing their life cycle. Steelhead trout follow a 
similar life cycle but of a less definable pattern and do not necessarily die 
when spawning. 

Dams without adequate fish passage facilities obstruct the natural passage- 
ways of the migrants. A structure developing a hydraulic head of only a few 
feet may create a barrier to upstream adult migrants, and reservoir im- 
poundments often submerge their spawning grounds. Downstream migrants 
may be exposed to severe physical stresses in passing through hydraulic 
structures, 

UPSTREAM MIGRANTS 


In the past much stress has been placed on the question of the ability of 
upstream migrant adults to negotiate fish ladders. However, during recent 
years the question of upstream passage has been resolved to the extent that 
any one of several methods appear to be feasible. This advance has been 
largely due to the Fisheries Engineering Research Program financed and 
sponsored by the Corps of Engineers, and the activities of other agencies. It 
appears that the primary problem associated with upstream passage is the 
ability to intercept or attract the migrants through the entrances by providing 
proper locations and hydraulic conditions, rather than the ability of the mi- 
grants to negotiate a ladder once found and entered, 


DOWNSTREAM MIGRANTS 


Young fingerling migrants move to the sea with the current. At dams they 
are involuntarily swept through spillways, turbines, or outlet works. In the 
Columbia River, fingerlings vary from 2 in. to5 in. long, except for steelhead 
trout which vary from 5 in to 9 in. Under natural conditions they migrate 
without being exposed to major pressure changes. In being swept through 
hydraulic structures at a dam, they are atthe mercy of the prevailing forces, 
The extreme turbulence and pressure changes of spillways and turbines in- 
flict severe stresses on fish passing through them. Some conditions unques- 
tionably cause injury and death, but for the most part, specific correlation of 
these conditions to distress has not yet been defined. 

There is considerable evidence derived from various research programs 
which indicates or substantiates the following statements: ; 


1. Static pressures associated with hydraulic heads of several hundred 
feet are not fatal, Although under certain conditions distress is evident, no 
lasting injury is apparent. . 
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2. The apparent mortality in passing over or through spillways at dams 
of relatively low head is minor. High free-fall spillways discharging into 
a plunge pool also appear to provide relatively safe passage. 

3, Negative pressures cause distress and may or may not cause permanent 
injury and death. In turbine passage and in laboratories negative pressures 
apparently have caused distress, and in some cases, death. Most observations 
indicate that, under conditions prevailing at dams, lasting ill effects are not 
apparent except in areas of extreme hydraulic turbulence or areas of pressure 
at or near the point of vaporization. The length and nature of exposure to 
both temperature and pressure changes and the degree of vacuum, may all be 
factors of significance. 

4, Exposure to an area of vaporization and cavitation causes death. 

5. Mechanical injury may be experienced by fish passing through turbines 
with small clearances between blades in relation to size of fish, The physical 
injury apparently varies with the relative dimensions and speed, Under proto- 
type conditions, indications are that if blade clearances are greater than the 
length of fish, the mechanical injury is not significant at normal turbine speeds. 


DOWNSTREAM MIGRANT PASSAGE THROUGH TURBINES 


The mortality tests which have been performed at McNary, Big Cliff, 
Baker, Elwha, and Bonneville dams were primarily designed to evaluate the 
extent of, and not to explore the cause of, death. Intensive review of those 
test results and the research data, failed to reveal factual evidence of corre- 
lation between mortality and either high static hydraulic pressures or negative 
pressures normally experienced throughout the major portion of hydraulic 
passageways in a power plant. Because of this, the Walla Walla District, 
Corps of Engineers, initiated investigations toward evaluating mortality in 
relationship to physical conditions. The purpose was to find a means of 
safely passing downstream migrants through turbines. The possibility of 
studying mortality in passing fish through small-scale turbine performance 
models was investigated, and findings were promising. 


MODEL TURBINE TESTS 


In June 1959, exploratory tests were made at the S. Morgan Smith (now 
Allis-Chalmers) Hydraulic Turbine Test Laboratory at York, Pennsylvania, 
o-.determine the feasibility of investigating mortality in model turbines. 

A low-head turbine performance test stand (Fig. 1) was used to conduct 
he tests. The closed hydraulic system was capable of providing turbine 
yperating heads up to 50 ft. Hydraulic control and performance calibration 
ff the model is excellent, as this test stand has been used extensively for 
urbine design research and acceptance tests for many years. The McNary 
2-in diameter Kaplan model runner (Fig. 2), intake, scroll case, wicket 
ates and draft tube were used for the tests. 

A 4-in diameter pipe 42 in. long, with a flap valve at the lower end and 
, screw cap and air line fitting at the upper end, served as a fish lock to 
nove fish into the model. A 4-in pipe, extending downward from the lock, 
ed the fish to a release point approximately midway in thevertical opening 
f the center intake opening. Nylon bobbinet material was used to construct 
net to capture the fish after they passed through the turbine (Fig. 3). 
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FIG, 1 


FIG, 2 
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In nearly all tests, two species of fish were used, generally a hundred 
fish of each species, A similar group of each species, utilized as a control, 
was subjected to the same water, transportation, and handling conditions, but 
did not pass through the turbines. Primary species used were largemouth 
bass (Micropterus salmoides) and fathead minnows (Pimephalus notetus), 
both ranging from 1.5 in to 2.4 in. long. 

Fish and various services were obtained from the Pennsylvania Fish 
Commission. The fish were held at the Commission’s Pleasant Gap hatchery 
and- delivered each testing day to the hydraulic laboratory at York, Pa. At 
the finish of each day’s testing, both test and control fish were returned to 
the hatchery and placed in hatching troughs for holding and recording of 
delayed reactions. A pathologist of the Pennsylvania Fish Commission ex- 
amined test fish at the hydraulic laboratory to record types of injuries sus- 
tained by the test fish under different test conditions. Injuries detected in 
turbine-passed test fish seemed to fall into two categories: strictly mechan- 
ical injuries caused by turbine blades, and internal injuries caused by rapid 
pressure changes. Mechanical injuries ranged from a slight nick some place 
on the body, to complete decapitation or severance of body. Pressure-caused 
injuries were characterized most often by internal hemorrhages, protruding 
eyeballs, and deflated air bladders. 

In general, with the runner setting at tailwater level and normal operational 
conditions free from cavitation, the mortality increased indirect relationship 
to speed. The mortality experienced was largely due to physical injury. 
Artificial braking, or operation at abnormally low speed for a given net head, 
resulted in less total mortality than at the normal high speed. The mortality, 
however, was greater than that experienced at the same speed under normal 
design conditions. Similar tests at the same net heads and speeds, but with 
the turbine runner set at a high level in reference to tailwater, resulted in 
high mortality. The high negative draft tube pressures and associated ex- 
treme cavitation conditions caused extensive pressure injuries fatal to 
nearly all survivors of mechanical injury. A summary of the test results is 
given in Table 1(a) and shown graphically in Fig. 5. In September 1959 the 
Martin power plant Francis-type runner (Fig, 4) was substituted for the 
McNary runner, No significant differences in pressure or death due to 
mechanical injury were experienced at normal speed. However, abnormal low 
speeds resulted in a greater increase in mortality than was experienced with 
the Kaplan runner, Apparently, the Kaplan blades simulated conditions in 
closer proximity to normal speed than could be simulated with the Francis- 
type runner, A summary of the test results using the Francis runner is given 
in Table 1(b) and shown graphically on Fig. 6. 

The tests were highly successful in demonstrating that the research on 
fingerling mortality in passing through turbines is practicable, In addition, 
the tests appeared quite conclusive in respect to the following: 


1. With a small model turbine, 12 in, in diameter, having blade clearances 
less than the length of test fish specimens, wide variations in survival could 
be achieved, depending on speed and tailwater elevation, % 

2. Mechanical injury could be distinguished from pressure injury, and 
was apparently directly proportional to speed, that is, exposure to physical 
contact with blades. 


3. Sufficient survival from mechanical injury was achieved to allow anal- 
ysis of pressure injury, 
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TABLE 1.—TURBINE MODEL TESTS 


Tailwater Speed© In | Percent 

Experiment No. Headwater Elevation Net Head Revolutions | Total 
Elevations | + Centerline In Ft Survival 

Of Runner (Test 


(1) 


A1(K) 12 OCT 59 


+ 0.80 76 
31(K) 13 OCT 59 + 0.75 85 
)2(K) 25 JUN 59 + 0,60 78 
93(K) 14 OCT 59 + 0,95 76 
42(K) 12 OCT 59 + 0.80 96 
33(K) 13 OCT 59 - 0.30 49 
33(K) 14 OCT 59 + 0.30 28 
31(K) 23 JUN 59 + 0.70 21 
32(K) 13 OCT 59 - 0.5 59 
1(K) 14 OCT 59 0. 4? 
31(K) 12 OCT 59 17 
32(K) 13 OCT 59 11 
)1(K) 25 JUN 59 19 
)2(K) 14 OCT 59 38 
t1(K) 15 OCT 59 1 
32(K) 15 OCT 59 10 
03(K) 15 OCT 49 3 4 
24(K} 15 OCT 59 ; ; ; 1300(n) 
\1(F) 20 OCT 59 : 2 313(n) 
32(F) 21 OCT 59 : : 313(n) 95 
\2(F) 20 OCT 59 : : : 939(n) 46 
13(F) 21 OCT 59 ‘ , 939(n) 49 
\3(F) 20 OCT 59 ‘ 939(h) 6 
33(F) 21 OCT 59 ; ‘ ; 939(n) 2 
31(F) 20 OCT 59 : : : 322(r) 21 
J1(F) 21 OCT 59 322(r) 30 
)2(F) 21 OCT 59 : 322(r) 23 
32(F) 20 OCT 59 ; d 315(r) 6 
:1(F) 21 OCT 59 : ; 313(r) 1 
(F) 22 OCT 59 : 3 . 768(n) 69 
'2(F) 22 OCT 59 , 768(n) 23 
\3(F) 22 OCT 59 313(n) 95 
1(F) 22 OCT 59 939(n) 49 
3(F) 22 OCT 59 ‘ F 939(n) 13 
2(F) 22 OCT 59 : 313(r) 15 


_4(McNary Model, 12 in, Diameter. 6 Adjustable Blades, Size Relationship of Model 
> Prototype = 1 to 24) 
_) (Martin Dam Model. 12 in, Diameter, 15 fixed Blades. Size Relationship of Model 
) Prototype = 1 to 15) 

© n - normal speed 
rr - reduced speed 
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4, Extremely high mortality prevailed where adverse hydraulic conditions 
resulted in cavitation, 


The marked increase in mortality associated with increases in negative 
pressures at low heads, together with lack of any evidence that the positive 
head was influential in mortality, encouraged extending the research to a 
high-head prototype installation. 


HIGH-HEAD TURBINE TEST 


The Fishery agency members of the Corps of Engineers’ Fisherie: 
Engineering Research Technical Advisory Committee expressed an interes 
in the results of the model tests and gave support to continued study, whic! 


© NORMAL SPEED - HIGH TAILWATER 
© REDUCED SPEED - HIGH TAILWATER 
4 NORMAL SPEED - LOW TAILWATER 
4 REDUCED SPEED - LOW TAILWATER 


SMALL SYMBOL, INDIVIDUAL TEST 
LARGE SYMBOL, TEST AVERAGE 


% SURVIVAL 


0-1 2 3 4 59 9.6 ia 8 9 10 1 12 13 a 9 
R.P.M. IN HUNDREDS 


FIG. 5.—KAPLAN RUNNER 


led to prototype experiments at a high head power plant. After a compre 
hensive study of power installations in the Pacific northwest, it was decide 
that the city of Tacoma’s (Washington) Cushman No. 2 hydroelectric plan 
was the most suitable project at which to investigate the effects of tailwater 
partial loadings, and the possibility of successful passage of fingerling 
through high-head turbines. Following this investigation a project under th 
Corps of Engineers, North Pacific Division, Fisheries Engineering Researc 
Program was established and conducted by the Walla Walla District. Fig, 
is a view of the Cushman penstock, looking down toward the power plar 
from a point just below the surge tanks. 

The Cushman Project has a hydraulic head of approximately 460 ft. Effi 
ciency of the units is within the range of normal practice and cavitatio 
characteristics are considered normal or possibly subnormal. The turbin 
runner diameter is 76 in, developing 37,500 hp at 300 rpm with a 440-ft heac 
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The runner has fifteen blades. The blade clearance varies from approxi- 
mately 15 in. at the inlet side to 3 in. to 3 1/2 in. at the inner discharge side. 
The draft tube is vertical, discharges intoa rectangular chamber, and is pro- 
vided with a discharge hydraucone, The turbine used in the tests is served 
by a 10.5-ft diameter penstock at anet hydraulic head of approximately 450 ft. 
The turbine discharge is approximately 800 cfs at design head and rated output. 

A series of marks (tattoos) for the fish used in the tests were set up to 
provide a certain mark for each ofthetwelve combinations of tailwater eleva- 
tion and wicket gate opening so that, regardless of when a fish was captured, 
the conditions under which it passed through the turbine would be known. 

Fingerling test fish were inserted in a penstock at El.-402 msl under ap- 
proximately 70 ft of static head through a pressure lock (Figs. 8 and 9). The 
fest specimens were collected by a large net at the draft tube exit, which 
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FIG, 6.—FRANCIS RUNNER 


icreened the entire discharge exclusive of minor leakage around the net 
nchorage frame (Fig. 10). The average recovery of all test specimens placed 
n the penstocks was 90%, The average recovery of all control fish released 
t the center of the net entrance was 82%. Pertinent power plant data are 
iven in Table 2, and test results are summarized in Table 3, Although test 
rocedures preclude refined evaluation, variations in results of duplicate tests 
ppear to be relatively small. Where variations of individual test results ap- 
ear significant, it is attributed to the fact that catches from one test some- 
mes overlapped those of similarly marked fish from an identical test. The 
sh were furnished, marked, handled and examined by personnel of the Wash- 
ton Department of Fisheries and the U.S. Fish and Wildlife Service. 
Under the more favorable tailwater and hydraulic conditions of normal 
irbine operation at best efficiency, and average survival of more than ue 
as experienced, as adjusted by mortalities experienced with control fish, 
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TABLE 2.—PHYSICAL PLANT DATA CUSHMAN NO, 2 
HYDROELECTRIC PLANT UNIT NO, 33 


Plant Location 
North Fork of the Skokomish River near Potlatch, Wash, 


Turbine: 

TY DC a 9 9 a nn Francis 

‘Manufacturer ------------------------------------ Allis-Chalmers 

Normal Head -----------------------------------~ 450 ft 

Rated Output ----=---=—----=-===---—-—--— =~ -=----—— 37,500 hp 

Maximum Head ---------------------------------- 485 ft 

Discharge at Normal Head and Rated Output--------- 800 cfs 

Maximum Tailwater Elev.------------------------- +12,00 msl 

Minimum Tailwater Elev, -------------------------- 0.0 msl 

Draft Tube -------------------------------------- Vertical with Hydraucone 

Speed -------------------- SS Sea Se SS SSE 300 rpm 

Turbine Runner Inlet Diameter -------------------- 83 in, 

Turbine Runner Discharge Diameter --------------- 76 in, 

Maximum Guide Vane Opening --------------------- 5 1/4 in, 

Minimum Clearance Between Runner 

Vanes Varied from ---------------------------- 3 in, - 31/2 in. 

Penstock 

IDR VIN SGI er ea Se te Varies from 10'6" to 8'0"' 

Length - Total ----------=---------===<====---=--~- 1,310 ft 

Length - Lock to Scroll --------------------------- 1,260 ft 


TABLE 3.—TEST FISH SURVIVALS CUSHMAN NO, 2 HYDROELECTRIC PLANT 
UNIT NO, 332 


Gate Opening High Tailwater Medium Tailwater Low Tailwater 
(1) (2) (3) (4) 


2 Average of all tests, 
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FIG, 11.—1960 CUSHMAN NO. 2 TESTS, UNIT NO. 33 
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Under the most adverse conditions of partial turbine loading and low tail- 
water, the survival was reduced to approximately 45%. Figs. 11, 12, and 13 
show relationship of survival to turbine characteristics and tailwater eleva- 
tions. The only apparent significant variables in test conditions were turbine 
Output and elevation of tailwater. The tailwater designations (Table 3) are 
averages of generalized categories associated with high, low, and medium 
tides. Tailwater elevations were further influenced by the discharge of other 
units; consequently, detailed correlation of survival to each of the variables 
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FIG, 12,—1960 CUSHMAN NO, 2 TESTS, UNIT NO, 33 
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cannot be made. Variations in survival of 10% or more were experienced 
with each of the variables individually. 


SUMMARY 


The exploratory tests on survival rates of fish passing through model and 
rototype turbines have definitely demonstrated that turbine performance 
sharacteristics are important, and possibly have a greater influence on sur- 
rival than all other factors. In all test programs, high mortalities have been 
served in passing of fish through turbines operating at low efficiencies and 
through cavitating conditions, Since the extent of cavitation zones in the water 
assages can be highly influenced by design and operational considerations, 
he prospects of improving survival are promising. 
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The exploratory research to date (1961) has indicated definite avenues of 
future study, and also gives much promise that passage through turbines may 
be made relatively safe. The investigations in the future should be aimed at 
three problem areas: 


1. Structural design of turbines to reduce mechanical injury. 

2. Physical setting of the runner in relationship to tailwater conditions to 
reduce negative pressures, 

3. Scheduling of turbine operation to obtain optimum hydraulic conditions 
throughout the water passages, 


Walla Walla District plans to conduct additional turbine model tests with 
fish in October 1960. Additional tests are also being planned for the spring of 
1961 at Cushman No, 2 to explore further both mechanical and pressure in- 
juries and, possibly, to achieve direct correlation between model and proto- 
type results as related to both pressure and mechanical injuries. The primary 
purpose of these tests will be to evaluate the significance of fish size, turbine 
blade clearance, and turbine speed as related to fish survival, 
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FOREWORD 


The purpose of the Task Force on Aerated Flow in Open Channels is to 
investigate the status of present (1961) knowledge of the phenomena of aerated 
flow in open channels and to prepare a report on the subject. The personnel 
of the Task Force and a brief record of its activities are presented under 
the heading “Acknowledgments.” 

The reason for the interest in the phenomena of aerated flow as a mechan- 
ical mixture is involved with the proper design of steep channels or chutes. 
As air is entrained in the flow, the mixture increases in volume of “bulks.” 
The bulking of aerated flow requires higher side walls than non-aerated flow. 
Thus, the phenomenon has important economic aspects in the design of chute 
spillways for earth dams, chutes for wasteways of irrigation canals, and high 
overflow spillways. 

The process of air entrainment in chutes is inextricahly involved in the 
mechanism of the generation of turbulence. The actual presence of large 
juantities of air in the form of bubbles, no doubt, affects the decay of turbu- 
lence. The complete problem is related to the problem of suspended sediment, 
xcept that solid particles tend to settle toward the bottom where turbulence 
is generated and air bubbles rise to the surface. 

_ Although much has been learned of the mechanism of the phenomenon by 
sxperimentation, a rational solution has not been accomplished, in the knowl- 
sige of Task Force members. The objective of this Report is to include 
zmpirical relationships which should be helpful to engineers in preparing an 
sconomical design of a channel to carry aerated flow. Reliable empirical 
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design criteria are not known to have appeared in the technical literature 
published in the United States. 


HISTORICAL RESUME 


A full account of the history of development of the knowledge of aerated 
flow is not within the scope of this Report, nor would it serve any useful 
purpose except for historians. A brief historical resume, however, is of 
value to afford the engineer aperspective ofthe problem, as well as to recog- 
nize outstanding contributions to the advance ofthe general information on the 
subject. A short list of references appears as Appendix I, followed by a more 
extended bibliography in Appendix II. 

Three significant events have come to the attentionof the Task Force mem- 
bers. The first important event concerns pioneer researchin Germany, about 
1925. The next two events involve the initiation of important research in the 
United States and a published symposium on entrainment of air in flowing 
water. The occurrence of the last two events were stimulated in no small 
measure by the interest of the original Committee on Hydraulic Research 
which preceded the organization of the Hydraulics Division of ASCE. 

The Committee on Research of the Hydraulics Division inherited the 
functions of the old committee, upon organization of the Division. In 1957, 
the name of the Committee on Research was changed to the Committee on 
Hydromechanics and a new Committee on Research was formed. It is indeed 
appropriate, therefore, that the sponsoring committee of this task force should 
maintain a continuing interest in the difficult problem of bulking caused by 
aerated flow. 

Instruments to measure the amount of air mixed with the water form an 
important aspect in the advance of the technical knowledge. R. Ehrenberg- 
er (1)1 described a simple, yet effective, device for measuring the concentra- 
tion of air entrained in the flow of chutes constructed in the laboratory at 
Versuchsanstalt fur Wasserbau, Vienna, Austria, His apparatus was simply a 
tube, pointing into the flow, and communicating through an elbow to a vertical 
glass tube of much larger diameter, The Pitot-like tube was rated in unaerated 
flow with a known velocity by measuring the time to fill a known volume of the 
vertical tube. The operation was repeated in aerated flow where the air bubbles 
rose to the surface in the vertical tube. By the increased time to fill the tube 
in aerated flow, Ehrenberger calculated the air concentrationin the air-water 
mixture, 

The Ehrenberger apparatus produced results which have been questioned 
because of the lack of similitude of flow into the Pitot-like tube in aerated and 
in unaerated flow. The apparatus was definitely crude compared with the 
modern electronic method, but ingenious, nonetheless. 


Ehrenberger published empirical relationships defining air concentration 
as follows: 


-0.05 


Py = 0.42 R (sin a) -0.26 for sin a < 0.476 


-0.05 -0.74 


ae =0.30R (sin a) for 0.476 < sin a < 0.707 


I Numerals in parentheses refer to corresponding references in Appendix I, 
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in which Py is the percentage of water in the toal air-water mixture, R denotes 
the hydraulic radius, and @ is the angle of the channel slope with the horizontal, 
Ehrenberger considers his formulas valid for R < 0.30 (0.98 ft). The width of 
his wood chute was 0.25 (0.82 ft). Nodoubt a substantial portion of the air was 
entrained by turbulence generated at the side walls. Therefore, the formula- 
tion of his results would not be applicable to a chute of different width. 

Experimental work was initiated in the United States about 1939. During 
that year J. C. Stevens, F. ASCE, was chairman of the original Committee on 
Hydraulic Research. Upon the recommendation of ASCE, the Engineering 
Foundation made a small grant to Lorenz G, Straub, F. ASCE, Director, St. 
Anthony Falls Hydraulic Laboratory. ASCE had intended that this would be 
“seed” money, to stimulate other sponsorship, which indeed it proved to be. 
Later the funds were supplemented by support from the University of Minne- 
sota, the United States Navy, and other sponsors. 

The early adjustable metal flume at St. Anthony Falls was 1 ft wide. In an 
effort to achieve true two-dimensional flow, with no sidewall effect on the air 
entrainment, the flume was rebuilt to a width of 18 in. in 1950, 

The most important advance in the science of aerated flow research at 
st. Anthony Falls Hydraulic Laboratory was the development of better methods 
of measuring the concentration of entrained air. The early method of measur- 
ing the entrained air volumetrically was to give way to an electronic method 
ia 1950. This apparatus was described in a paper by O. P. Lamb and J. M. 
Killen (2). Essentially it involved a pair of electrodes set in the line of the 
flow direction. The electrical resistance of the laboratory water without 
entrained air was measured, in effect. The concentration of air in aerated 
flow could then be calculated by measuring the increase of resistance across 
the electrodes. 

Two series of experiments were conducted with the wide flume and the 
electrical resistance apparatus. Results of a series with the smooth steel 
wall were in an unpublished report by Straub and A. G. Anderson (3), for 
the United States Navy in 1955, The results of the test series on walls 
roughened with granular particles were reported by Straub and Anderson (4) 
in 1958. Significantly, Straub was awarded a Research Prize by ASCE in 
1958 for his work on aerated flow. 

A Symposium on the Entrainment of Air in Flowing Water was published 
in the Transactions of ASCE in 1943. The paper by L. Standish Hall (5) re- 
sorted on field tests of three chutes operated by the Pacific Gas and Electric 
company. He also included field test data on the Kittitas chute by the Bureau 
of Reclamation, United States Department of the Interior (USBR). For reasons 
yf large width and a mimimum amount of vertical and horizontal curvature, 
he Kittitas results proved the most desirable for later analysis. 

The tests on the Kittitas chute were made by C. W. Thomas and reported 
o the USBR in 1938 (6). Thomas measured the water-surface profiles as 
yell as the mean velocity by the salt-velocity method, With a known water 
lischarge, the quantity of entrained air could be estimated. 


PRESENT STATUS OF KNOWLEDGE 


Narrow Chutes.—For purposes of this Report, a narrow chute is defined 
iS one which has a width less than five times the depth. Air entrainment can 
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ordinarily be considered to be a relatively simple two-dimensional problem 
for wide chutes, Three boundaries, the bottom andtwo sides, contribute to the 
generation of turbulence which causes air entrainment in narrow chutes. It 
is probale that the effect from one side wall reaches nearly to the opposite 
wall, thus compounding the air-entrainment effect in narrow chutes. Fig. 1 
demonstrates development of the side-wall effect of the Uncompaghre Chute 
very clearly. Much of the data reported by Hall (5) fall in the category of 
narrow chutes, The discussions of this paper embraced attempts to generalize 
the results for the calculation of air concentration and the resultant bulking of 
the flow. 


FIG. 1.—AIR ENTRAINMENT FROM SIDE- 
GENERATED TURBULENCE 


FIG, 2,~AIR ENTRAINMENT CAUSED BY SIDE-WALL 
EFFECT IN NARROW CHUTES 
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In later years, other attempts at generalization of the air-entrainment 
phenomena were based on the Hall paper. Although these empirical conclusions 
were roughly applicable to narrow chutes, they lead to over-design of spillway 
chute side walls when applied to wide chutes. Fig. 2 shows the side-wall effect 
of the USBR Kittitas Chute. 

Wide Chutes.—Maintaining the definition for a narrow chute expressed 
previously, a wide chute is considered to be one which has a width greater 
than five times the depth. An example of a wide chute may be seen in Fig, 3, 
which shows the Ft. Peck spillway. Close scrutiny of the water surface during 
operation led to the conclusion that there was a negligible amount of bulking 
from air concentration. Such an observation was impressive when contrasted 
with a high air concentration indicated by previous empirical criteria based on 
narrow chute operation. 


FIG, 3.—LOW AIR ENTRAINMENT FROM BOTTOM-GENERATED TURBULENCE 
IN A WIDE SPILLWAY CHUTE 


Hydraulic design engineers have been concerned that the use of design 
sriteria based on unselected field data for narrow chutes had influenced the 
endency toward costly over-design of side walls. In recent years, engineers 
of the United States Army Waterways Experiment Station, Vicksburg, Miss., 
(WES), analyzed the smooth-channel data obtained at St, Anthony Falls (3) in 
sombination with selected Kittitas field observations, The Kittitas data were 
selected on the basis of the observations with depths less than one-fifth of the 


vidth. 
' . 
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Straub and Anderson had found the air concentration, c, in two-dimensional 
smooth channel to be a function of S/ q2/ 3, in which S is the sine of the slope 
and q.is the unit discharge. 

The WES issued a hydraulic design criteria chart in June, 1957, based on 
the foregoing analysis of the Straub and Anderson smooth-channel results. 


The curve of best fit was 


c . 


il} 
oO 
1) 
co 
_ 
(2) 
oq 
ren 
oO 
Pea 
alin 
a 
bd 
Sa 
oo 
+ 
oO 
a 
-~] 
~-~ 
my 
— 


This chart presented a much more economical basis for the design of wide 
chute spillway walls than did previous criteria from unselected narrow chute 
results. 

As mentioned previously, Straub and Anderson (4), presented the results 
of tests on rough channels in December, 1958. They concluded that for an 
artificially roughened channel, air concentration was a function of S/ql 5, 
This function of slope and unit discharge was based on an empirical rela- 
tionship between ¢ and V*/ d, in which V* is the shear velocity and dt isa so- 
called transition depth as defined therein. The new Straub and Anderson data 
on rough channels were also analyzed in comparison with selected Kittitas 
field results. The large unit discharge of the prototype chute at Kittitas pro- 
duced small values of S/ qi 5, and consequently small concentrations of air. 
The curve of best fit for the combined laboratory and field data was deter- 
mined by use of a digital computer andis shown in Fig, 4. The standard error 
for air concentration was 0¢ = 0,061. The Task Force believes that this rela- 
tionship offers a realistic basis for the design of wide chute spillways within 
the limitations of the present status of knowledge. 

Overflow Spillways.—The relationship represented by the curve in Fig. 4 is 
believed applicable to wide ungated overflow spillways, although no prototype 
observations are yet available for verification. A laboratory study by W. J. 
Bauer (7) made at the State University of Iowa, Ames, Iowa, covered an 
investigation of the development of the turbulent boundary layer on an over- 
flow spillway. His observations were not carried beyond the point on the pro- 
file where the turbulent boundary layer reached the surface. It is at this part 
where the phenomenon of air entrainment begins. ’ 

The Task Force is not aware of any research which deals with the effect 
of spillway crest piers on air entrainment. Examination of the operation of 
overflow spillways in the field indicate that the turbulence generated at the 
pier noses can cause substantial air entrainment. Indeed, it is believed that 
for modern spillways with a large design head compared with the width of the 
gate bay, the pier disturbance has a greater effect on air entrainment and its 
resuitant bulking, than does the bottom turbulence alone. The effect of piers on 


air entrainment at the Chief Joseph Dam on the Columbia River may be seen 
in Fig. 5. 


NEED FOR FUTURE RESEARCH 


This Report has presented an empirical criterion as a guide to the judg- 
ment of the engineers in the design of wide steep channels. Nevertheless, 
much remains to be learned of the basic mechanics of the generation of turbu- 


lence, its energy as the water surface is approached, and of the actual 
phenomenon of air entrainment. 
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STRAUB AND ANDERSON RANGE, 74 LABORATORY 
OBSERVATIONS ON A ROUGH CHANNEL 


KITTITAS RANGE 
18 SELECTED 
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FIG, 4,—AIR ENTRAINMENT BASED ON TESTS OF A ROUGH 
LABORATORY CHANNEL AND FIELD TESTS 
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Theory.—The problem is related to the more fundamental one of the genera- 
tion of turbulence at the boundary, its diffusion upward through the liquid, and 
the action of the vertical components of turbulence at the water surface. More 
needs to be known of the reflection or dissipation of the vertical components at 
the surface with various degrees of turbulence and the actual piercing of the 
surface by turbulence cells. It is undoubtedly the latter phenomenon which gives 
rise to air entrainment. 

The present (1961) knowledge of turbulence is based primarily on the great 
mass of observations with the flow of air. The work in aeronautics is not 


FIG, 5,—EFFECT OF SPILLWAY PIERS ON AIR ENTRAINMENT 


concerned with the phenomenon of turbulence in the vicinity of a liquid with 
a free water surface. Unfortunately, the hot wire anemometer, which is the 
principal measuring device in wind tunnel turbulence studies, is not readily 
adaptable to use with air-water mixtures. Nevertheless, theoretical studies 
of turbulence in open channels should start with a study of the theory and test 
results of aerodynamics, 

Laboratory Equipment.—Although the combined effect of both side walls 
and bottom are important in certain structures, the more simple two-~ 
dimensional problem of a wide channel should be approached first. A dis- 
advantage of studying wide channels with steep slopes is that large discharges 
are needed to obtain substantial depth. The fundamental problem of turbulence 
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in open channels can be studied with fairly light slopes, Practically uniform 
flow can be obtained with a reasonable channel length, if the bottom is artifi- 
cially roughened. A fairly wide channel with a rough bottom and smooth side 
walls would then have very little side-wall effect, 

More laboratory work is needed to perfect instruments which will measure, 
at least, the longitudinal and vertical components of turbulence velocities in 
water. The Hubbard-Ling hot-film anemometer (8)(9) is one approach and 
several laboratories have experimented with pressure transducers mounted 
on stream-lined supports. 

The three-dimensional problem could well be investigated by first studying 
the side-wall effect, alone. A wide laboratory flume is again indicated to 
minimize reflection from the opposite wall. The experimenter would probably 
choose to roughen only one wall and maintain the bottom and opposite wall 
smooth, 

The behavior of the turbulence from discrete wall protuberances, whether 
knobs or vertical strips, as the disturbance reaches the surface needs to be 
studied, Little is known of the effect of joint offsets in channel walls. The effect 
of piers on overflow spillways belongs to this class of problems. The Task 
Force believes that the fundamental research could best be performed in the 
laboratory, where the geometrical variables may be changed with little 
expense, 

Field Tests.—More observations on existing prototype structures should 
be made and publicized, if only to alert the profession to the important design 
problems involved. The problem of holding measuring instruments in high 
velocity flows of the prototype is severe in itself. Provisions for measure- 
ment of such a rudimentary variable as velocity are badly needed in the pro- 
totype. Velocity observations would yield information for friction-loss studies 
as well as for basic information onhigh velocity turbulent flow. The provision 
for high-velocity measurement must normally be done during construction, 
by embedded metal which canholda support at the bottom or a span for attach- 
ing the support from above, Embedded metalhas been installed for the attach- 
ment of instrument piers to the bottom of the Ft, Randall chute spillway. Ob- 
Servations at this installation must await the availability of test flows, which 
is a common problem with all field tests. 

Air entrainment caused by the disturbances of crest piers on overflow 
spillways may be a critical condition for the design of side walls of this type 
of spillway. General field observations of the pier effect may well serve as 
a useful guide in the design of future spillways, until more basic work can 
be accomplished in the laboratory. 
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FORMS OF BED ROUGHNESS IN ALLUVIAL CHANNELS# 


By D. B. Simons, ! M. ASCE, and E, V. Richardson,2 A.M. ASCE 


SYNOPSIS 


Laboratory experiments in a large recirculating flume, in addition to field 
studies, have established that resistance to flow and sediment transport in 
alluvial channels are related to the form of bed roughness. The forms of bed 
roughness that have been observed in the laboratory flume and in natural 
alluvial channels can be divided into two regimes of flow, on the basis of their 
Shape, resistance to flow, and sediment transport. These divisions are the 
lower flow regime and the upper flow regime. Between the two regimes there 
is a transition zone. The bed forms are ripples and dunes in the lower flow 
regime and plane bed, standing waves, and antidunes in the upper flow regime. 
The bed form in the transition zone is variable, ranging from that which is 
typical of the lower flow regime to that for the upper regime. 

Resistance. to flow in the lower flow regime is relatively large and sedi- 
ment transport is relatively small. Conversely, in the upper flow regime the 
resistance to flow is relatively small and sediment transport is large. In the 
‘ransition zone both resistance to flow and sediment transport vary appreci- 
ably. 

The variables that determine the form of bed roughness are the slope of 
he energy grade line, the depth of flow, fall velocity or effective fall diameter 
of the bed material, and the shape of the channel. The fall velocity is a func- 
ion of the size of bed material, the shape factor of the sediment, the mass 
lensity of the sediment, the mass density and the viscosity of the water- 
sediment mixture. 


i a ea ee 
Note.—Discussion open until October 1, 1961, To extend the closing date one month, 
1 written request must be filed with the Executive Secretary, ASCE. This paper is part 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
siety of Civil Engineers, Vol. 87, No. HY 3, May, 1961. 
@ Publication authorized by the Dir., U. S. Geol. Survey, Dept. of Interior. 
_ 1 uHydr, Engr., U. S, Geol. Survey, Colorado State Univ., Fort Collins, Colo. 
_ 2 yYydr, Engr., U. S, Geol, Survey, Colorado State Univ., Fort Collins, Colo, 
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Laboratory and field studies indicate that the change in resistance to flow 
that occurs when the form of the bed roughness changes from dunes to plane 
bed or standing waves, accounts for the discontinuity observed in stage dis- 
charge relationships on those streams that flow in both regimes during a 
runoff event. 


INTRODUCTION 


Flow in alluvial channels is complicated by the fact that the form of the 
roughness elements and, hence, the resistance to flow and the sediment trans- 
port are functions of the characteristics of the fluid, flow, bed material, and 
channel, Because of the limited scope of most studies and the interdependency 
and the number of variables involved there are many conflicting statements in 
the literature. This was brought out by Norman H. Brooks,3 A.M. ASCE and 
those who discussed his paper. However, a clearer insight into the alluvial 
channel problem results when the various forms of bed roughness that can 
occur in an alluvial channel, the factors that affect the form of bed roughness, 
and the variation of resistance to flow and sediment transport with forms of 
bed roughness are known, 

Herein is described the bed configurations that have been observed ina 
large recirculating laboratory flume with a sand bed and in alluvial channels 
in the field; the magnitude of resistance to flow and sediment transport as- 
sociated with these bed forms; the variables that influence the bed configura- 
tion and, where known, how these variables influence the bed form, A possible 
relationship that may prove useful in predicting the configuration of the bed 
is introduced and one field situation for which the change in bed form is of 
prime importance is described, 

This analysis is an outgrowth of a continuing study of alluvial channels 
by the United States Geological Survey, Dept. of Interior (USGS) at Colorado 
State University, Fort Collins, Colo. The investigation is conducted primarily 
in a recirculating flume 150 ft long, 2 ft deep, and 8 ft wide. Slope can be 
varied from 0 ft per ft to 0.01 ft per ft and discharge from 2 cfs to 22 cfs. To 
supplement the large flume a recirculating flume 2 ft wide, 3 ft deep, and 60 
ft long is also used. 

The depth of sand bed in the flumes was approximately 0.6 ft. With shal- 
lower depths of sandthe dunes could not fully develop and portions of the flume 
floor were exposed, The 0.45 mm data quoted in this report were previously 


3 “Mechanics of Streams with Movable Beds of Fine Sand,” by Norman H, Brooks, 
Transactions, ASCE, Vol, 123, 1958, pp. 526-549, 
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|published4 and the 0.28 mm data are included herein (Table 1). A more detailed 
description of the facilities and methods of collecting the data are available.° 


| REGIMES OF FLOW 
: 


It is commonly recognized that four regimes of flow occur in open channels 
iwith rigid boundaries.6 Two of these regimes frequently occur in natural 
jalluvial channels. These are the tranquil-turbulent flow and rapid-turbulent 
jflow regimes. Tranquil and rapid flow are adequately defined by either the 
|Specific energy diagram or the Froude number, F. When F < 1, the flow is 
|tranquil and when F > 1, the flow is rapid, The Froude number is defined as 


Bhi ahd SRE I A ane RC SA (1) 


ve D 
‘in which V is the velocity of flow, in feet per second; g denotes the gravita- 
tional acceleration, in feet per second squared; and D refers to the depth in 
feet. 

Normally the mean velocity V and average depth D are used to compute 
the Froude number, However, it is important to recognize that with the ex- 
treme variability of flow conditions that occur in the cross section of an 
alluvial channel it may be advantageous to consider local values of velocity 
and depth and the corresponding Froude number to help explain observed 
phenomenon, It is not uncommon to be confronted with a situation where the 
Froude number, based on average V and D, is less than unity and yet in the 
same cross section, the local Froude number may exceed unity. 

The absolute magnitude of the Froude number in either the tranquil or the 
rapid flow regime depends onthe scaleof the system and is only quantitatively 
significant for the system under consideration. For instance, in the large re- 
circulating flume which is 8 ft wide, the dunes only occurred when 0.3< F<0.6, 
whereas in a large deep river, dunes can occur whenF << 0,3. Also in the 
large flume, the beginning of motion occurs at aF ~ 0.15, whereas, in a 
very small flume using the same bed material and width-depth ratio, the be- 
ginning of motion may occur at F > 1.0. 

If F < 1 the water accelerates over the artificial or natural humps on the 
stream bed and decelerates over the depressions or troughs. This is illustrated 
in Figs. 1(a) and 1(b). Fig. 1(a) also illustrates a separation zone and the 
existence of strong circulation in the trough area of a dune bed and the fact 
that boils appear on the water surface just downstream of the crests of the 
sand waves. The separation and the turbulence generated dissipates consid- 
erable energy which increases the resistance to flow. 

When F > 1, the water decelerates over the humps and accelerates in the 
troughs as illustrated in Figs. 1(c) and 1(d). Fig. 1(c) illustrates symmetrical 
sand and water waves of sinusoidal form that have been described mathemati- 
cally with limited accuracy by the writers, In flumes these sand and water 


4 “Resistance to Flow in Alluvial Channels,” by D. B. Simons and E, V. Richardson, 


Beoccedings, ASCE, Vol. 86, No. HY 5, May, 1960, 
_ 9 “Studies of Flow in Alluvial Channels; Flume Studies Using Medium Sand,” by D, 
B. Simons, E. V. Richardson, and M, L. Albertson, Water Supply Paper 1498A, U. S. 
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6 “On the Four Regimes of Open-Channel Flow,” by J. M. Robertson and Hunter 
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waves are commonly observed in the rapid flow regime. There is little sepa- 
ration and circulation in the rapid flow regime when plane bed and standing 
waves exist. In these cases, the dissipation of energy reflected by resistance 
to flow is primarily the result of shear on the bed and the formation of waves. 


TABLE 1.—BASIC VARIABLES FOR RUNS WITH 0.28 mm SAND IN THE 8 FT FLUME 


Viscosity | Tem- 
in ft per] x 105 ft2 


Run Slope x | Dis- 
10“per | charge 
ft cfs 
(1) (2) (3) 
6 0.005 3.79 
if 0.007 6.61 
8A 0.011 7,76 
8B 0.023 7.76 
10 0.041 4,16 
5 0.045 10.73 
13B 0,062 13,46 
4 0.069 10,73 
thik 0.073 4.92 
33 0.090 15,74 
12 0.100 12.70 
12 0,108 TA9 
14 0,116 8.61 
20 0.120 18,14 
ie 0.131 15,19 
21 0,131 20,39 
19 0,134 9,90 
162 0.134 17,23 
23 0.134 22,02 
174 0.136 10,01 
34 0.136 15,28 
18 0,141 11.96 0.61 
30 0,142 15,68 0.64 
34 0,150 5.50 0.44 
22 0.153 14,92 0.60 
15@ 0,158 L2KS8M a) MORTS 
24 0.172 21,98 0,82 
25 0,199 21.85 0.72 
28 0,229 15,72 0.55 
29 0.278 15.70 0,52 
26 0,328 15.51 0.50 
32 0.470 21.76 0.58 
27 0.533 15.47 0.43 
31 0.593 21,34 0,56 
35 0.815 21,33 0,54 
37 0,820 8.34 | 0.30 
38 0.930 15,26 0.40 
36 1,007 21.38 0.57 


2 Slope poorly defined. 
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With antidunes the resistance to flow is related to the shear on the bed, the 
‘formation of waves, and the energy dissipated in the breaking waves. 


: 
FORMS OF BED ROUGHNESS 
_ The forms of bed roughness observed inthe flumes and in alluvial streams 
are illustrated in Fig. 2. These forms of bed roughness were described and 
divided into the tranquil and rapid flow regimes," Additional laboratory 
istudies by the writers and a field investigation, ‘ indicate that it is more logi- 
cal to classify forms of bed roughness into a lower and upper regime con- 
mected by a transition zone than to work with the tranquil and rapid flow re- 
gimes that are defined by the Froude number. 

The lower and upper flow regime classification isbased upon magnitude of 
resistance to flow, bed, and water surface configuration, and mode and magni- 
tude of sediment transport. In the lower regime the bed roughness is ripples 
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“Rapid Flow, Alluvial Channel Rapid Flow 
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FIG, 1.—-RELATION BETWEEN WATER SURFACE AND BED CONFIGURATION IN 
THE TRANQUIL FLOW AND RAPID FLOW REGIMES 


or dunes or both, the resistance to flow is large and the bed material load is 
relatively small. The majority of the bed material load moves close to the 
ed over the backs of the ripples or dunes and at the crest of these roughness 
slements some of the bed load avalanches down the faces of the ripples or 
junes where it is temporarily stored causing the ripples and dunes to advance 
lownstream. The other part of the bed material load is carried onward, some 
© the surface in the boil areas. 

When the bed is plane or when standing waves and antidunes develop the 
‘esistance to flow is small, the bed material transport is relatively large 
ind is nearly continuous in motion. With ripples and dunes bed load is closely 
‘elated to the velocity of the roughness elements. With plane bed, standing 


7 Personal communication with D, R. Dawdy. 
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waves, or antidunes the bed load is closely related to the velocity of the in- 
dividual grains moving in contact with the bed. In the transition region which 
connects lower and upper flow regimes there is a discontinuity in the re- 
sistance to flow and sediment transport relations. 

There are reports in the literature of aplane bed with bed material move- 
ment which develops before ripples in the lower regime of flow. Based upon 
experimentation, this seems to be a small flume phenomenon. In the large 
flume whenever general movement of bed material occurred (d< 0.5 mm) 
ripples or dunes ultimately formed. 

Ripples are usually spaced on the order of 0.5 ft to 1.5 ft and have ampli- 
tudes of 0.01 ft to 0.1 ft. The spacing of the dunes is usually greater than 2 ft, 
and their amplitudes range from 0.10 ftto many feet depending on the depth of 


(a) Typical ripple pattern, F «| (e} Plone bed, F-<1.and.d¢0.4mm 
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FIG, 2.—FORMS OF BED ROUGHNESS IN ALLUVIAL CHANNELS 


flow and the characteristics of the bed material. For example, in small al- 
luvial channels the dunes may be 0,5 ft to 1 ft high with a spacing of 5 ft to 
20 ft, whereas, in the Mississippi River sand waves of dune form have been 
recorded by W. C. Carey, F. ASCE andM. D, Keller,8 A.M. ASCE which have 
lengths of several hundred feet and amplitudes as large as 40 ft. Dunes as 
defined here are frequently referred to as bars. 

In the transition zone for the same boundary shear the bed configuration 
may range from plane bed to fully developed dunes depending primarily on 
antecedent conditions. Generally, the bed consisted of long dunes of small 
amplitude that decreased in amplitude as the boundary shear was increased, 
However, when the bed form consisted of dunes, slope or depth could be in- 
—-)YgXKs 


8 “Systematic Changes in the Beds of Alluvial Ri ca 
vers,” by W.C. Car 
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' 
creased to relatively large values before the bed form changed to the plane 
bed or standing waves of the upper flow regime. Conversely if the bed was 
plane, the slope or depth could be decreased significantly before dunes de- 
veloped. 

__ A plane bed developed in the flume experiments following the transition 
zone with the 0.28 mm sand, With the 0.45 mm sand a plane bed did not de- 
velop following the transition. Instead, standing waves developed. The standing 
sand and water waves that occurred with the 0.45 mm sand and that did not 
occur with the 0.28 mm sand were of sinusoidal form, in phase, and did not 
move either upstream or downstream with time to any great extent, The 
amplitude and spacing of these waves depended upon the size of the flow 
system and the characteristics of the fluid and the bed material. The ampli- 
tude of the water waves increased with depth and were 1.5 to 2 times the 
amplitude of the sand waves. 

The antidunes are very similar to standing waves except some of the sand 
and water waves continued to grow in amplitude until they broke, Antidunes 
have been observed in natural streams in which the median diameter of the 
bed materials have ranged in size from fine sand to coarse gravel. In the 
flume studies the breaking of the antidune waves usually occured when the 
amplitude of the water waves were about twice the amplitude of the corre- 
sponding sand waves. At this time, the water surface in the trough was at 
approximately the same elevation as the crests of adjacent sand waves, These 
sand and water waves moved upstream prior to breaking, Either one or two 
waves broke at onetime oratrainof several waves broke more or less simul- 
taneously. After the antidune waves broke, anewtrain of waves developed and 
the antidune cycle repeated itself or the waves subsided without developing 
0 their break point and then reform to break or abate again. As indicated by 
W. B, Langbein,? F, ASCE the Froude number, where antidunes first form, 
nereases with increasing size of bed material. 

In the large flume at the steeper slopes, the antidune activity was in the 
orm of chutes and pools. Flow was rapid and accelerated in the chutes and 
low was tranquil in the pools, as in Fig, 2(h). 

When antidunes break they resemble the hydraulic jump and can be analyzed 
vith some success as such, particularly with two-dimensional flow. Large 
juantities of sediment are carried into suspension in the breaking antidune 
vave, 

In the breaking wave and immediately downstream, the velocity is rela- 
ively small. The storage of water in the short reaches where the waves are 
yreaking causes the discharge to vary with time and this action tends to set 
ip slug flow. In Mendano Creek, which is located in the San Luis Valley of 
solorado, antidunes set up slugs of water thathad a depth about twice average 
lepth and traveled down the channel spaced at intervals of about 350 ft at 
early 10 fps. At the time of this field observation, the average discharge 
yas approximately 120 cfs, the median diameter of the bed material was 0.3 
am, and the channel slope was 1.67%. 


VARIATION OF RESISTANCE TO FLOW WITH BED ROUGHNESS 


4 


_ The variation of resistance to flow, as measured using the Darcy-Weisbach 
, with the form of bed roughness and the influence of size of bed material 


9 : . . 
_ 9 “Hydraulic Criteria for Sand Waves,” by W. B. Langbein, Transactions, Amer. 
eophysical Union, 1942, pp. 615-618, 
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and slope on the form of the bed roughness, resistance to flow, and sediment 
transport are illustrated in Fig. 3 and Table 2. 

The. resistance to flow is relatively small for plane bed prior to beginning 
of motion, However, resistance to flowis larger for this case than it is for the 
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plane bed or standing wave case where there is motion of bed material. This 
is understandable because, in the first case, resistance results from grain 
roughness that is fixed; whereas, in the second case the grain roughness = 
not fixed but moving and a form of rolling friction exists, 
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As general motion of the bed material was initiated by increasing the 
boundary shear, ripples formed that caused a large increase in resistance 
to flow. With sufficient increase in slope, and depth constant, dunes formed, 
and with further increase in slope transition conditions developed. The ripples 
formed from the 0.28 mm bed material had a slightly larger resistance to 
flow than the ripples formed of 0.45 mm bed material, Conversely, the dunes 
formed from 0.45 mm bed material had a greater resistance to flow than the 
dunes formed of 0.28 mm bed material (Table 2). 

In the transition zone, because the bed form was unstable, the resistance 
to flow was extremely variable. In general, the magnitude of the resistance 
to flow was between that for fully developed dunes and a plane bed. The vari- 
ation in resistance to flow in the transition zone was larger for the 0.45 mm 
bed material than for the 0.28 mm bed material. Also note that the transition 
zone occurred at a much flatter slope when using the 0.28 mm sand, 

With still further increase in boundary shear the upper flow regime was 
reached. In this regime of flow the bed roughness was plane bed, standing 
waves, antidunes or some combination of these. Resistance to flow was a 
minimum for a given bed material throughout the plane bed or standing wave 
range or both. As antidunes develop, the resistance to flow increased slightly 
and continued to increase with further increase in boundary shear. The rate 
of increase in resistance to flow was much more rapid for the finer sand and 
was reflected directly by increased antidune activity. 


VARIATION OF BED MATERIAL LOAD AS A FUNCTION OF BED FORM 


The large flume studies indicated that the bed material transport varied 
with form of bed roughness and size of sand, as shown in Fig. 3 and Table 1. 
The transport for the ripples reached a maximum of 150 ppm for the 0.28 
mm bed material and 100 ppm for the 0.45 mm bed material. With dunes as 
the form of bed roughness, the maximum concentration of bed material trans- 
port was about 800 ppm and 1200 ppm, respectively, for the 0.28 mm sand 
and the 0.45 mm sand, 

For the 0.28 mm sand, the bed material load in the upper flow regime 
ranged from 1500 ppm for the plane bed to 42,000 ppm with antidunes. The 
range in concentration of bed material discharge for the 0.45 mm sand was 
from 4,000 ppm with standing waves to 15,000 ppm for the antidunes. 

The bed material load in the transition zone was unpredictable but varied, 
for both sands, from slightly greater than the maximum concentration for 
the dunes to a concentration that was characteristic of that for the plane bed 
or the standing waves. 


FACTORS INFLUENCING FORM OF BED ROUGHNESS 


The variables that fix the form of bed roughness in alluvial channels are 
Form of bed roughness = @ [D, S, do, Pp, g, w, Ss, fs| Re es Sok (2) 

in which 
word’ = $[d, Pg, P, g, Syu]..... Pg ~(z 


and D refers to the depth, S denotes the slope of the energy grade line, 
defines the median diameter of the bed material, o is a measure of the size 
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te, P refers to the mass density of the water-sediment mixture 
Ormally taken as that for water unless considerable fine sediment is pres- 
mt), p g is the mass density of the bed material, g denotes the acceleration 
f gravity, w is the fall velocity of the bed material, f, describes the seepage 
ree, Ss is a shape factor of the cross section, Sp is a shape factor of the 
ticle, 4 defines the viscosity of the water-sediment mixture, and d’ refers 
) the effective median fall diameter of the bed material. In addition to tem- 
rature, 4 also depends on the concentration of fine sediment, 
_ The full significance of all of these variables has not been revealed (as of 
961) because the investigations have not been broad enough to include all 
onditions. However, the effect of some of the variables are at least quali- 
atively known and in the following paragraphs they are summarized, Re- 
erring to Eq. 1, velocity or bed material load can be substituted for bed 
oughness., 

Slope S.—The effect of slope was illustrated in Fig. 3. With depth constant 

change in slope increases the boundary shear. With sufficient increase in 
hear the form of the bed roughness changes, To develop a given form of bed 
oughness in alluvial channels that has the samebed material, the slope must 
e large if the discharge is small if the discharge is large. 

Depth D.—The effect of depth is not well defined because the range in 
epth which has been studied is small, However, in the experiments conducted 
1 the large flume the variation in depth has been adequate to indicate that a 
hange in depth can change the bed form. In natural alluvial channels in which 
lope is essentially constant bed form changes have been reported that have 
esulted from changes in depth,10,11,12 However, there appears to be combi- 
ations of slopes and sizes ofbed material where the bed form will not change 
ith a large change in depth. 

In addition to the changes in bedform that can occur with changes in depth, 
ere are the uncertain effects of changes of depth on bed material load and 
=sistance to flow when the bed form does not change. With ripples there is a 
lative roughness effect. The resistance to flow decreases with increase in 
spth. Whereas with dunes, the flume data indicate that resistance to flow 
creases with increase in depth. However, it has not been verified that re- 
stance to flow with a dune bed form increases continuously over a large 
unge of depths. 

Bed material load is also affected by changes in depth of flow. Based upon 
recent study of available field data, B. R. Colby!3 has shown that because 
locity increases with depth and concentration decreases with depth, a 
culiar but logical relationship exists between sediment load, depth of flow, 
id velocity. Within a certain range of velocities, if velocity is held constant, 
e bed material load increases with increasing depth, within another range 

velocities, bed material load decreases with increasing depth. 


~ 


10 “Discontinuous Rating Curves for Pigeon Roost and Cuffawa Creeks in Northern 
ississippi,” by B. R. Colby, ARS 41-36, U. S, Dept. of Agric. 

1 Discussion by J. K. Culbertson and C.F, Nordin of “Discharge Formula for Straight 
luvial Channels,” by Hsin-Kuan Liu and Shoi-Yean Hwang, Proceedings, ASCE, Vol. 
, No. HY 6, June, 1960. 
12 «studies of Flow in Alluvial Channels; Depth-Discharge Relation of Alluvial 
reams,” by D.R. Dawdy, Water Supply Paper 14896, U.S. Geol, Survey, 1961, in press. 
13 personal Communication, “ 
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The effect of the combination of depth and slope in the shear term, Y DS 
has been studied by I. A. Shields,14 R, A. Bagnold15 and Hans A, Einstein 
F. ASCE, and Ning Chien, M. ASCE, 16 to mention a few. All have indicatec 
that shear is an important parameter, whereas Brooks? stated “that for the 
laboratory flume it was found that neither the velocity nor the sedimen 
discharge concentration could be expressed as a single-value function of th 
bed shear stress, or any combination of depth and slope... .” It is apparen 
in Fig. 3 that what Brooks states is true for the transition zone, However, ii 
the upper and lower flow regimes shear is closely related to velocity ani 
sediment concentration. 

Size of Bed Material.—The effects of the physical size of bed material o. 
the bed form are (1) its effect as a grain roughness, and (2) its influence o: 
the fall velocity or the effective fall diameter, which is a measure of the mo. 
bility of the particle. The mobility of the particle, when related to the flui 
forces, is of the utmost importance in determining whether the bed form i, 
ripples, dunes, plane bed, standing wave, or antidunes. 

In the second effect the physical size of the bed material, as measured b 
fall diameter!7 or sieve diameter, is important only as it is related to th 
other variables in Eq. 3. However, fall diameter has the advantage of elimi 
nating the shape factor and density of the particle as variables. Knowing th 
fall diameter, the fall velocity or effective fall diameter of the particle, i 
any fluid at any temperature, can be computed if the fluid’s density and vis 
cosity are known. 

The importance of size of bed material on the form of bed roughness wa 
indicated by the experiments using the 0.45 mm and the 0.28 mm sands. Th 
length of the dunes (crest to crest or trough to trough) observed while ex 
perimenting with the 0.28 mm sand were much longer than the dune length 
for the 0.45 mm sand. Although the amplitude of the dunes for the two sand 
were the same, their shapes were different. The dunes formed of 0.45 mr 
sand were more angular. The large spacing betweenthe dunes and the smalle 
separation zone and turbulence in the troughs accounts for the fact that som 
ripple runs exhibited a larger f value than the dune runs where.the 0.28 mr 
sand was used. The flow over the relatively long plane backs of these dune 
was quite efficient, reducing the over-all resistance to flow to a value les 
than that computed for the ripple case. The dunes, using the 0.45 mm sanc 
had much larger f values than the ripples in the same sand or the dunes i 
the 0.28 mm sand because of the closer spacing of the dunes, the large sepa 
ration zone in the dune troughs, and the greater magnitude of turbulence. Th 
increase in resistance to flow that occurred with an increase in sand size 


a a ee ee eee eee 

14 «application of Similarity Principles and Turbulence Research to Bed-Load Mov 
ment,” by I, A, Shields, 1936, Trans, from German by W. Poh and J, C. van Vehelin, 
S. Soil Conservation Service, Coop. Lab., Calif. Inst, of Tech., Pasadena, Calif, 

15 “The Flow of Cohesionless Grains in Fluids,” by R. A, Bagnold, Philos. Tran: 
Royal Soc, of London, Vol. 249, No. 964, 1956, pp. 235-297. x 

16 Discussion by H. A, Einstein and N, Chien of “Mechanics of Streams with Movak 
Beds of Fine Sand,” by N. H. Brooks, Transactions, ASCE, Vol, 123, 1958, pp. 553-5€ 

17 “Measurement and Analysis of Sediment Loads in Streams, Some Fundamentals 
Particle Size Analysis,” Inter-Agency Report 12, U. S. Govt, Printing Office, Washin 
ton, D. C., 1957, / 

is “The Effect of Fine Sediment on the Mechanics of Flow in Alluvial Channels,” 
W. L. Haushild, thesis presented to Colorado State Univ., at Fort Collins, Colo. in 191 
in partial fulfillment of the requirements for the degree of Master of Science. 
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. different from the trend reported by L. B. Leopold, F, ASCE, and T. 
Maddock, F. ASCE,19 

_ The change from lower to upper flow regime occurred at a much smaller 
| value, and the rate of change was faster for the finer sand, Also, using 
the fine sand the form of bed roughness changed from dunes to a plane bed 
when F = 0.4, and the bed remained plane as the boundary shear was in- 
eee until F = 1, then antidunes developed. Using the coarser sand, the 
lunes decreased in amplitude with increasing boundary shear when F =~ 0.6 
and transition conditions persisted until F ~ 1, then standing waves developed. 

That a plane bed developed with tranquil flow with the 0.28 mm sand and 
not with the 0.45 mm sand is attributed to the greater mobility of the finer 
sand and the limitations of the experimental equipment. The boundary shear 
required to eliminate the dunes was sufficiently small for the more mobile 
finer sand so that the change could occur at a depth and slope such that the 
F < 1. Whereas, with the coarser sand the discharge capacity of the flume 
system was not sufficient to eliminate the dunes at a depth and slope where 
the change would occur with F < 1. Field observations indicate that streams 
with 0.45 mm or coarser bed material, but having large depths and flat slopes 
may have a plane bed form of bed roughness when F < 1. 

Fall Velocity and Effective Fall Diameter.—As indicated in the preceding 
section, the fall velocity or the effective fall diameter of the bed material is 
9 major importance in determining the bed form, resistance to flow, and 
sediment transport. Because the size concept is so entrenched in the litera- 
ure, effective fall diameter is used in preference to fall velocity, Effective 
‘all diameter of a sediment particle is defined as the diameter of a sphere 
laving a specific gravity of 2.65 and a fall velocity in quiescent distilled 
water of infinite extent at a temperature of 24°C equal to the fall velocity of 
he particle falling alone in any quiescent stream fluid at stream temperature. 
This definition neglects the unknown effect of turbulence on fall velocity. 

The effect of the various factors given in Eq. 3 on fall velocity are well 
“nown although they are often overlooked. For instance, it is known that vis- 
‘osity and mass density of a fluid varies with temperature and affects the fall 
elocity, but the fact that fine sediment affects the viscosity18,20,21 and the 
mass density of the fluid are often overlooked. An increase in concentration 
if suspended fine sediment (wash load) or adecrease in temperature increases 
he viscosity and specific density of the liquid and decreases the fall velocity 
ind effective fall diameter of the bed material. Particles of different physical 
liameter have the same effective fall diameter if their fall velocities are the 
ame. With the same effective fall diameter the forms of bed roughness, re- 
stance to flow, and bed material load should be approximately the same in 
pite of the difference in physical size. 

_ With respect to temperature effect there have been conflicting statements 
n the literature. Vito A. Vanoni, F. ASCE, and Brooks22 reported that with 
n increase in temperature there is a decrease in resistance to flow; whereas, 
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19 “The Hydraulic Geometry of Stream Channels and Some Physiographic Implica- 
ions,” by L. B. Leopold and T. Maddock, Jr., Prof. Paper 252, U. S. Geol, Survey. 

20 “Fluidity and Plasticity,” by E.C, Bingham, Mc Graw-Hill Book Co, Inc., New York, 
922. 

21 «Flow Properties of Disperse Systems,” J. J. Hermans, ed., North-Holland Pub- 
shing Co,, Amsterdam, The Netherlands, : 
29 “Laboratory Studies of the Roughness and Suspended Load of Alluvial Streams,” 
y V. A. Vanoni andN.H. Brooks, Report No. E-68, Calif, Inst. of Tech., Pasadena, Calif, 
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D. W. Hubbell23 reported that with an increase in temperature there is an 
increase in the resistance to flow. These two apparently contradicting state- 
ments, that show the effect of temperature variation, are easily explained by 
considering the forms of bed roughness which existed in the alluvial channels 
in question, and the effect of temperature change on the fall velocity and, 
hence, effective fall diameter of the bed material particles. 

A decrease in temperature inereases the viscosity of the water and de- 
creases the fall velocity of the sand. Thus, a decrease in temperature de- 
creases the effective fall diameter of a given sand. Consequently, if a sand 
bed is covered with ripples and the temperature of the water is decreased, 
the mobility of the particles is increased due to the decrease in effective fall 
diameter of the sand, larger ripples form, and resistance to flow increases. 
On the other hand, if the form of bed roughness is near to or in the transition 
zone and there is a reduction in temperature, the reduction in effective fall 
diameter causes the dunes to shear out to a greater extent or perhaps causes 
the bed to become plane, which (ineither event) is accompanied by a decrease 
in resistance to flow. Both of these phenomena are reversible. These two 
variables (temperature and fine sediment) have been studied by the USGS. 

Seepage Force.—In an alluvial channel there is usually either inflow or 
outflow from the channel through the bank and bed material that causes seep- 
age forces because there is a loss of head. Although the effect of seepage 
forces on bed form has not been studied it is known that if there is inflow, 
the seepage force acts to reduce the effective weight of the sand, If there is 
outflow from the channel, the seepage force acts in the direction of flow and 
increases the effective weight of the sand. With inflow in the extreme case, 
the seepage forces can be large enough to set up a quick-sand condition and 
the weight of the material is in equilibrium with the seepage forces. The 
seepage force reduces or increases the effective size of the bed material in 
contact with the bed by changing the effective weight of the sand on the bed, 
Once a particle of sand loses contact with the bed, the seepage forces no longer 
have a direct effect on it. 

Gradation of the Bed Material.—The comments concerning the effect of the 
size of the material, fall velocity, and seepage force indicate that the bed 
material can be adequately described by a fictitious diameter. This obviously 
is not the case, Although, extensive studies into the effect of gradation of bed 
material on form of bed roughness, and resistance to flow have not been made, 
it would appear that an alluvial stream with a bed material composed of only 
one size would not react to fluid forces the same as if the bed material had 
the same median diameter but a large range of sizes. | 

Vanoni and Brooks22 investigated the effect of gradation of bed material 
on bed material load and concluded that for a given median diameter of bed 
material the load increases as the range of sizes or gradation increases, 
Einstein24 also recognized the importance of the gradation of bed material, 

Shape Factor of the Cross Section.—The effect of channel shape is rather 
nebulous. However, wide shallow channels have been simulated and under 
certain conditions two or more forms ofbed roughness have developed simul- 
taneously, for example, small dunes on part of the cross section and plane 
ee —ee—eeeeeeSSeeSESeeee Se 

23 “Investigations of some Sedimentation Characteristics of a Sand Bed Stream,” by 
D, W. Hubbel et al., Open-File Report, U. S, Geol. Survey, Lincoln, Nebr. 


24 «The Bed Load Function for Sediment Transportation in Open Channel Flows,” by 
H, A. Einstein, Tech, Bulletin 1026, U.S, Dept. of Agric., 1950, 70 pp. 
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bed, standing waves or antidunes on the other. The resistance to flow is vari- 
: able and difficult to evaluate whenthereis multiple roughness. The resistance 
to flow is dependent on the percent of the bed covered with each bed form. 
Multiple roughness and the variable resistance to flow associated with this 
condition has been observed in alluvial streams by the writers, 

| Current (1961) studies of resistance to flow and sediment transportation 
indicates that the forms of bed roughness for given alluvial bed material may 
‘be quite different in flumes of different size, other conditions being essen- 
tially the same. 

Using the large flume, the forms of bed roughness observed in the two 
regimes of flow were as illustrated in Fig, 2. These forms of bed roughness 
Seem to agree with field conditions. Using the 0.45 mm bed material in the 
small flume the observed bed configurations were different. When bed move- 
ment began, ripples did not form in the small flume. In fact, ripples never 
developed in this flume with this size of bed material. When bed shear in the 
small flume was increased sufficiently, the plane bed was replaced by dunes. 
These dunes were more two dimensional than those which occur in the large 
flume and these dunes never, at any time, had a ripple roughness superposed 
on them as was the case in the large flume. 

In the rapid-flow regime, conditions were quite similar in both flumes., 
However, the resistance to flow in the smallflume was much greater because 
flow was more nearly two dimensional which caused the total width of the 
flume to be occupied by the waves and antidunes that developed. 

The results of this comparison indicate the many problems that the ex- 
perimentor faces who tries to effectively utilize the data collected by various 
investigators from their flumes of different size and design. 


PREDICTION OF FORM OF BED ROUGHNESS 


Thus far, no completely adequate method of predicting form of bed rough- 
ness has been developed. M. L. Albertson, F, ASCE, and others presented25 
a method suitable for predicting bed forms for the laboratory case, but it 
was not completely satisfactory for the field case. Similarly R. J. Gardi2§ pre- 
sented a figure relating Shield’s parameter YDS/A 7d, the Froude number, 
F, and the form of bed roughness, This figure was not as satisfactory for the 
laboratory case but was superior for the field case. As another possibility 
Fig. 4 was developed which is a modification of the figure by Hsin-Kuan Liu, 
M. ASCE for ripple formation.27 Fig. 4 separates the bed forms for the vari- 
ous sand sizes for both streams and flumes, fairly well. Because the fall 
velocity is used, it also takes into account all viscous and gravitational forces 
that appear to affect the bed form. However, Fig. 4 may not separate the bed 
forms when larger depths are studied. Until additional data are available that 
cover an adequate range of depth, it will be difficult to develop the criterion 
necessary to predict form of bed roughness in alluvial channels with confi- 
dence. 


- 25 Discussion by M. L. Albertson, D.B. Simons, and E, V, Richardson of “Mechanics 
of Sediment-Ripple Formation,” by Hsin-Kuan Liu, Proceedings, ASCE, Vol. 84, No. HY 1, 


February, 1958. 
26 “Total Sediment Transport in Alluvial Channels,” by R.J. Gardi, thesis presented 


Colorado State Univ., at Fort Collins, Colo., in 1960, in partial fulfillment of the re- 
uirements for the degree of Doctor of Philosophy. { 

27 “Mechanics of Sediment Formation,” by Hsin-Kuan Siu, Proceedings, ASCE, Vol. 
, No. HY 5, May, 1957. 
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In the meantime, the possibility of predicting the bed configuration based 
upon the mechanics of surface waves and the appearance of the water surface 
should be considered. This method is not applicable to design but is ex- 
tremely useful for analysis. However, some training is required on the part 
of the individual before he can apply it effectively. 

Whether the flow is rapid or tranquil can be determined by observing the 
direction of travel of an artificially induced surface wave using the wave 
celarity concept. If the disturbance generated to oppose flow moves upstream 
with respect to a fixed point of observation, flow is tranquil, and if it is 
swept downstream, flow is rapid. 

It is also usually possible to determine form of bed roughness and whether 
the lower or upper flow regime exists in an alluvial channel by observing the 


FIG, 4.—FORMS OF BED ROUGHNESS IN ALLUVIAL CHANNELS 


condition of the water surface. In the upper flow regime, the water surface 
is likely to appear relatively smooth and fast or there will be some evidence 
of standing waves and antidunes [as illustrated in Figs. 2(e), 2(f), 2(g), ance 
2(h)] which immediately fix the general form of bed roughness that exists. 
In the lower flow regime the velocity will be relatively slow and the major 
rises in water surface (usually boils) will be out of phase with the dunes. For 
a rippled bed, the water surface will be quite plane and placid except for very 
small depths, In this case small waves willbe generated on the water surface 
by the sand ripples. With a dune bed configuration, there will be turbulence 
generated at the water surface in the form of boils downstream of the dunes 
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Usually the color of the water is different in the boil area due to the large 
_ Concentration of suspended sediment carried to the water surface within the 
boil. The strength of the boil activity is dependent on the magnitude, spacing 
and shape of the dunes, and the depth of flow. In the transition zone the water 
surface will be rather plane with minor boils appearing on the water surface. 
The classification of a stream into the proper regime of flow and form of 
bed roughness by observing the water surface combined with some knowledge 
of the bed material would give some information about the transport of sedi- 
ment and resistance to flow. This information could then be applied to the 
interpretation of stage-discharge relations and other flow phenomena, 


FIELD PROBLEMS 


The conditions observed in the field and in the large flume agree quite 
well, However, in the field the slope of energy grade line through a given reach 
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FIG. 5.—TYPICAL QUALITATIVE STAGE—DISCHARGE CURVES 
FOR ALLUVIAL CHANNELS 


of channel varies within narrow limits and the depth varies greatly as the 
discharge changes; but in the flume studies the slope can be varied widely 
and depth variations are small due to the physical limitations of the pumping 
plant and the flume. There is ample evidence that a change in depth can shift 
the form of bed roughness from dunes to transition, plane bed, and antidunes, 
The reduction in resistance to flow which occurs as dunes are eliminated 
accounts for the discontinuity in stage-discharge relations which have been 
‘observed in the laboratory and for many natural streams.10,11,12 Typical 
qualitative stage-discharge curves which illustrate the effect of clenge of bed 
roughness on stage are presented in Fig. 5, Einstein and Chien in their 
discussion of Brook’s paper3 recognized that the depth-discharge curves 
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illustrated in Fig, 5 could occur, However, contrary to present evidence, they 
thought these curves would only rarely if ever occur in nature, 

In Fig, 5(a) the break which is shown is caused by a change in bed rough- 
ness from dunes to plane bed or standing waves. The break occurs on the 
rising stage at a larger discharge than on the falling stage. There is usually 
considerable scatter in the stage discharge relation with ripples and dunes 
on the bed due to the wide variation of resistance to flow. Under certain cir- 
cumstances, this variation can be explained, For example, in Fig. 5(b) the 
magnitude of resistance to flow lags the actual discharge, that is, the change 
of bed roughness lags the change of discharge. This results in a smaller re- 
sistance to flow and smaller depth than would normally occur for equilibrium 
flow. The reverse occurs on the falling stage. In this instance, at the peak 
discharge, assume the form of bed roughness is large dunes. As discharge 
decreases, the large dunes are not reduced inamplitude and spacing at a rate 
in keeping with the reduction in discharge, and the resistance and depth are 
larger than for equilibrium flow on the falling leg of the hydrograph. This 
results in a loop type of stage-discharge curve which resembles a hysteresis 
curve. 

The loop curve of Fig. 5(b) may be reversed, Then the falling stage points 
would plot below the rising stage points or the two curves may coincide with 
the falling stage points plotting on the same curve as the rising stage points, 
or the rising and falling stage curves may cross, The form of the curve which 
results depends upon the rate of change of discharge with time, the magnitude 
of the discharge, the characteristics of the bed material, the effect of the fine 
sediment load, and the characteristics of the channel. Under certain circum- 
stances, temperature could play the same role as the fine sediment. 


SUMMARY 


The forms of bed roughness that exist in an alluvial channel because of 
similarities in form, resistance to flow, and sediment transport are divided 
into categories of lower flow regime, transition zone, and upper flow regime. 
Based upon laboratory and limited field investigations, the major forms of 
bed roughness in the normal order of occurrence with increasing slope with 
depth, bed material, temperature and concentration of fine material constant 
are as follows: 


Lower flow regime 


1. Ripples, 
2. Dunes with ripples superposed, and 
3, Dunes. 


Transition (Bed roughness ranges from dunes to plane bed or standing 
waves) 
Upper flow regime 


1. Plane bed, - 
2. Standing waves, and 
3. Antidunes. 


Resistance to flow and sediment transportation vary greatly as form of bed 
roughness changes. Consider the two sands with median diameters of 0.45 
mm and 0.28 mm which have been thoroughly investigated in a large recircu- 
lating flume, For a dune bed configuration f varies from 0.0489 to 0.1490, and 
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the corresponding bed material load ranges from 100 ppm to 1,200 ppm. For 
the antidune condition, f varies from 0.0247 to 0.0672 and bed material load 
ranges from 5,000 ppm to 42,000 ppm. 

The form of bed roughness depends primarily on the slope of the energy 
grade line, depth, fall velocity or effective fall diameter of the bed material 
and shape of the channel, In the transition zone the bed form also depends on 
the antecedent form of the bed. There are several other variables whose 
effect on the bed form may only be of minor importance. They are the seepage 
force, the gradation of the bed material, and the cohesive properties of any 
fine sediment in the bed, and the physical size of the sand as a grain rough- 
ness. 

It has been determined that fall velocity or effective fall diameter has a 
significant effect on the bed form. Sands with a large difference in physical 
size have been found to have the same bed form, resistance to flow, and 
sediment transport when their effective fall diameters were the same, The 

factors affecting fall velocity have an important effecton the bed form. These 
are the size of the bed material, viscosity, and mass density of the water- 
sediment mixture, mass density of the bed material, and the shape of the sand 
particles. 

Knowledge of the forms of bed roughness that develop in alluvial channels 
and the variables which influence bed roughness is essential to the theory of 
resistance to flow in alluvial channels, the design of stable channels and the 
refinement of sediment transport theories. 

To further improve existing concepts of flow in alluvial channels it is 
important to proceed with the following studies: 


The effect of depth of flow; 
the effect of size of bed material using larger sand sizes than 0.5 mm; 
the effect of gradation of bed material; 
the mechanics related to the meandering and braiding of streams; 
the effect of width-depth ratio; 
. the mechanics of flow around bends in alluvial channels to establish 
stability criterion and magnitude of the energy losses; and 

". the mechanics of diffusion, dispersion, and turbulence in alluvial 


channels. 


The development of new instruments is an essential part of future work, 
Developing methods of measuring turbulence, instantaneous velocity profiles 
and boundary shear in open channels is necessary for the continued develop- 
ment of the theory of flow in open channels, 
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ESTIMATING POTENTIAL EVAPOTRANSPIRATION@ 


By W. Russell Hamon,! A. M, ASCE 


SYNOPSIS 


Methods of computing potential evapotranspiration by analytical procedures 
have been based on the application of the turbulent-transport and energy- 
balance concepts. Empirical formulas, correlating some temperature function 
and adjusting for day-time hours, have proved valuable in practical utilization. 
This latter approach has been used to formulate a simple computational pro- 
cedure whereby average daily potential evapotranspiration is represented as 
proportional to the product of day-time hours squared and the saturated water 
vapor concentration (absolute humidity) at the mean temperature. The day- 
time factor was determined from a consideration of the disparity between net 
radiation and temperature, latitudinally, and the fact that transpiration is re- 
stricted during darkness since the leaf stomata are closed. 

Computed values of potential evapotranspiration obtained by the new pro- 
cedure are compared with those obtained by the more complex Thornthwaite 
method (24)2 and other methods currently employed, General applicability 
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seems justified from comparisons between observed and computed values of 
potential evapotranspiration, both on a yearly and seasonal basis. 


INTRODUCTION 


The determination of evapotranspiration is of professional interestto civil 
engineers, hydrologists, meteorologists, and agriculturists. To these and 
others, the problems of water and the conservationof water supplies—in view 
of the “population explosion” and the increasing demandon water resources— 
is more than academic. 

Since evapotranspiration represents the transport of water from the earth’s 
surface back to the atmosphere—the reverse of precipitation—it assumes an 
important role in the hydrologic cycle and mustbe given careful consideration 
in assessing the water resources of a particular region. Given estimates of 
evapotranspiration then, through an appropriate accounting system of soil 
moisture, it becomes possible to indicate the water balance of an area. 

Evapotranspiration is controlled by meteorological factors-radiation, wind, 
temperature, and humidity; soil moisture tension; soil moisture conditions; 
and plant physiological factors. Its direct measurement is difficult and climatic 
data have been used in a number of relations (2,7,20,24) to obtain useful esti- 
mates, These procedures rely either on the correlation between water evapo- 
rated from an open water surface and that from a vegetative area or an empir- 
ical relationship with climatic parameters and day-time hours, When soil 
moisture is not limiting, evapotranspiration is almost entirely dependent upon 
the amount of solar energy received by the surface and the resulting tempera- 
tures (27). Such a condition of optimum soil moisture results in “potential” 
evapotranspiration, 

The concept of potential evapotranspiration was introduced by; .C..We 
Thornthwaite (24) and is understood to represent the rate at which water 
would be removed from an extended area covered by an actively growing 
green crop, completely shading the soil, and with a non-limiting supply of 
water, The rate of evapotranspiration is not influenced by the advection of 
air that has been heated over drier areas, and the ratio of energy utilized in 
evaporation to that heating the air remains constant, M. A. Kohler, F. ASCE 
(13) has suggested that the evaporation from a large water surface must be 
equivalent to potential evapotranspiration. In this respect, the several methods 
adaptive to computing evaporation from an open water surface seem applicable 
in the estimation of potential evapotranspiration, provided the necessary data 
may be obtained, On the other hand, soil and plant factors may significantly 
alter the water lost by transpiration so that simpler methods based on tem- 
perature and day-time hours may be effective, The purpose of this paper is 
to formulate such a functional relationship. 


ANALYTICAL METHODS 


The turbulent-transfer approach to the determination of reservoir evapora- 
tion has received considerable attention. This method views evaporation as 
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equal to the flow of water vapor through the surface layer of air and deter- 
mined by the product of the vertical gradient of vapor pressure and the rate of 
Mixing of the air. The rate of mixing, expressed as a coefficient of diffusion 
or austauch coefficient, depends upon the rate of change with height of wind 
Speed. Such equations as derived (22), (19), (25) could be used for the deter- 
mination of evaporation over any kind of surface where the difference in sur- 
face properties is expressed through a roughness constant (22). The practical 
application of this approach is severely restricted due to the accuracy re- 
quirements in measurement of wind and humidity at two levels to determine 
the rate of change in wind and humidity with height (26). 

In the extensive Lake Hefner study (29) of lake evaporation, a simple 
empirical formula 


E = k (ey, * eq ) £(u) re aye Sosierte aabiheShtac weaken (1) 


relating evaporation, E, to the difference in vapor pressure-saturated at the 

surface temperature, ey, and the atmospheric vapor pressure, eg (as Dalton 
formulated in 1802); adjusted by a wind factor, f(u), and a proportionately con- 
stant, k, gave results comparable to those obtained by the complex mass- 
transfer formulas. Surface temperature is required for this simplified ap- 
proach so that it is of little use since pan evaporation is usually measured 
Whenever the water temperature is observed. The Lake Hefner study also 
confirmed that the energy-budget approach, when appropriate data are avail- 
able, may be used successfully in estimating lake evaporation. 

The problem in the energy balance approach to estimating evaporation is 
tO measure or assess all other sources and sinks of energy and to leave evap- 
Oration as the only unknown, The complete energy balance of a crop volume 
has been expressed by C. B. Tanner (23) but it may ordinarily be considered 
as consisting of net radiation, Ry; heat flux to the soil by conduction, S; heat 
flux to the air by convection, A; and the heat of vaporization of water, E; ex- 
pressed as 


BPO URL Ba ety cea 5. he as es Pucca py alee (2) 


Usually the soil heat flux is neglected but may be measured directly such as 
the net radiation while the sensible and latent heat terms can be separated 
using the Bowen ratio method (23). 

There has been considerable discussion of the validity of equating the 

constants of the heat and vapor transport equations in the Bowen ratio and 
apparently such an assumption is only valid for conditions over a water sur- 
face. In using the energy budget under conditions of restricted soil moisture, 
as forced in estimating evapotranspiration, Ray K. Linsley, F. ASCE (15) 
states that the Bowen ratio is no longer useful as a means of eliminating the 
sensible heat exchange from consideration, 
- In determining the evapotranspiration from an irrigated plot, Maurice H. 
Halstead (8) has illustrated that a combination of both the turbulent-transfer 
and energy approach is required, To evaluate the “oasis-effect,” Halstead 
recast turbulent-transfer equations so that they employed temperature and 
humidity differences between the surface and one height, and introduced 
another aerodynamic equation to indicate the height to which the influence of 
the oasis extends (8). Assuming conditions of net radiation, vertical wind 
gradient, roughness, temperature, and humidity such that the evaporation 
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from a 6-ft tank, when surrounded by a completely dry area, was 0.45 cm 
per hr; then, by increasing the size of the “oasis” to a 1-mile field, the com- 
puted evaporation rate was reduced to less than1/3 that of the 6-ft tank, This 
extreme example demonstrates that actual evaporation can exceed that re- 
sulting from net radiation by a considerable amount, depending upon the size 
of the field and the difference between moisture conditions of the field and its 
surroundings, Even for a field of cotton extending 10 miles upwind, observa- 
tions confirm that energy was gained from some drier area upwind (14). 

H. L. Penman (20) has successfully combined the energy budget and aero- 
dynamic approaches to compute open water evaporation. For the energy bal- 
ance estimate, Penman neglected net heat flux to the water and considered 
the Bowen ratio to remain constant; and by assuming identity in the heat 
transfer and vapor transfer coefficients in the aerodynamic relations, the 
need for knowing the surface temperature was eliminated, Conversion fac- 
tors have been applied to these hypothetical open water evaporation values 
to estimate what is termed potential evapotranspiration, obtained by com- 
paring tanks with turf and water; and for a 2-yr period of observation they 
were found to range from 0.6 during the winter to 0.8 in the summer, at 
Rothamsted, England. The dominant factor in determining the seasonal cycle 
of the conversion factor is the range in day-time hours which leads to bigger 
summer values and smaller winter values in more northerly latitudes but to 
only a small annual range at low latitudes (21). Such coefficients only convert 
pan evaporation to an equivalent area of well watered vegetation and, therefore, 
cannot represent potential evapotranspiration due to the “oasis” effect that 
increases with the aridity of a locality. Besides, the Bowen ratio takes on 
larger values in arid climates (23) and the empirical aerodynamic relation 
varies with climate, Even so, important generalizations about potential evapo- 
transpiration still hold, as stated by Penman (21): 


1. “For complete crop covers of different plants having about the 
same colour, i.e., the same reflection coefficient, the potential trans- 
piration rate is the same, irrespective of plant or soil type.” (Physi- 
ological responses of the plant, particularly at high temperatures, may 
alter this response.) 

2, “This potential transpiration rate is determined by prevailing 
weather.” 

3. “The transpiration of a short-green cover cannot exceed the 
evaporation from an open water surface exposed to the same weather,” 


EMPIRICAL RELATIONS 


A number of empirical methods have been advanced to estimate actual and 
potential evapotranspiration, all employing some temperature function. Early 
work by Hedke (11) produced an equation wherein the seasonal estimated con- 
sumptive use (evapotranspiration) was considered proportional to the “effec-. 
tive” heat expressed in day degrees and taken as the difference between the 
monthly temperature and germinating temperature for each crop. Using yearly 
data in the same fashion, Robert L. Lowry and Arthur F, Johnson, F, ASCE. 
(17) found a high correlation between consumptive use, under conditions where 


aye moisture was not limiting, and accumulated degree days with a base of 
32 -F, 


a 
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H. F. Blaney, F. ASCE, and W. D. Criddle, M. ASCE (2) correlated actual 
measurement of consumptive use with the product of monthly mean tempera- 
ture and the monthly day-time hours asa percentage of yearly day-time hours, 
to obtain applicable coefficients for different irrigated cropping conditions in 
the West. The average coefficient varies from around 0.6 for field crops to 
1.2 for rice over the growth period under arid and semi-arid conditions, The 
coefficients, however, show a marked seasonal variation and for grapefruit 
in the vicinity of Phoenix, Ariz., they range from 0.55 in January to 0.75 in 
August (2). 

A widely used method for estimating potential evapotranspiration (P.E.) is 
that devised by Thornthwaite (24) through correlating meanmonthly tempera- 
ture with evapotranspiration as determined from the water balance for valleys 
where sufficient moisture was available to maintain active transpiration, The 
equation, when obtaining 30-day month values in centimeters, takes the form 


P.E, = 1.6 L (20¢) rage) ped Roa wt. eae a (3) 
J 

in which L represents day-time hours expressed in units of 12 hr; t is the 

mean monthly temperature in degree centigrade; J, a heat index obtained 

by the summation of the twelve monthly indices (i), in which 


ee AAS DR OE, Aare Senin ers (4) 
and the value of a is obtained from a cubic function of the heat index (J), as 
a=0.000000675 J3 - 0.0000771 Tot 001100.) 0-4900 39 ope (5) 


In commenting on the formula, Thornthwaite (24) stated that its work-ability 
requires the use of nomograms and that the chief obstacle in developing a 
rational equation is the lack of understanding of the reason why potential 
evapotranspiration corresponding to a given temperature is not the same 
everywhere. 


DEVELOPMENT OF METHOD 


In seeking any new relationship to represent average potential evapo- 
transpiration as a function of temperature, due consideration must be given 
to the geographical variation between temperature and radiation. Average 
daily temparatures in July over the eastern United States vary from 80-85° in 
the southern states to 65-70° over the northern tier states. In contrast, the 
average incoming solar radiation exhibits only a small variation (6) and for 
clear skies the incoming radiation is a constant quantity at the summer sol- 


-stice for latitudes 25° to 50° over North America (10). 


Since net radiation comprises the heat source of an area, aside from ad- 
vected heat in certain localities, a relationship between temperature and net 
radiation is desired. Hamon (9), upon examination of limited net radiation ob- 
servations and adjusting observed incoming radiation by computing total long- 
wave nocturnal radiation, found that for monthly data net radiation could be 
represented as a linear function of the product of day-time hours and tempera- 


ture, as 
Re DSP ye ee ee (6) 
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in which R, is the average daily net radiation in langleys (calories per square 
centimeter); D isthe average duration of day-time hours in units of 12 hr; and 
T is the average daily mean temperature in degrees fahrenheit. 

The need for considering wind in average estimates of potential evapo- 
transpiration may be reckoned from the effect of wind on open water evapora- 
tion and the response of transpiring plants. For daily wind movement of 50 to 
100 miles, the percent error in evaporation with percent change in wind is in 
the order of 0.25% in dry regions (16). Penman states (21), in reference to his 
formula (developed in a humid region), that it is possible to have 100% over- 
estimate in evaporation by the aerodynamic term leading to only 10% over- 
estimate in open water evaporation. Perhaps more conclusive in substantiating 
the elimination of wind as a variable under natural conditions is the response 
of plants. Y. Yamaoka (31), in a laboratory study of transpiration from plants 
exposed to constant temperature and relative humidity, found that the trans- 
piration rate remains constant under a particular radiation intensity for wind 
speeds above 1 m per sec (about 2 miles per hr), 

It then appears justifiable to consider wind as a constant factor in any 
climatological estimate of average potential evapotranspiration. For a condi- 
tion, however, where an “oasis” situation exists any estimate of actual evapo- 
transpiration would be highly dependent on advected heat and, therefore, the 
wind. 

Adhering to the concept of potential evapotranspiration and accepting a 
constant Bowen ratio, we may assume that anaverage relative humidity exists 
at some height near the surface which would result in a proportionaility 
between the saturated vapor density at the mean surface temperature, and the 
average difference between the saturated vapor density at the surface and the 
vapor density at the reference height. Therefore, other things being equal, and 
referring to the turbulent transport equations, average potential evapotran- 
Spiration may be considered as proportional to the saturated water vapor 
density P; at the air temperature near the surface, represented as 


Pier Rei e'fG gels on Rae kee (7) 


Considering that transpiration is a daytime phenomenon due to the function- 
ing of the leaf stomata (open only in daylight) and that turbulence occurs prin- 
cipally during the daytime, an adjustment for day-time hours is required, 
Also, the adjustment enters again to account for the disparity between net 
radiation and mean temperature, latitudinally. 

From the preceeding considerations, it seems feasible to formulate the 
final simplified expression for potential evapotranspiration as 


PIB APG D* Pes adnne pad stacy ak orig sane 


in which P. E. represents the average potential evapotranspiration in inches 
per day; D is the possible hours of sunshine in units of 12 hr; P; is the satu- 
rated water vapor density (absolute humidity) at the daily mean temperature 
in grams per cubic meter; times 10-2; C = 0.55, chosen to give appropriate 
yearly values of potential evapotranspiration as indicated by observations re- 
ported in the literature (24), (17). 


EVALUATION 


Yearly values of potential evapotranspiration as obtained from the new 
relation (Eq. 8) for which C = 0.55, are practically identical to those obtained ' 
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by the highly complicated Thornthwaite method (Fig. 1) with a consistent de- 
parture of about 5% noted only for latitudes of southern Florida, A further 
evaluation on yearly data has been made by comparing average computed 
values of potential evapotranspiration to consumptive use data for valleys 
covering a wide variationin climate as obtained by Lowry and Johnson (17) and 
G, R. Williams, F. ASCE (30), as listed in Table 1 and shown in Fig. 2. 
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FIG. 1.—COMPARISON OF COMPUTED AVERAGE ANNUAL 
POTENTIAL EVAPOTRANSPIRATION 


lation (Eq. 8) is whether 

An important test of the reliability of the new re 

_ a aimplitied procedure is adequate in representing the seasonal varia- 
. Several comparisons are made, 

: ieee ares values of water loss from an irrigated valley, Mesilla 


Vally in New Mexico, are reported inthe study of Lowry and Johnson (17). The 
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FIG, 2.—CONSUMPTIVE USE FOR BASINS COMPARED WITH 
COMPUTED POTENTIAL EVAPOTRANSPIRATION 


TABLE 1.—DATA OF OBSERVED WATER LOSS AND COMPUTED AVERAGE 
POTENTIAL EVAPOTRANSPIRATION FOR FIG. 2 


Length of 
Record, in 
years 


(3) 


Valley and State Mean Water| Estimated 


(1) 


1, Mesilla, N. Mex, 109,000 
2, Pecas, N. Mex, 37,850 36.3 
3, Sangamon River, Ill. 1,640,000 29.3 
4, Green River, Ky. 5,000,000 31.4 
5. Tallapoosa River, Ala.-Ga, 1,060,000 34.0 
6. Mad River, Ohio 307,000 27.3 
7, Skunk River, Iowa 1,890,000 27.6 
8. N. Fork of White River, Mo. 755,000 Sabai! 
9. N. Platte, Wyo.-Neb. 462,000 24,1 
10. Black River, Wisc. 494,000 23.9 
11, Cypress Creek, Tex. 545,000 38.350 ue 
12, Wagon Wheel Gap, Colo. 222 15.0 
13, Mich,-Il,, River, Colo. 43,000 17.2 
14, West River, Vt. 23.9 


Lake Cochituate, Mass. 
. Swift River, Mass, 
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correspondence of these observed values of actual water loss to computed 
average potential evapotranspiration is satisfactory as seen in Fig, 3. 

Estimates of the average water loss for periods of adequate water supply 
has been obtained by W. F. Fox, F. ASCE (5) for basins in Wisconsin by the 


NY 
. COMPUTED AVERAGE 
SAE ae N POTENTIAL EVAPO. 
N TRANSPIRATION 
ats 
E NN 
= CONSUMPTIVE USE N N 
S _ NN 
" NNSA 
Y NNSA 
= YNNNA 
: NNN A 
2 NNNA 
NNN A 
N NWN A 
NNN A 
> Db LD > 
LR ReMis tacM Da scpeaAn es 0. ND 
MONTHS 
FIG, 3.—SEASONAL DISTRIBUTION OF CONSUMPTIVE 
USE AND POTENTIAL EVAPOTRANSPIRATION 
Ta a aT ee cl ee 
‘e & eB fRaT FORD, WISCONSIN) 
(e—*) TREMPEALEAU RIVER 
20 a (DODGE, WISCONSIN) 


as (»—x) BAD RIVER 
18 = (MELLEN, WISCONSIN) 
16 —— 


=F COMPUTED NORMAL POTENTIAL 
EVAPOTRANSPIRATION 
a (LaCROSSE, WISCONSIN ) 


OTRANSPIRATION, INCHES. 


(AVERAGE DAILY EVAP 
ao 


04 
0 O 2 — = 
eT ———— = toa 
5 1 15 15 i 5 1 5 5 5 
JAN. FEB. MAR, APR. MAY JUNE JULY AUG. SEPT. oct. Nov. DEC. 


G. 4.-POTENTIAL EVAPOTRANSPIRATION FROM WATER BALANCE FOR BASINS | 


er-balance method. His determined average values of water loss and the 
mputed curve of average potential evapotranspiration compare favorably 


shown in Fig. 4. 
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Both Penman and Blaney (4) have computed estimates of evapotranspira- 
tion by their respective methods for Asheville, N. C., and for the climate of 
this location should approximate potential evapotranspiration. In addition, 
computations by the Thornthwaite method and by Eq. 8 have been listed in 
Table 2. Results from the new procedure are fully satisfactory in comparison 
with the other widely used methods and succeeds in correcting a short-coming 
of the Thornthwaite method wherein winter-time estimates are considered too 
low and mid-summer values excessive (18). 

An extensive investigation has been made of the yearly water loss and 
lake evaporation in New England and New York (12). The difference between 
the yearly lake evaporation in this study, obtained by an adaptation of the 
Penman equation (16), and values of average potential evapotranspiration, 
obtained by Eq. 8, show small departures as seen in Fig. 5. This close corre- 


TABLE 2.—COMPUTED AVERAGE MONTHLY EVAPOTRANSPIRATION FOR 
ASHEVILLE, N. C., BY SEVERAL METHODS 


(Totals) 


Methods: Blaney (a), Penman (b), Thornthwaite (c), and by Eq. 8 (d). 


obtained by Eq. 8, show small departures as seen in Fig. 5. This close corre- 
spondence is only expected inhumid regions where the “oasis” effect is small. 
Since Eq. 8 has shown every indication of representing a unique and simple 
procedure for obtaining useful estimates of average potential evapotranspira- 
tion, computations have been made using monthly mean data for the regior 
eek ae Rockies. A map of generalized average yearly values is presentec 
The conditions required for potential evapotranspiration eliminates the 
advection of air heated over surrounding drier areas and allows for equilibrium 
conditions where evaporation is proportional to the energy available. Evapora- 
tion losses from pans or soil tanks are always in excess of potential evapotran- 
spiration—ranging from small differences in humid regions to excesses 0! 
"100% in hot-arid regions. For large irrigated areas or valleys this discrepancy 
is practically eliminated. Under non-irrigated conditions, the actual evapo- 
transpiration depends upon the available soil moisture and the potential evapo- 
transpiration may be adjusted for its estimation. 

The decrease in evapotranspiration in relation to soil moisture is viewec 
differently by various investigators. Thornthwaite (18) has adopted a systen 
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that decreases evapotranspiration on a direct function of the remaining soil 
moisture on the basis of field data that showed the fraction of available energy 
used in evaporation as proportional to soil moisture. From data of C. H. 
Wadleigh (28), the suppression of plant growth andthe soil moisture above the 
permanent wilting point are seento vary inthe same fashion with soil moisture 
stress, This indicates that the adjustment to potential evapotranspiration to 
obtain actual evapotranspiration should be made proportional to the plant 
available moisture under conditions where soil moisture is limiting, Whether 
or not soil moisture is “limiting” depends upon the soil and soil aeration, 
the plant, and meteorological factors. 


CONCLUSIONS 


The equation developed to estimate average potential evapotranspiration 
(Eq. 8) adheres to the desirable features of requiring only readily available 
, mean temperature data-converted to saturation water vapor density and day- 
time hours, expressed in an easily computable fashion. General applicability 
seems justified in view of the close correspondence between observed and 
computed values of average potential evapotranspiration, both ona yearly and 
seasonal basis for widely scattered localities. 


APPENDIX,—READING REFERENCES 


1, “Monthly Consumptive Use Requirements for Irrigated Crops,” by H. F. 
Blaney, Proceedings, ASCE, Vol. 85, No, IR1, March, 1959. 


2. “Determining Water Requirements in Irrigated Areas from Climatological 
and Irrigation Data,” by H. F, Blaney and W. D, Criddle, U.S. Dept. of 
Agric., Soil Conservation Service, Washington, D.C., 1950. 


3. “Evaporation and Evapotranspiration Investigations in the San Francisco 
Bay Area,” by H. F, Blaney andD, C. Muckel, Journal, Amer, Geophysical 
Union, Vol. 36, No. 5, October, 1955, pp. 813-820. 


4. Discussion by H, F. Blaney of “Estimating Evaporation,” by H. L. Penman, 
Transactions, Amer, Geophysical Union, Vol. 37, No. 1, February, 1956, 
pp. 46-48. 


5. “Computation of Potential and Actual Evapotranspiration,” by W. E. Fox, 
U.S, Weather Bur., Washington, D.C., 1946. ng j 


“Average Solar Radiation in the United States,” by S. Fritz and T. H. Mac- 
Donald, Heating and Ventilating, July, 1949, 


7. “The Derivation of an Equation for Potential Evapotranspiration,” by M. H. 
Halstead, Johns Hopkins Univ. Pub. in Climatology, Vol. 4, No. 5, Sea- 
brook, N. J., 1951, pp. 10-12. 


= 


HY 3 EVAPOTRANSPIRATION 119 


10. 


1. 


12. 


13. 


14, 


15. 


16. 


L7. 


18. 


19, 


. “Some Meteorological Aspects of Evapotranspiration,” by M. H. Halstead 


and W. Covey, Proceedings, Soil Sci, Soc. of America, Vol. 21, No. 5, 
September-October, 1957, pp. 461-468. 


. “Estimating Potential Evapotranspiration,” by W. R. Hamon, Unpublished 


Thesis, Massachusetts Inst. of Tech., Cambridge, February, 1960. 


“Insolation as an Empirical Function of Daily Sunshine Duration,” by W. 
R. Hamon, L. L. Weiss, and T. W. Wilson, Monthly Weather Review, Wash- 
ington, D.C., June, 1954. 


“Consumptive Use of Water in Irrigation,” Progress Report of the Duty of 
Water Committee of the Irrigation Division, Transactions, ASCE, Vol, 94, 
1930, pp. 1349-1399. 


“Average Annual Runoff and Precipatation in New England-New York 
Area,” by C. E. Knox and T. J. Nordenson, U.S, Geol. Survey, Hydro. 
Invest. Atlas HA7, Dept. of Interior, Washington, D.C. 


“Meteorological Aspects of Evaporation Phenomena,” by M. A. Kohler, 
Transactions, Internatl. Assoc. of Science and Hydrology, Vol. II, pp. 
421-436. 


“Some Aspects of the Relationship of Soil, Plant and Meteorological Fac- 
tors to Evapotranspiration,” by R. R. Lemon, A. H, Glaser, and L. E. 
Salterwhite, Proceedings, Soil Sci. Soc. of America, Vol. 21, No. 5, 
September, October, 1957, pp. 461-468. 


“The Common Ground of Meteorology and Hydrology,” by R. K. Linsley, 
Bulletin, Amer. Meteorological Soc., Vol. 41, No, 8, August, 1960. 


“Hydrology for Engineers,” by R. K. Linsley, M. A. Kohler, and J. L. H. 
Paulhus, McGraw-Hill Book Co., New York, 1958. 


“Consumptive Use of Water in Agriculture,” by R. L. Lowry and A. F. 
Johnson, Transactions, ASCE, Vol. 107, 1942. 

“The Measurement of Potential Evapotranspiration,” by J. R. Mather, ed., 
Publ, in Climatology, Vol. VII, No. 1, Johns Hopkins Univ., Lab. of Climat- 
ology, Seabrook, N. J., 1954. 

“Some Further Consideration of the Measurement and Indirect Evaluation 
of Natural Evaporation,” by F. Pasquill, Q. J. Royal Meteorological Soc., 
Vol, 76, 1950, pp. 287-301. 


_ “Natural Evaporation from Open Water, Bare Soil and Grass,” by H. L. 


Penman, Proceedings, Royal Soc. of London, Series A, Vol. 193, 1948, 
pp 120-148. 


. “Evaporation: An Introductory Survey,” by H. L. Penman, Netherlands 


Journal Agric. Sci., Vol. 4, No. 1, 1956, pp. 9-29. 


. “Eddy Diffusion of Momentum, Water Vapor and Heat near the Ground,” by 


N, E. Rider, Philosophical Transactions, Royal Soc. of London, Series A, 
1954, pp. 481-501. 


. “Energy Balance Approach to Evapotranspiration from Crops,” by.CaB: 


Tanner, Proceedings, Soil Sci. Soc. of America, Vol. 24, No. 1, January- 
February, 1960. ‘ 


120 May, 1961 HY 3 


24, “An Approach toward a Rational Classification of Climate,” by C. W. 
Thornthwaite, Geographical Review, Vol. 38, 1948, pp. 55-94. 


25. “The Determination of Evaporation from Land and Water Surfaces,” by 
C. W. Thornthwaite and B, Bolzman, Monthly Weather Review, Vol. 67, 
1939, pp. 4-11. 


26. “Note on the Variation of Wind: with Heightin the Layer near the Ground,” 
by C. W. Thornthwaite and M. H. Halstead, Transactions, Amer. Geophysical 
Union, 1942, pp. 249-255. 


27. “The Water Budget and its Use in Irrigation,” by C. W. Thornthwaite and 
J. R. Mather, Water, the Yearbook of Agriculture, U. S. Dept. of Agric., 
Washington, D.C., 1955, pp. 346-357. 


28, “Soil Moisture in Relation to Plant Growth,” by C. H. Wadleigh, Water, 
the Yearbook of Agriculture, U.S. Dept. of Agric., Washington, D.C., 1955, 
pp. 358-361, 


29. “Water Loss Investigations: Lake Hefner Studies,” Tech. Report, U.S. 
Geol, Survey, Professional Paper No. 269, 1954. 


30. “Natural Water Loss in Selected Drainage Basins,” by G. R. Williams, 
et al., U.S. Geol. Survey, Water Supply Paper 846, U.S. Dept. of Interior, 
Washington, D.C., 1940. 


31. “Experimental Studies on the Relation between Transpiration Data and 
Meteorological Elements,” by Y. Yamaoka, Transactions, Amer. Geo- 
physical Union, Vol. 39, No. 2, 1958, pp. 249-266. 


2823 May, 1961 HY 3 


Journal of the 


HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


ROUGHNESS SPACING IN RIGID OPEN CHANNELS 
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Maurice L, Albertson,2 F, ASCE 


SYNOPSIS 


A series of experiments to determine the effect of roughness spacing on 
open channel flow is reported. The experiments were conducted in an 8-ft 
wide by 72-ft-long adjustable-slope flume. Roughness elements consisting of 
sheet metal baffles were placed on the bedof the flume at various longitudinal 
and transverse spacings. 

The data are analyzed in terms of the von Karman- Prandtl concepts of 
turbulent flow near a rough boundary. Logarithmic flowformulas are developed 
in which the roughness factor is relatedto the size and spacing characteristics 
of the roughness elements. 

The experimental results are compared with those of some other investi- 
gators. In addition, some observations are made onthe von Karman turbulence 
coefficient and the problem of depth determination in rough channels. 


INTRODUCTION 


Notation.—The letter symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, in the Appendix, 


_Note.—Discussion open until October 1, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 3, May, 1961. 
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Many important phenomena involved in the flow of water in open channels 
are only partially understood. One of the most basic of these is the effect of 
boundary roughness on the flow. 

Over the years, a number of empirical roughness formulas have been 
developed. Almost all are based on the Chezy formula 
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in which V is the mean velocity of flow, C denotes the Chezy resistance func- 
tion, R defines the hydraulic radius of the channel cross section, and S is the 
slope of energy gradient. Formulas developed by Manning, Ganguillet and 
Kutter, Bazin, and others relate C tothe boundary roughness. These formulas 
have little or no rational foundation, They give fair results when applied over 
the fairly narrow range of flow conditions on which they were based, but fre- 
quently lead to results that are seriously inerror when applied outside of this 
range. Experiments by Sayre%, A. R. Robinson* M. ASCE, and Sayre and 
Albertson®, show conclusively that the roughness coefficients inthese formulas 
do not necessarily remain constant for a particular boundary roughness. 

In the last few decades, the development of the boundary layer theory has 
opened the way toward a rational understanding of some of the flow phenomena 
involving boundary roughness, 

Using the von Karman-Prandtl concepts of turbulent flow near a rough 
boundary in addition to data from his own experiments with pipes having sand- 
coated interiors, J. Nikuradse& developed an artificial roughness standard for 
pipes, that is based on the diameter of the sand grains kg constituting the 
roughness, 

The application by G. H. Keulegan? and others of the Nikuradse roughness 
standard to open channels has been quite successful in describing grain-type 
roughness in wide, open channels. However, it has been found inadequate for 
describing certain other types of roughness in which the relative spacing in 
addition to the relative size of the roughness elements is an important bound- 
ary characteristic. Dune and ripple formations onthe beds of alluvial channels 
are considered typical of the latter type of roughness. 

Using square strips extending across the bed of the channel as roughness, 
Ralph W. Powell8,9, F. ASCE, studied the effect of the longitudinal spacing of 
the strips. Robinson10 and Robinson and Albertson4 reported experiments in 
which the size of geometrically similar roughness baffles was varied, but the 


3 “Artificial Roughness Patterns in Open Channels,” by W. W. Sayre, thesis present- 
ed to Colorado State University, at Fort Collins, Colo., in 1957, in partial fulfillment of 
the requirements for the degree of Master of Science, 

4 “ Artificial Roughness Standard for Open Channels,” by A. R. Robinson and M, L. 
Albertson, Transactions, Amer. Geophysical Union, Vol. 33, No. 6, December, 1952, pp. 
881-888, 

5 “The Effect of Roughness Spacing in Rigid Open Channels,” by W. W. Sayre and M. 
L. Albertson, Dept. of Civ. Engrg., Colorado State Univ., Fort Collins, Colo., CER59 
WWS31, October, 1959, 78 pp. | 

6 “Strémungsgesetze in rauhen Réhren,” by J. Nikuradse, Forschungsarbeiten auf 
dem Gebiete des Ingenieur Wesens, Heft 361, Vol. 4, 1933, Berlin, pp. 1-22. 

“Laws of Turbulent Flow in Open Channels,” by G.H.Keulegan, Journal of Research, 
U.S, Natl. Bur. of Standards, Vol. 21, December, 1938, pp. 707-741. 

8 “Flow in a Channel of Definite Roughness,” by R. W. Powell, Transactions, ASCE, 
Vol, 111, 1946, pp. 531-566, 2 

9 “Resistance to Flow in Rough Channels,” by R. W. Powell, Transactions, Amer. 
Geophysical Union, Vol. 31, No. 4, August, 1950, pp. 575-582. 

10 “Artificial Roughness in Open Channels,” by A. R. Robinson, thesis presented to 
Colorado State Univ., at Fort Collins, Colo., in 1950, in partial fulfillment of the require- 
ments for the degree of Master of Science. 
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ratios of longitudinal and transverse spacings to baffle height were held con- 
stant. For a particular roughness pattern they demonstrated that the Chezy 
resistance function depends only on the relative roughness (ratio of flow depth 
to baffle height), assuming rough boundary conditions. 

These and other studies have contributed toward the establishment of a 
reproducible artificial roughness standard that will be applicable to open 
channels characterized by a “dune”-type roughness, Itis hoped that this study 
will stimulate further interest toward this end. 


ANALYSIS 


Theoretical Analysis.—The von Karman-Prandtl equation for velocity dis- 
tribution near rough boundaries is of the form 


Vv 2.30 

= = =— log i + Cy cee cueLahels te ceeeoe eee (2) 
= 
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in which v is the velocity at distance y from the boundary, Kk defines the von 
Karman turbulence coefficient, y is the distance above the bed of the channel, 
k denotes the measure of roughness size, C; defines a constant of integration, 
andy/to /p denotes the shear velocity at theboundary. Its derivation appears 
in many standard fluid mechanics texts. It has been verified by experiment 
over a wide range of conditions in pipes, in fluids moving past a flat plate, and 
in open channel flow provided that secondary circulation is not an important 
consideration. 

An expression for the mean velocity can be obtained by integrating Eq. 2 
over the depth of flow in the following manner 
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that results in the expression 
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in which yp is the normal depth (the depth of flow occurring when slopes of 


energy gradient, water surface, and bed are equal), and Cg defines an experi- 
mental constant dependent on the roughness spacing. From the Chezy and 


Darcy-Weisbach equations 
V 
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in which C is the Chezy coefficient, R denotes the hydraulic radius, S defines 
the slope of the energy gradient and f is the Parcy-Weisbach resistance co- 
efficient. Eq. 4 can thus be written in the form of a genéral resistance equa-~ 


tion such as 
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or 


= ee 


0.81 Yn 
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Equations proposed by Keulegan, Powell, and Robinson and Albertson are of 
the form of Eq. 6(a). The von Karman-Prandtl equation for rough pipes is of 
the form of Eq. 6(b). 

Because k in Eq. 3 is defined in terms of roughness size, it is logical to 
conclude that Cg must be a factor dependent upon the arrangement and shape 
of the roughness elements, An alternate form of Eq. 6(a) would be 

y 
EC naieelll- loge A. WS ra iet eae vent ty ee (6c) 
tag 
in which X appears as a roughness parameter dependent on the size, shape 
and spacing of the roughness elements. In other words, X should completely 
describe the boundary roughness, 

Two fundamental assumptions are implicit in the foregoing analysis. The 
first is that the boundary is hydrodynamically rough (roughness protruding 
above the laminar sublayer) and that consequently viscosity has negligible 
effect. The second is that the channel is of sufficient width, or that the bed 
roughness is so great relative to the side-wall roughness or both, as to elimi- 
nate any appreciable sidewall effect. 

Dimensional Analysis.—The variables used in the dimensional analysis 
were chosen so as to represent the experimental conditions. They can be 
grouped under three headings. 


Boundary Characteristics.—The boundary characteristics are as follows: 


width of channel; 

height of roughness baffle; 

width of roughness baffle; 

longitudinal spacing of roughness baffles; 

transverse spacing between roughness baffles; 

thickness of roughness baffle; 

shape factor of roughness elements; 

roughness parameter describing size, spacing and shape of rough- 


ness elements, and having the dimension of length; and 
dX = shape of channel. 
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Flow Characteristics.—The flow characteristics are: 


V = mean velocity of flow; 
R = hydraulic radius; and 
To = boundary shear, 


Fluid Characteristics.—The following comprise the fluid characteristics; 


p = mass density of fluid; 
p. = dynamic viscosity of fluid; and . 
Ay = difference in specific weight of water and air at the fluid interface. 


If X completely defines the roughness as indicated in the theoretical analy- 
sis, the relationship existing among the variables describing the roughness 
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can be stated 


The general relationship among the remaining variables is 
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Choosing V, R, and /p as repeating variables, the Buckingham Pi-theorem 
yields six7 -terms among which the following relationship exists: 
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In a wide channel, y, canbe substitutedfor R, and if there are no side-wall 
effects B/R can be eliminated. Because the shape of the channel remains 


rectangular and wide throughout the experiments, can also be eliminated. 
Eq. 9 can now be restated as 
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in which R and F are the Reynolds and Froude numbers respectively. 
By use of 
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Eq. 10 can be simplified and rewritten as 
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If the height of the roughness elements is large relative to the thickness of 
the laminar sub-layer, viscous effects are negligible, and R can be eliminated. 
The Froude number assumes importance only when appreciable surface waves 
or disturbances are present. Consequently, when the boundary is hydrodynam- 
ically rough and whenthere are no surface eaves, Eq. 12 reduces to the simple 
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Turning now to the problem of roughness geometry, Eq. 7 stated that 
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e ratio b/a and the shape factor of the roughness elements 0 were held con- 
t and, therefore, canbe eliminated. Because tis very small (approximately 
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0.05 in.) relative to all other roughness dimensions, t/a can be considered in 
effect equal to zero, and hence of negligible importance. Because the dimen- 
sions of the individual roughness baffles remain constant the ratio X/a is 
independent of baffle size, andcanbe replacedby y , a parameter that appears 
to be dependent only on the baffle spacing. Eq. 7 can now be written in a sim- 
pler form 


y= %- oy = a gi 6 A aoa peaks ee (14) 


By reasoning, the roughness size and spacing dimensions can be combined 
into a single dimensionless density parameter a b/x(e+b). This parameter is 
numerically equal to the ratio between the combined area of all the roughness 
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FIG, 1.—FLUME AND CIRCULATION SYSTEM 
elements projected perpendicularly to the direction of flow as the numerator 
and the total area of the flume bed as the denominator, 


ab  _ 2 roughness area 
x(e+b) ~ area of flume bed ‘*°°cccc ccc: 15, 


Another expression for the spacing parameter is then 
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DESCRIPTION OF EXPERIMENTS 


The experiments were conducted in an 8-ft wide by 2-ft deep rectangulai 


flume which measured 72 ft in length. The flume and circulation system art 
shown in schematic form in Fig. 1. 
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The bed and walls of the flume were constructed of 1/2-in. plywood rein- 
forced by 2-in. by 4-in. ribs and cross members. The flume rested on a pair 
of 6-in. I-beams which were supported by screw jacks at 10-ft intervals. 
Changes in the longitudinal slope of the flume could be accomplished by ad- 
justing the jacks. 

At the flume entrance flow was distributed and excessive local turbulence 
eliminated by a manifold diffuser. 12 

The water-surface level in the flume could be controlled by adjusting the 
vertical-slat tailgate at the downstream end of the flume. The tailbox, a 30 
hp single-stage propeller pump, and a 14-in. pipeline completed the major 
elements of the circulation system. 

The discharge was controlled by a valve in the 14-in. line. Discharge was 
measured by means of a 10 1/2-in. orifice plate located in the 14-in. pipe, 24 
ft downstream from the pump. 

Water-surface slope and flow depth were determined by five stationary 
point gages which were spaced over a length of 40 ft along the flume, The 
point gages were read with the aid of an engineer’s level. An 11/32-in dia- 
meter Prandtl-type pitot tube with a tilting manometer was used to obtain 
velocity profiles in the vertical, The water temperature was measured with 
an ordinary laboratory mercury thermometer. 

Sheet metal baffles measuring 6-in wide and 1 1/2-in high were used as 
roughness elements. As shown in Fig. 2, the baffles were placed at definite 
longitudinal and transverse intervals, thus forming symmetrical repeating 
patterns. Each baffle was centered on the space between baffles in the rows 
immediately upstream and downstream, hence alternate rows were the same. 
All patterns were symmetrical about the centerline of the flume, Half-width 
baffles were placed adjacent to the side walls in order to maintain the same 
geometry there as along the centerline. The photograph in Fig. 3 shows one 
of the roughness patterns. All of the patterns are shown in schematic form in 
Fig. 4. 

In performing the experiments, the discharge, slope, and roughness pattern 
were the independent variables, and the normal flow depth was the primary 
dependent variable. Once apredetermined combination of independent variables 
was established, the procedure during any given run consisted essentially of 
manipulating the tail-gate setting until the water-surface slope equaled the 
slope of the flume bed. This could be done quite accurately by intentionally 
establishing and then interpolating between slight M, and Mg backwater curves. 
The condition of equal slopes of the water Surface and the bed defined normal 
flow. The depth of normal flow was then measured, and velocity profiles were 
Obtained in at least three verticals across the flume. 

Tests were conducted at discharges varying from 2 cfs to 6 cfs and at 
slopes of 0.001 and 0.003 for each of the six roughness patterns. A total of 
twenty-seven runs was made. 


PRESENTATION AND DISCUSSION OF RESULTS 


_ Data describing the overall flow characteristics and the roughness are listed 
in Table 1. The relationships developed in the theoretical and dimensional 
analysis are largely substantiated by the experimental results in the form of 


12 The Manifold Stilling Basin,” by G. R. Fiala and M, L. Albertson, Dept, of Civ. 
on Colorado State Univ, Research Foundation, October, 1958, 32 pp. 
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dimensionless graphs that are presented and discussed in the following sec- 
tions. In addition, information pertaining to the von Karman turbulence coef- 
ficient, and to the problem of depth determination in rough channels is also 
presented. 

Chezy Coefficient as a Function of Relative Roughness.—Fig. 5 is a plot 
of the Chezy resistance function against relative roughness based on the baffle 
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FIG. 2,—DEFINITION SKETCH OF ROUGHNESS BAFFLE INSTALLATION 


height. The data fall into a family of parallel lines when plotted on semi- 
logarithmic paper. Each line represents data obtained from experiments with 
a particular roughness pattern. Robinson’s datal0 which are superimposed on 
the plot, represent experiments in which the same roughness pattern, but two 
different baffle sizes, were used. The roughness baffles used by Robinson 
were geometrically similar (b= 4a) to those used by the writers. 
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Data obtained by H.J. Koloseus,!1M. ASCE, for conditions of stable flow 
are also plotted on Fig. 5. These tests were conducted in a 2-ft-wide flume 
with 3/16-in. brass cubes as roughness elements. Itis interesting to note that 
Koloseus’ data were obtained in the super-critical (F>1) flow range, whereas 
Robinson’s and the writers’ data all represent flow conditions well within the 
sub-critical range. The data are all described by the equation 

Beenie tee er Cy a eae ty ee (17) 
VE a 
which is similar to Eq. 6(a) given in the theoretical analysis. The constant 
C2 for each line, that is each roughness pattern, is defined by extrapolation 
as the value of C/./g at which y,/a = 1. The writers’ and Koloseus’ data show 


FIG, 3.-EXAMPLE OF A ROUGHNESS PATTERN 


clearly that Cg is a function of the longitudinal and transverse roughness 
spacing, and Robinson’s data show that itis independent of absolute roughness 
size. ‘ 

+ Fig. 6isa plot of the Chezy resistance function against the relative rough- 
ness based on the roughness parameterX . Description of the boundary rough- 
ness by the single parameter X groups all the data about a single curve 
P 

} 
- 11 «The Effect of Free-Surface Instability on Channel Resistance,” by H. J. Koloseus, 
th sis presented to the State Univ. of Iowa, at Iowa City, Iowa, in 1958, in partial fulfill- 
ment of the requirements for the degree of Doctor of Philosophy. 
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FIG, 4.—SCHEMATIC DIAGRAM OF THE SIX ROUGHNESS PATTERNS 
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FIG. 5.—VARIATION OF RESISTANCE FUNCTION WITH RELATIVE ROUGHNESS Yn 
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described by the equation 
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that is similar to Eq. 6(c) in the theoretical analysis. The parameter X has 
dimension of length and isa constantfora particular boundary roughness. The 
value of y, extrapolated to C/,/g = Odefines xX for a particular combination of 
baffle size and spacing. The properties of a and Cg in Eq. 17 are combined in 
X according to the relationship 


ieee o0CHGE eo 7 beteiooh «5 (19) 
The ratio X/a was suggested as a baffle spacing parameter ? in the dimension- 
al analysis. 

Eqs. 17 and 18 are essentially equivalent. However, definition of the relative 
roughness in terms of X as in Eq. 18 appears to have a number of advantages. 
Grouping all of the data about a single curve is the most obvious advantage. 
Description of the roughness by a single term greatly facilitates the inclusion 
of the roughness effect in various dimensionless parameters. In addition, 
Eq. 18 possesses the same advantages in simplicity of application as does the 
Manning formula in which the single termndescribes the boundary roughness, 

Should it be demonstrated, however, that the coefficient 6.06 varies appreci- 
ably with roughness form or channel shape but not with roughness spacing. 
Eq. 17 would seem to be more significant than Eq. 18. 

The empirically-determined coefficient 6.06 in Eqs. 17and18 is equivalent 
to a von Kdrman turbulence coefficient of 0.38. This compares favorably with 
the coefficient 0.40 determined by Nikuradse® for turbulent flow in smooth and 
rough pipes, and the coefficients determined by a number of investigators in 
open channels. 

In Figs. 5 and 6 some of the points representing observed flow conditions 
for the two sparsest of the writers’ roughness patterns are seen to fall below 
the curve established by the other points. Inthese runs a combination of small 
depths and relatively high velocities occurring directly over the crests of the 
roughness baffles caused the formation of breaking standing waves immediately 
downstream from the baffles. In effect, these baffle-generated waves con- 
stituted energy sinks not accounted for by either Eq. 17 or Eq. 18. This ex- 
plains the low C/./g values obtained for these runs, and also agrees with 
Koloseus’ results!1 for unstable flow which clearly demonstrate that channel 
resistance is increased in a systematic manner by free-surface instability. 

The tendency of Robinson’s data to drift below the curve at larger values 
of relative roughness is attributed to an increase in sidewall effect accompany- 
ing the increase in depth. The flume usedby Robinson was only 9 in. wide, and 
the depth of flow in several runs approached the same magnitude. 

Comparison of Chezy and Manning Formulas.—Fig. 7 compares Kq, 17 
with the Manning formula. Basically it is the same plot as Fig. 5, except that 
it is extended over two additional log cycles of relative roughness, The straight 
dashed lines are extrapolations of the experimental data, each line represent- 
ing a particular degree of boundary roughness. The solid curved lines repre- 
sent constant values of Manning’s n thatwere plotted onthe C//g versus y,/a 


a 
ee) The line Cg = 10 compares Eq. 
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0.125 | 0.186 
0.125 | 0.120 


0.125 |.0869 
0.125 |.0637 
0.125 |.0198 


Robinson 


O47 |.0390/0.937 
-0833 | .0780/0.937 
7 -0156 |.00186 | 0.119 
15 |.0156 |000574|.0368 
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FIG, 6.—VARIATION OF RESISTANCE FUNCTION WITH RELATIVE ROUGHNESS Yp/% 
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FIG. 7.—COMPARISON BETWEEN RESISTANCE BASED ON MANNING FORMULA AND ON LOGARITHMIC EQUATION 
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17 with the Manning formula overarange of roughness commonly encounterec 
in practice, 

Fig. 7 demonstrates conclusively that application of the Manning formule 
to the range of roughness included in the experiments would lead to serious 
error. It is seen that the boundary described by the baffle spacing ratios 
x/a = 2 and e/a = 3.9 wouldhave values of Manning’s n ranging from approxi- 
mately 0.060 to 0.100 over the relatively small two-fold range of depth covere¢ 
by the experiments. However, in the region where the Manning n curves are 
approximately tangent to lines having a slope of log 6.06, the Manning formulz 
and Eq. 17 should give results which are in fairly close accord. This is the 
range of conditions within which the Manning formula is commonly applied tc 
open channel flow. 

Influence of Roughness Pattern.— 

Resistance to Flow.—Fig. 8 is a plotofthe roughness spacing parameter 7 
against the roughness density ratio Sa 5 2. that was suggested in the dimen- 
sional analysis. Data from several additional sources are superimposed on the 
plot. Robinson’s and Koloseus’ data have been discussedin the preceding sec- 
tion. H. Schlichting’s datal3 were obtained in a closed rectangular conduit 
measuring 4-cm high by 17-cm wide. Baffles measuring 0.3-cm high by 0.8-cm 
wide were fixed to the top plate. Powell’s8 1946 roughness consisted of 1/4-in 
and 1/8-in square strips extending continuously across the bottom and up the 
sides of an 8-in-wide flume. W. Rand’s datal4 were obtained in a 2-ft-wide 
flume in which 1/4-in, square strips extending across the bottom were usec 
as roughness, Bazin’s experiments!5 were conducted in a 2-m-wide flume ir 
which strips 1-cm thick and 2.7-cm wide extended across the bottom and ur 
the sides. In the Waterways Experiment Station experiments16 triangulaz 
ridges 0.04-ft high with a vertical upstream face and the downstream face 
sloped at 2 to 1 served as roughness. The ridges extended across the bottom 
and up the sides of the 2.08-ft-wide flume, The brushed-concrete surface upor 
which the ridges were superimposed was rougher than the base surface in the 
other experiments. 

Considering the differences in roughness and channel forms it seems that the 
various curves should be compared only in a qualitative sense. There is < 
general tendency in Fig. 8 for the spacing parameter (¥ = X/a), that is propor- 
tional to the hydraulic roughness by a factor of roughness height, to increase 
with increasing roughness density. For the continuous strip roughness 7 


reaches a maximum at about _2.)_ = 9. 
ju ian) 0.1 after which the boundary in effect 

becomes smoother as the roughness density increases. For the triangulaz 
ridges, maximum roughness occurs at about _2.)__-9.9. . 
dges, ug ju ee) 0.2. Evidently, for com. 


parable spacings, the ridges are somewhat rougher than the square strips 
particularly at close spacings. At the more distant spacings the comparitively 


13 “Experimental Investigation of the Problem of Surface Roughness,” by H, Schlicht 
ing, Memorandum No, 823, Natl. Advisory Committee for Aeronautics, Tech. 1937, 

14 Unpublished Data, by W. Rand, State Univ. of Iowa, Iowa City, Iowa, 1952, 

15 “Recherches Hydrauliques: Memoirs Divers Savants,” by H.E, Bazin and H. Darcy 
Science, Mathematiques, et Physiques, Vol. 19, Paris, France, 1865, 

16 “Roughness Standards for Hydraulic Models Report No, 1, Study of Finite Bound 
ary Roughness in Rectangular Flumes,” Waterways Experiment Sta., U. S, Corps c 
Engrs., Vicksburg, Miss., Tech. Memorandum No, 2-364, 1953, 29 pp. 
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rough brushed-concrete surface that served as abase for the triangular ridges 
evidently came into play. 

For the continuous roughness, it is interesting to note that the roughness 
density parameter reduces to a/x, because the transverse spacing e between 
roughness elements is zero. It seems likely that the curves representing the 


Robinson— baffles, /=4 
Schlichting—baffles,b/4=2.67 

Koloseus— cubes, =! 
Powell — continuous square strips, bho 
Rand —continuous square strips, bio 


Bazin—contnuous wide strips, %4—>@ 
U SWE S—continuous triangular ridges, D/g—>@ 
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FIG. 8. VARIATION OF SPACING PARAMETER WITH ROUGHNESS DENSITY 


isolated roughness elements (baffles and cubes) also should reach maximum 
values of 7 , but at greater roughness densities that extend beyond the range of 


the data. 


Fe Ke = 
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As would be expected, the rising limbs of the curves in Fig. 8 are arranged 
in order of increasing relative width (b/a) of the roughness elements as the 
spacing parameter wy increases. 

A plot somewhat similar to Fig. 8 which contains more complete data for 
continuous-strip roughness has been given by Joe W. Johnson, M. ASCE, and 
E. A. LeRoux, 17 


FIG, 9.—DIMENSIONLESS VELOCITY PROFILES 


Mimi... SSE eee 
17 Discussion by J. W. Johnson and E. A, LeRoux of “Flow in a Channel of Definite 
Roughness,” by R. Ww. Rowell, Transactions, ASCE, Vol. 11, 1946, pp. 555-559, 
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Velocity Distribution.—Profiles of velocity distribution in the vertical are 
plotted in dimensionless form on Fig. 9, Each point represents a weighted 
average of at least three readings taken at the same elevation above the bed, 
but at different verticals across the width of the flume. The weights were 
assigned according to the position of the vertical with respect to the rough- 
ness elements. 

The nature of the velocity profiles Cleasly seems to have been affected by 
roughness pattern, particularly by the longitudinal spacing of the baffles. 
Velocity profiles for the four denser roughnesspatterns are quite well repre- 
sented by the equation 


Vv 
Oo 
T— 
p 
‘or which the constants were determined empirically. 
For the sake of comparison the von Karm4n-Prandtl equation 


= 5.75 log we SOG ak co tatets 8. fin ng ee (21) 


= 6.06 log ¥ roe ere oe ee Be Oe (20) 


72 


or velocity distribution in turbulent flow near a rough boundary is plotted as a 
oroken line. The Nikuradse sand-grain roughness kg is related to X by the non- 
inear relationship 

12.2 X 1.053 


% ahead Shee cane (22) 
n 
hat corresponds approximately to kg/X= 11 over the range of conditions 
ested. 

The deviation of the points representing velocity readings taken for the two 
jparsest roughness patterns is attributed to the non-uniformity of the rough- 
less and hence the turbulence in the sparse-pattern experiments. The raw 
elocity-profile data are available.9° 

General Resistance Diagrams for Open Channels.—Fig.10and Fig. 11 were 
lotted for the purpose of determining the order and range of hydrodynamic 
‘Oughness covered in the experiments. Some available open-channel data rep- 
esenting other regions of the roughness spectrum are included in these plots 
Ss a means of testing the applicability of the Prandtl-von Kdrman-Nikuradse 
ype resistance diagrams, that were developed for flow in pipes, to flow in 
ride, rigid-boundary open channels. Assumption of wide-channel conditions, 
nat is y, = R, is justified for all data shown on Fig. 10 and Fig. 11 with the 
xception of the smooth-boundary data in which case R was substituted for yp. 

Fig. 10, which is an open-channel adaptation of the well-known general 
esistance diagram for flow in pipes, illustrates the effect of the Reynolds 
umber, that is viscosity, on the resistance function, 

_ The laminar-flow equation 

SSE By 8 0 cere a ia i atas eeeoe C5) 


3 derived directly from the following equations describing oo uniform, 
aminar flow 
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and see pe waite ae ee ( 24b) 


The smooth-boundary formula 


ge ae. SO Baye as 8" Gee. 


is a direct adaptation of Blasius’ empirical relationship for flow in smooth 
pipes. The other smooth-boundary equation 


LCs) be oie ee (26) 
ve Cc 
vgs 
is a modified form of Keulegan’s equation” 
See ae See ed eee See (27) 


forflowin a smooth channel of infinite width. The coefficient 5.75 was change: 

to 6.06 for the sake of conformance with Eq. 17 and Eq. 18, and the constan 

was changed to compensate for the alteration of the coefficient. 
Experimental data for the six roughness patterns fall well within th 


Ryf 
rough-boundary region that begins at ¥, 7x 733. At this condition the rati 
n 


X/5' of the roughness to the thickness of the laminar sublayer is one (see Fig 
11). The data are described by afamily of horizontal lines, each having a con 
stant value of relative roughness yy /y. Hence, in the region covered by th 
data, the resistance function is dependent only on relative roughness and i 
therefore independent of the Reynolds number or viscous effects. 

The laminar-flow and smooth-boundary data shown onthe plot were obtaine 
from a variety of sources. Laminar-flow data were obtained by W. M. Owen! 
in a rectangular, 1.50-ft-wide brass-plate flume, and by the Waterways Ex 
periment Station19 in a rectangular, 2.30-ft-wide flume witha smooth cemer 
finish. Smooth-boundary data were obtained by H. J. Tracy, F. ASCE, an 
C. M. Lester,20A. M. ASCE, in a rectangular, 3.50-ft-wide steel flume andb 
the Waterways Experiment Station29 in the previously described cemer 
flume. The transition from laminar to turbulent flowfor the smooth-boundar 
data occurs in the range 700<R < 1,200. 

Fig. 11 is an open-channel adaptation of the transition function for pipe: 


By means of the abscissa a that is proportional to the ratio X/6" , the dat 
n/X 


are arranged in order of hydrodynamic roughness. 

18 «laminar to Turbulent Flow in a Wide Open Channel,” by W. M, Owen, Tran 
actions, ASCE, Vol, 119, 1954, p. 1157. 

19 éStudy of River Bed Materials and their Movement with Special Reference to t 
Lower Mississippi River,” Waterways Experiment Sta., U.S. Corps of Engrs., VioksGu 
Miss., Tech. Memorandum No, 17, 1935, 161 pp. 

20 «Experiments in a Smooth, Rectangular Open Channel,” by H. J. Tracy and C. 
Lester, Administrative Report, ue S. Geol. Survey, 1959. 
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The smooth-boundary equation 


Cc Nhe R 
- 6.06 log he 6.06 log yw SEERA Ec Rae (28) 
g 


Ve =o 


ve X 


is obtained merely by subtracting 6.06 log y, /x from both sides of Eq, 26, 
The transition equation 


Sati) B. 
Yn x 

ee OG log <= = -6.06 log | 1+ ALB Sh thei (29) 
g 1.58 R 


was obtained by combining the smooth and rough-boundary equations in the 
‘Manner proposed by C. F. Colebrook and White?! for pipes, and is asymptotic 
to both the smooth-boundary and rough-boundary equations. The broken-line 
curves, indicating the extremes ofthe transition region, approximate the posi- 
tions of the uniform and non-uniform sand roughness curves on the trnasition 
plot for pipes. 

Data for the six roughness patterns appear at the rough end of the scale. 
Although most of the data representing Koloseus’ two patterns also fall in the 
rough-boundary region, a few points representing the sparser of the patterns 
lap over slightly into the transition region, The Waterways Experiment Station 
datal9 for beds consisting of 4 mm gravel, 0.525 mm sand, 0.310 mm. sand 
and smooth cement, and the steel-flume data of Tracy and Lester20 are also 
plotted on Fig. 11. The alluvial-bed data all represent flow conditions before 
the beginning of movement, so that inallcases the boundary is in effect rigid. 
In plotting the smooth-boundary data, it was assumed that X = 0.00004 ft for 
the smooth-cement finish, and that X = 0.00001 for the steel surface. Eq. 22 
was used in computing X-values from the median diameters of the 0.525-mm 
sand (X = 0.00019 ft) and the 0.310-mm sand (xX = 0.00011 ft). The 4-mm gravel 
data fall mostly in the rough-boundary region, overlapping into the transition 
region. The data for the 0,525-mm sand fall in the transition region, for the 
0.310-mm sand in the transition and smooth regions, and for the smooth- 
cement and steel-flume surfaces, along the smooth-boundary curve. 

. Seattering of the data in the smooth-boundary and transition regions is not 
as great as it appears to be when it is remembered that the ordinate C/ye - 
6.06 log y,/x is a direct difference rather than a percentile deviation index. 
For example, considering a hypothetical point located at C/ye - 6.06 log 


Yn/X = -3 and ae = 0.320, the deviation on the ordinate scale from the 
n 

smooth-boundary curve is one unit or 25%. However, based upon a C/s/g value 

of 15, which is about average for the data in the smooth and transition regions, 

deviation by one unit amounts to less than 7%. 

_ The ratios X/5' and k,/6' are plotted along the abscissa scale of Fig. 11. 

The ratio kg/X = 10 is an average over the approximate range of relative 


H SL 
_ 21 «Turbulent Flow in Pipes, with Particular Reference to the Transition Region Be- 
tween the Smooth and Rough Pipe Laws,” by C. F, Colebrook, Journal, Inst. of Civ. Engrs., 
Vol, 11, 1939, pp. 133-156. 
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roughness 10,000 > y,/xX >1 included in the data plotted on Fig. 11. The fact 
that the rough-boundary region begins at approximately X /6' = 1 appears to 
be significant, This compares with k,/65' = 6 indicated22 by Hunter Rouse F. 
ASCE, as the approximate demarcation between the transition and rough- 
boundary regions for pipes. 

The wide band within which data may fall in the transitionregion can be 
explained on the basis of roughness spacing or size distribution or both. 
Where the spacing between roughness elements is large or the size of rough- 
ness elements is not uniform, isolated elements even atX/5' << 1 may project 
above the laminar sublayer in sufficient numbers to upset the turbulence 
pattern generated by the smooth boundary, but not in sufficient numbers to 
cause rough-boundary turbulence. However, when the roughness elements are 
densely-packed and uniform in size, smooth-boundary flow may persist at 
X/65' approaching one and the transition from smooth to rough boundary con- 
ditions would be much more abrupt. 

The von Karman Coefficient.—One of the purposes of the study was to throw 
some additional light on the question of how the von Karman k might vary, if 
in fact it does vary, with boundary roughness as suggested by Rand.23 

The coefficient kK initially appeared in the literature as a “universal con- 
stant” of turbulent exchange in von Karman’s hypothesis24 of turbulent simili- 
tude. On the basis of Nikuradse’s work® with smooth and rough pipes, the 
value of K has commonly been accepted as equal to 0.40, although Vito A. 
Vanoni25 F, ASCE, in analyzing Nikuradse’s datafound K-values ranging from 
0.32 to 0.42 with 0.37 as anapproximate average. Evidence has been mounting 
steadily to show that K, as expressed in the logarithmic velocity distribution 
equations, is not a constant, but is systematically influenced by such factors 
as suspended-sediment concentration, secondary circulation, and perhaps 
boundary roughness, 

In the experiments. reported herein two methods of determining K were 
used. In the first method K was determined from analysis of the velocity- 
profile data through use of the relationship. 


See Sper Gi ea ot 2:8? 98. eee 


in whichV70/p = 8 yn S and Av is the change in velocity over one logarith- 
mic cycle of depth. The portions of the velocity profiles that followed most 
closely a logarithmic distribution were used in determining x, and velocity 
measurements close to either the baffle crests or the water surface were given 
only secondary consideration. 

The second method of determing K was based on Eq. 6a from which it is 
seen that the slope of each of the lines drawn in Fig. 5 is defined by the quantity 


a a a ee Be eS eee 
22 «Elementary Mechanics of Fluids,” by H. Rouse, John Wiley and Sons, Inc., New 
York, 1946, 376 pp. ; ; 
23 Discussion by W. Rand of “Wind-Tunnel Studies of Sediment Movement,” Proceed- 
ings, Fifth Hydr, Conf., State Univ. of Iowa, Bulletin 34, 1953, pp. 132-133, 3 
4 “Mechanische Ahnlichkeit und Turbulenz,” by Theodor von Karman, Proceedings, 
Third Internatl, Congress for Applied Mechanics, Stockholm, Sweden, 1930, 
; 25 “Some Effects of Suspended Sediment on Flow Characteristics,” by V. A. Vanoni, 


Proce ngs, Fifth Hydr. Conf., State Univ. of Iowa, Studies in Engrg., Bulletin 34, 1953, 
pp. ~1os, 
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2.30/K. Best-fit lines were drawn through the data representing each rough- 
ness pattern in Fig. 5 and an average or integrated K-value for each rough- 
ness pattern was obtained from the relationship 


in which AC/+/g is the change in the resistance function over one logarithmic 
eycle of relative roughness. 

In Fig. 12 K is shown as a function of roughness density. Values of K as 
determined by both methods are shown on this plot. The curves compare the 
average K-values as determined by each method, each point on the curve 
representing an average value of K determinedfrom all the runs for a particu- 
lar roughness pattern. 

Fig. 12 indicates that the roughness pattern is instrumental in causing the 
apparent variation in x . The principal significance of Fig. 12, however, lies 
in the comparison between the average Kk -values as determined by the 
velocity-profile method and the integrated method. Agreement is remarkably 
good for the three denser patterns, but for the two sparsest patterns there is 
little semblance of agreement. It is probable that the integrated method pro- 
vides a more reliable index of x than does the velocity-profile method, and 
that the large apparent xk -values obtained by the velocity-profile method for 
the sparse patterns are indicative of large-scale local vortices set up by the 
isolated roughness baffles, rather than such large variations in x. This is to 
be expected, however, because the kK - theory was originally developed for 
cases of wake-interference flow in pipes, where the turbulence, close to the 
boundary, is essentially random in nature, Extension of this theory to the 
isolated-roughness case, where the flow is three dimensional or the rough- 
ness is non-uniform, or both, is, therefore, neither logical nor justified. 

This being the case, the portion of Fig. 12 drawn as a solid line appears to 
justify the adoption of an average value of x for the range of conditions tested. 
In effect, selection of the constant 6.06 as the coefficient in Eqs. 17 and 18 
amounted to adopting the average value K = 0.38. 

It cannot be said that the experimental results provided an entirely satis- 
factory: answer to whether the von Karman x is in fact a constant for clear 
water in open-channel flow, or how and why x might vary if it is not a con- 
stant. For the type of roughness studied, the velocity distribution method of 
determining x proved deficient inanumber of respects for the sparse patterns. 
Without taking a large number of velocity profiles in verticals at various 
positions relative to the roughness baffles, it was impossible to ascertain 
whether or not a truly representative picture of the velocity distribution had 
been obtained. Limitations in the velocity-measuring equipment and deviations 
from the logarithmic law close to the baffles and near the water surface added 
to the complications. Nevertheless, based onthe experimental results, indica- 
tions are that, provided the conditions of essentially wake-interference flow 
and uniform roughness are met, x is a constant for steady, uniform flow of 
clear water. 

- Depth Determination in Rough Channels.—Various investigators, such as, 
Rand 6 have suggested that when roughness elements are closely spaced, the 
¢ : 3 


A i ici hannels,” 
_ 26 Discussion by W. Rand of “Artificial Roughness Standards for Open C ; 
PE ransactions, Amer, Geophysical Union, Vol. 35, No. 4, August, 1954, pp, 649-650. 
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LEGEND 

° -— Value of A for a single run. 
Determined from velocity profiles. 

e@ - Average value of A fora 
roughness pattern. Determined from 
velocity profiles. 

@® - Average value of A for a 


roughness pattern. Determined from 
x/O = 18 
e/a= 1.8 


CANg vs Yn/a_ relationship. 
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FIG, 12.—VARIATION OF VON KARMAN COEFFICIENT WITH ROUGHNESS DENSITY 
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jepth of flow should be measured from the water surface to some reference 
evel between the crests and the troughs of the roughness elements, 

As implied in the section describing the experiments, the depth of flow as 
measured from the water surface to the flume bed was used in the analysis. 
However, quite an extensive study was made of the velocity profiles obtained 
‘or the three densest roughness patterns before it was decided not to apply 
any depth corrections. 

A detailed description of the depth-correction analysis is beyond the scope 
af this paper. The negative results of this analysis appear justified for quite 
jense roughness patterns, however, when the roughness elements consist of 
hin isolated baffles such as usedin these experiments. This type of roughness 
zyeometry permits flow through the transverse spaces between roughness ele- 
ments. In addition, because the proportion of the total bed area actually taken 
ip by the roughness elements— <Sch) ~ 0.008 for the densest pattern—is 
sxtremely small, the turbulence down close to the bed level of the flume is 
orobably not greatly inhibited. The results of turbulence flume studies by 
Albertson2? reinforce these hypotheses. In these studies a staggered pattern 
of individual roughness baffles proved extremely effective in maintaining 
arge sediment concentrations in suspension without appreciable deposit, 
whereas extensive deposits occurred at comparable concentrations when 
he roughness consisted of baffles extending continuously across the width 
0% the flume. These studies have application here because of the basic 
similarity between the mechanisms of mass transfer andmomentum transfer. 


CONCLUSIONS 


Analysis of the experimental data from this and other studies gives rise 
0 the following conclusions with respect to evaluating the roughness of rigid- 
Oundary, rough open channels: 


1. The experimental results give additional evidence that the variation of 
he Chezy resistance function with the relative roughness is logarithmic in 
lature. 

2. The equations 


y 
= 6.06 log =~ 4: Cg sie dite, ees (TT) 


rae Sia 


y 
= 6.06 log esse cece eee eee (18) 


re found to apply over the range of roughness conditions encountered in the 
ix arrangements of roughness baffles tested. In these equations, the constant 
yy) depends only onthe relative spacing of roughness elements of height a, 
ndthe general roughness parameter X is a function of both the relative size 
nd relative spacing of the baffles. 

3. Eq. 17 is shown to be considerably more accurate than the Manning 
ormula over the range of roughness and flow tested. 


we 
27 «Turbulence Flume to Measure Bed Load,” by M. L, Albertson, Transactions, 
mer, Geophysical Union, Vol, 32, No. 6, December, 1951, pp. 909-911. 
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4, The ratio ay which is equivalent to the combined area of all rough- 


ness baffles projected perpendicularly to the direction of flow divided by the 
total area of the bed, is found to be an adequate definition of the roughness 
density. 

5, Distribution of velocity in the vertical direction when the roughness 
pattern is relatively dense is found to be similar to that described by the von 
Karman-Prandtl equations for rough boundaries. The equation 


Vv 


26.00 log 294" 26 BUSES Des (20) 

> x 

o/p 
adequately describes the velocity distribution in such cases, A pattern con- 
sisting of sparsely-spaced, isolated roughness elements, however, gives rise 
to persistent horizontal and vertical vortices, and consequently considerable 
local variations in the turbulence pattern which cause marked deviations from 
the von Kdrmdn-Prandtl type distribution. 

6. The general resistance diagram, in which the resistance function is 
plotted against the Reynolds number, and the Colebrook-White type transitior 
function are found to be applicable toproblems of uniform flow in wide, rigid- 
boundary open channels. 

7, Indications are that the von Karman turbulence coefficient x is inde- 
pendent of roughness pattern, provided that conditions of wake-interference 
flow and uniform roughness prevail. In the final analysis xk is in effect con- 
sidered a constant, equal to 0.38, for the range of conditions tested in this 
study. 

8. Negative depth corrections are found to be unnecessary for roughness 
systems consisting of staggered individual baffles, at least when the trans- 


abt 
verse spacing ratio e/b 2 about 1, and the volume ratio sS(e+p 8 small 
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APPENDIX.—NOTATION 


The following symbols, adopted for use in this paper, conform essential] 
with “American Standard Letter Symbols for Hydraulics” (ASA Z10.2-1942) 
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repared by a committee of the American Standards Association with Society 
epresentation, and approved by the Association in 1942: 
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Symbol 
a 
B 
b 


Cy 


Q 
he) 


Pune 02 re = 
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imCmciee< ce mn Aa ty 6) 


< 
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Definition Dimensions 


Height of roughness baffle. 
Width of flume. 

Width of roughness baffle. 
Chezy discharge coefficient. 
Constant of integration. 


Experimental constant dependent on rough- 
ness spacing. 


Transverse spacing between roughness baffles. 
Froude number. 

Darcy-Weisbach resistance coefficient. 
Acceleration of gravity. 

Measure of roughness size. 

Equivalent sand roughness (Nikuradse). 
Roughness coefficient (Manning). 

Flow discharge. 

Hydraulic radius. 

Reynolds number. 

Slope of energy gradient. 

Thickness of roughness baffle. 

Mean velocity of flow. 

Velocity at distance y from boundary. 
Longitudinal spacing of roughness baffles. 
Distance above bed of channel. 


Normal depth, that is, depth of flow occurring 
when slopes of energy gradient, water surface, 
and bed are equal. 


Difference in specific weight of the two fluids 
involved (water and air) at fluid interface. 


Thickness of laminar sublayer. 


11.6v 
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von Karman turbulence coefficient. 


L1/2/T 


L 
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Channel shape factor. 

Dynamic viscosity of fluid. 

Kinematic viscosity of fluid. 

Fluid mass density. 

Shape factor of roughness elements. 
Shear intensity. 

Shear intensity acting on bed of channel. 
Roughness parameter. 


Spacing parameter. 
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EDDY DIFFUSION IN RESERVOIRS AND PIPELINES 


By Frank L. Parker,! A. M. ASCE 


SYNOPSIS 


Previous studies have shown that the ocean eddy diffusion coefficient is a 
function of eddy size. Reservoir tests with radioactive tracers confirm this. 
Tests of eddy diffusion in pipelines are analyzed and compared with Taylor’s 
theoretical solution. The pipeline coefficients are greater, as expected, be- 
cause of obstructions and curved sections in the pipe. 


INTRODUCTION 


Diffusion in ponds and pipelines is a combination of molecular and eddy 
liffusion.2,3,4,5 Molecular diffusion is the result of the heat motion of in- 
jividual molecules and is extremely slow in liquids and gases, because the 
mean free path of molecules is so short. Eddy or turbulent diffusion, while 
‘undamentally the same process as molecular diffusion, is the mass transfer 


ne  ———————————————————————————————————————_—————————————— 
Note.—Discussion open until October 1, 1961. To extend the closing date one month, 
. written request must be filed with the Executive Secretary, ASCE, This paper is part 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 87, No. HY 3, May, 1961. 
1 Engrg. Leader, Oak Ridge Natl. Lab., Union Carbide Corp., U.S. Atomic Energy 
yomm,, Oak Ridge, Tenn. 
“Modern Developments in Fluid Dynamics,” by S. Goldstein, Clarendon Press, Ox- 
ord, England, Chapter 5, 1938. 
3 “Turbulent Diffusion,” by G. K. Batchelor, Lecture Series, No. 4, Univ. of Mary- 
and, College Park, 1960. 
4 «tntroduction to Some Topics in Turbulence,” by F. N. Frenkiel, Lecture Series, 
lo. 3, Univ. of Maryland, College Park, 1960. 
5 “Statistical Theory of Turbulence,” by S, Goldstein, Lecture Series, No. 6, Univ. of 
Maryland, College Park, 1960. / 
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of molecules in eddies of various sizes. Essentially, molecular diffusion is 
submicroscopic in movement, while turbulent diffusion is macroscopic. Eddy 
diffusion has been found to vary from 1,000 to 1,000,000 times greater than 
molecular diffusion in open bodies of water.® Therefore, molecular diffusion 
will be ignored in the computations. 

Turbulent diffusion may be taking place in numerous ways within the same 
area. Large-scale eddies and small-scale eddies within them induce various 
rates of mixing, The complete randomness of the eddy diffusion process makes 
it a difficult phenomena to formulate, and it is only recently with the applica- 
tion of statistical methods to the problem that some ideas of the fundamental 
nature of the process are being uncovered. At present, (1961), it is only pos- 
sible to approximate the magnitude of the diffusion rate. 


HISTORY OF DIFFUSION STUDIES 


The earliest work on turbulent diffusion was that of Osborne Reynolds who 
differentiated in turbulent flow between the primary and secondary character- 
istics of the motion. The separation of turbulent velocity into two components 
is accomplished by assuming the mean velocity along the direction of travel 
as the primary component and a superimposed turbulent velocity having a 
temporal mean velocity of zero as the secondary component. This can be car- 
ried to an extreme by taking the time and length scales small enough so that 
even laminar flow would appear to be turbulent. 

Prandtl advanced turbulence theory by relating the intensity of turbulent 
flow to a parameter of distance. The “mixture length” of Prandtl represents 
the average distance a small mass of fluid will travel before it loses its 
increment of momentum to the region into which it comes.’ The work of 
Prandtl, however, bore no close relationship to the physical nature of turbu- 
lence. G. I. Taylor next applied the concept of the transfer of vorticity rather 
than momentum, but still retained the concept of mixture length which repre- 
sents the transverse path in the vorticity transfer. This formulation accorded 
better with experimental results. Theodor von Karman, Hon. M. ASCE, at- 
tempted to improve these efforts with his thesis of turbulent similitude based 
on the assumptions that the mechanism of turbulence is independent of visco- 
ity, except in the immediate neighborhood of the flow boundaries and that the 
pattern of secondary flow is statistically similar from point to point, varying 
only in time and length scales. These assumptions lead to equations of the 
same form as Prandtl’s and Taylor’s but withthe mixture length a function of 
the velocity distribution of the mean flow. A difference between momentum 
transfer and vorticity transfer is that in momentum transfer the eddy coeffi- 
cient is a function of a mixing length and the transverse velocity of the stream, 
while in vorticity transfer the eddy coefficient is independent of the secondary 
velocity. 

These methods lead to relatively simple and semi-empirical solutions 
which are widely used in engineering applications. For more adequate theo- 
retical treatment and closer representation of the physical reality, they are 


10a. 
6 «The Oceans,” by H.U. Brendan, M. W. Johnson, and R.H. Fleming, Prentice-Hall 
New York, 1942, p. 91. 


7 “Fluid Mechanics for Heavbalic Engineers,” by H. Rouse, McGraw-Hill Book Co, 
Inc., New York, 1938, p. 186. 
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not sufficiently accurate to enable fundamental research into turbulent flow 
te proceed. The fundamental problems of turbulence are best approached 
through statistical analysis. Taylor produced the first study of turbulence by 
statistical methods. The notions introduced by Taylor in the investigation of 
isotropic turbulence were the study of the correlation between two fluctuating 
quantities in the turbulent flow, the study of the decay of turbulence, and the 
study of the spectrum of turbulence (that is,the distribution of energy among 
eddies of different sizes). The foregoing devices are means to solve the main 
question of the joint probability distribution of the displacement of any “n” 
marked particles. 

The theoretical work in the field continues at an ever-eccelerating rate by 
such physicists and applied mathematicians as Heisenberg, Chandrasekhar, 
Onsager, Von Weizacker, von Karman, Batchelor, Townsend, Dryden, Gold- 
stein, de Feriet, Loitsiansky, Taylor, and Kolmogoroff, This research has 
been mainly concerned with laminar flow in the boundary layer, the initial 
stages of transition from laminar to turbulent flow, and the theory of non- 
viscous compressible and incompressible flow. As Dryden wrote® in 1952: 


“Progress of knowledge in fully developed turbulent flow processes 
has been at a more moderate rate, and our present methods (of solu- 
tion) are still semi-empirical in character.” 


The most frequent formulations used by engineers continue to be some form 
of the mixing length equations. 


MATHEMATICAL THEORY OF DIFFUSION 


Semi-empirical formulations of the mixing-length type will be used in this 
study. Recognizing that the diffusion rate is afunction of the size of the eddies 
90ssible, an average size of eddy andthe corresponding empirical is obtained. 
It should be noted that in the ocean or a reservoir there are many sizes of 
eddies operating within one another to increase the transfer area so that a 
sonstant concentration is attained. Because of the random nature of the 
addies, the concentration distribution in the ocean or reservoirs cannot 
follow an exact mathematical law, but must be treated in the statistical 
sense. That is, if enough tests are conducted the averages will tend to 
stabilize. Therefore, single tests must be approached warily. 

The diffusion equation used is of the form 


or the elementary case of no convection in one dimension, in which, c is the 
‘oncentration, a2 denotes the eddy diffusion coefficient, x is the length, and 
represents time. 

The coefficient a2 is the constant of proportionality between the rate of 
liffusion across a plane area and the concentration gradient normal to that 
rea; it is analogous to thermal conductivity in heat flow. Eq. 1 is equivalent 
o the general equation of conduction in heat flow. The vast amount of litera- 


Ng ee a ay |, eat eyaae crv on 
Fluid Mechanics Supports Tasks of Civil Engineering,’ by H. L. Dryden, Civil En- 
| hes Vol, 22, No. 5, 1952, p. 329, 
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ture on heat transfer is thus available for the more involved studies of! 
diffusion. The most important source of solutions of heat-conduction problems 
has been the text by H. S. Carslaw and J. C. Jeager.? 


PREVIOUS STUDIES OF EDDY DIFFUSION COEFFICIENTS IN OCEANS 


Among the early investigators, H. U. Sverdrup, M. W. Johnson, and R. H. 
Fleming have tabulated previous results for eddy diffusion coefficients, 1¢ 
Their values for horizontal diffusion are 2 x 103 to 4 x 105 ft sq per sec anc 
for vertical diffusion 10-3 to 1 ft sq per sec. This would imply that vertical 
diffusion can be ignored when compared to horizontal diffusion. They note that 
the eddy diffusion coefficient is a function of the size of the eddy. 

Later, Sverdrup, assuming diffusion to be a function of average velocity 
and of mixture length (that is, the dimension of the eddy) and independent o: 
the depth, hazarded a guess atthe values of diffusion coefficient.11 For eddies 
of radius of 3 x 103 ft, 3 x 104ft, 3 x 10° ft, and 3 x 108 ft, the diffusion coef- 


TABLE 1.—BURKE’S EDDY DIFFUSION COEFFICIENTS 


Diffusion coefficient, Diffusion coefficient, in ft sq per sec 


Eddy radius = 100 ft|Eddy radius = 200 ft 


(b) Taking into account the size of 
the eddies 


Test Number 


in ft sq per sec 


(a) Ignoring the size 
of the eddies? 


4 The diffusion coefficient equals a constant. 


ficients were estimated to be of the order of magnitude of 2.0, 5 x 102, 5x 103 
and 5 x 104 ft sq per sec, respectively. Sverdrup suspected that these value 
might be too large. 

After this, C. J. Burke, in a series of tests12 off the coast of Californi 
where he photographed dye-diffusion patterns from anairplane, arrived at th 
eddy diffusion values shown in Table 1. 

The values of the diffusion coefficients obtained in this test were used t 
estimate the spread of the contaminant for “Operation Crossroad” (the 194 
Atomic Energy Commission tests in the Pacific Ocean). 


————— 
9 “Conduction of Heat in Solids” by H.S. Carslow and J.C. Jaeger, Clarendon P 
Oxford, England, 1947, : ger, ndon Pres; 


10 “The Oceans,” by H. U. Sverdrup, M. W. Johnson, and R, H, Fleming, Prentice 
Hall, New York, 1942, p. 484, 4 
11 “Speculations on Horizontal Diffusion in the Ocean,” by H. U. Sverdrup, Scrip; 
Inst. of Oceanography, Oceanographic Reports No.1, January 7, 1946, p. 1 (unpublished 
“Horizontal Diffusion of Dyes inthe Ocean,” by C.J. Burke, Scripps Inst. of Oceanc 
graphy, Oceanographic Reports No, 2, September 1, 1946, p. 9 (unpublished), 
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W. H. Munk, G. C. Ewing, and R. R. Revelle measured the concentration of 
reef water at various parts of the Bikini lagoon and arrived at the following 
values:13 The vertical diffusion coefficient due to wind shear is 0.2691 ft sq 
per sec; the horizontal diffusion coefficient for an eddy size of 103 ft is 0.5 ft 
Sq per sec; and for an eddy sizeof 3 x 10° ft, itis 5 x 104 ft Sq per sec. They 
present as a primary constant the eddy diffusion coefficient divided by the 
radius of the eddy which they found equal to 0.0065 to 0.0164 ft per sec with 
0.013 ft per sec being the average. These values were obtained assuming one- 
dimensional flow in the analysis of the test results. 

Measuring the values of a2 for the dispersion of radioactivity following the 
underwater atomic blast, they arrived at a mean value of a2 for surface 
spreading of 1.6 x 102 ft Sq per sec for r equal to 9.83 x 103 ft. This yields a 
value of a2/r of 0.0164 ft per sec. The coefficient of diffusion at a depth of 
150 ft was about one-third the coefficient of diffusion at the surface. 

H. Stommel, by photographing sheets of paper scattered from an airplane 
onto the surface of the sea at Woods Hole, Mass., found!4 a series of values 
of a2 ranging from 2.252 x 10-3 to 8.256 x 102 ft sq per sec. A fairly straight 
line was obtained when the diffusion coefficient was plotted against the 4/3 
o0wer of distance separating the sheets of paper. The equation used is 

Rima ae oo ee ee (2) 
in which L denotes the distance between the sheets or particle, and € is a con- 
stant. This is the Richardson Lawof “neighbor diffusivity” and F(L) ~ 3.03 a2. 


FIELD TESTS FOR DIFFUSION COEFFICIENTS IN RESERVOIRS 


A series of tests were conducted to determine the diffusion rates in reser- 
foirs, Radioactive tracers, because of their lowcost, ease of detection, small 
sravimetric quantities required, and absence of density problems, were used 
0 measure the eddy diffusion coefficients.15 Various reaches of Reservoir 
No. 2 of the Metropolitan District Commission at Framingham, Mass. (Fig. 1) 
vere used to determine the eddy diffusion coefficient in quiescent bodies of 
vater. This reservoir is completely isolated from the potable water supply, 
und fishing and boating are prohibited. 

Method.—The tests were run in the following manner: (a) a small cove or 
asin or section of the reservoir was selected; (b) transit stations allowing 
ull view of the reservoir were established; (c) a dye and radioisotope mixture 
vere affixed to the end of a fishpole and deposited at a predesignated spot in 
he reservoir; (d) the dosing apparatus was removed, and a test tube was 
ttached to the end of the fishing pole; (e) the rowboat then circled the initial 
pot and water samples were taken in aseries of concentric circles about the 
nitial dosing point. Because of the action of the wind and water currents, the 
ate of dispersion of the tracer seldom exhibited radial symmetry about the 
dint of dosing. The dye provided a useful guide for the selection of sampling 


a 
13 “Diffusion in Bikini Lagoon,” by W. H. Munk, G, C. Ewing, and R. R. Revelle, 
'‘ransactions, Amer, Geophysical Union, Vol. 30, No. 1, 1949, p. 62. 
14 “Horizontal Diffusion Due to Oceanic Turbulence,” by H. Stommel, Journal of 
farine Research, Sears Foundation, Vol, 8, No, 3, December 30, 1949, pp. 205-214, 
15 “Radioactive Tracers in Hydrologic Studies,” by F. L. Parker, Transactions 
mer, Geophysical Union, Vol. 39, June 1958, pp. 434-439. ane a 
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points following the period of initial sampling. At the transit stations angula 
measurements were taken to locate each sampling point. By intersection c 
the angular measurements from two stations the exact points of samplin 
were plotted on the map. This system allowed greater flexibility than coul 
have been attained using a fixed set of predetermined sampling points. Wit 
the 12-ft fishing pole it was possible to approach the dye in the boat withot 
inducing too much turbulence in the area. The dye,sodium fluorescein, wa 
chosen only after a series of tests showed it to have the maximum visibilit 
in low concentration of all dyes tested (eosin, picric acid, Bordeau red, scarle 


INLET AT DAM 
RESERVOIR NO. 1 


DEAD INLET 


UPPER INLET POND 


ROUTE 135 SCALE IN MILES 


FIG, 1.—FRAMINGHAM RESERVOIR NO, 2 


4R conc., acid black conc., malachite green oxalate, ecio glavcine, fuch 
crystal). It is also the dye ‘chosen for use in air-sea rescue kits, A typig 
reservoir test is shown in Fig. 2 and Table 2. 

Analysis.—In the analysis of the test, the following mathematical mod 
was used: A “slug” of a pure radioisotope is added at a point in a lake ( 
infinite expanse and constant depth) and rapidly mixed from top to bottom | 
as to form an instantaneous line source of constant strength, The mixing 
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the water is assumed to be isotropic. Thereafter, the isotope is subject to 
two actions: (a) radial diffusion away from the source and (b) nuclear disin- 


TABLE 2,—SAMPLING POINTS OF AUGUST 14, 1951 TEST 


Time After 
~ sci is 


August 14, 1951 - Warm, sunny, no wind 


Comment 


al Extremity 
2 —-e ee ee 
3 Picked up bottom 
4 Picked up bottom 
5 Picked up wind shift from southerly to westerly 
6 Surface wind 
7 Surface wind 
8 Surface wind 
9 Surface wind 
10 Front of path 
11 Front of path 
12 Front of path 
13 Back of path 
14 Back bottom 
15 Back bottom 
16 Back bottom 
a7 Back bottom 
18 Front of path 
19 Densest path 
20 No dye 
21 Dense 
22 Densest 
23 Integrated depth sample 
24 Integrated depth sample 
25 Integrated depth sample 
26 Integrated depth sample 
27 Integrated depth sample 
28 Integrated depth sample 
29 Integrated depth sample 
30 Integrated depth sample 
HAN oe a ee) i a i oe 
feme ) 6Ste- ju’ be S9GAI Pe ee ane 
33 Grab samples in path 
34 Grab samples in path 
35 End of path 
36 In path 
37 In path 
38 In path 
39 Dosing point 
40 End of path bottom 
41 End of path bottom 
End of path bottom 


End of path bottom 


tegration ‘together with removal by self-purification processes (precipitation 
deposition, adsorption, biological uptake, etc.). 7 
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Consider (Fig. 3) an element of volumerdr dé dz in the plane of diffusion, 
in which r is the radius from the source, d@ the angular width of the element 
and dz the height. The rate at which an isotope enters the element by diffusion 
is 

= a” rdé dz pa 
or 


The rate at which it leaves by diffusion and die-away is 


2 8 8c 
- a" (r+dr) dé dz = eeerar) +krdrdédze 


in which a2 denotes the eddy diffusion coefficient reflecting the degree of 
mixing (length2/tim)E) and kis the effective decay constant of the radioisotope. 


ae 
ay 
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FIG, 2.—RESERVOIR TEST OF AUGUST 14, 1951 


FIG, 3,—ELEMENT 3 


Equating the difference between these two rates to the rate of accumulation 
of isotope within the element, one obtains 


1 k _ 1 ee 
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The appropriate solution of Eq. 3 for the boundary conditions stipulated is 


given by 
2 
N 
ea exp ( 5 xt) ee (4) 
4nat 4 t 
a 


in which c is the concentration at any radius and at any time after the isotope 
is added and N represents the amount of radioactive material initially added 
per foot of depth. An ideaofthe functional relation existing between c, r, andt 
is given by (Fig. 4). The concentration distribution takes the form of a sym- 


| 


FIG. 4.—DIFFUSION FROM AN INSTANTANEOUS LINE SOURCE 


metrical bell-shaped hill (Gaussian frequency distribution) with circular con- 
centration contours. With increasing time the hillspreads out and the concen- 
tration at the center decreases; also as time goes on the volume of the hill 
decreases because of radioactive disintegration and self-purification. 

It is instructive to consider the history of one particular concentration 
contour. Let r' be the radius of the circle having a concentration of c’ which 
Might correspond to a critical toxic level or tolerance concentration. From 
Eq. 4 it is found that 


N , 
ri=2Yyat ae Fa np aa iesdeae Cate deta 


a ne 
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Eq. 5 indicates that the critical radius at first increases rapidly, then passes 
through a maximum, and with progressive disintegration recedes to zero. The 
recession is rapid toward the end. At any time the concentrations within the 
critical contour are greater than c’,.and outside of it they are less than ck 
Of especial significance is the way the initial concentration, N, appears in 
the formula, A tenfold or hundredfold increase in N has a relatively small 
effect on r', 

Results.—The results of this type of analysis of the pond tests are shown 
in Table 3. These values of a2 /r are somewhat higher than those obtained by 
Munk, This may be explained, in part, by the circumstance that Framingham 
Reservoir No. 2 is relatively shallow, and uniform concentration distributions 
in the vertical direction were quickly established. In the ocean, however, 
steep concentration gradients persist for long periods following the addition 
of the tracer; the resulting continued downward transfer of the tracer causes a 
reduction in the dispersion inthe horizontal directions. From the data obtained 
previously, eddy diffusion coefficients in oceans and ponds versus the size of 
the eddies are plotted in Fig. 5. Ignoring Sverdrup’s estimates, it is possible 


TABLE 3.—RESULTS OF FRAMINGHAM RESERVOIR TESTS 


Date (1951) r, ft a”, sq ft per a2/r 
sec 
682 144 ot 


to plot a line of average eddy diffusion coefficients versus the radius of the 
eddy. The equation of the line is 


PL eT Ns 


The results of this investigation indicate that Richardson’s law of neighbor 
diffusivity obtains more than Munk’s suggestion of a constant for the eddy 
diffusion coefficient divided by the radius of the eddy. The use of a power 
function, rather thana linear relationship, would accord more with the physical 


system where with increasing radius the diffusing surface would increase 
faster than linearly. 


-4 (1.47 
© 


DIFFUSION IN PIPELINES 


Sir Geoffrey Taylor showed theoretically and experimentally that contam- | 
inants in pipelines diffused in the following manner:16,17 


16 “Dispersion of Soluble Matter in Solvent Flowing Slowly Through a Tube,” by Siz 
nc rps Seiya Proceedings, Royal Soc. of London, Sect. A, Vol. 219, No. 1137, August 

, , p. 186. 

17 “The Dispersion of Matter in Turbulent Flow Through a Pipe,” by G, I. Taylor, 
Proceedings, Royal Soc. of London, Vol. 223, No. 1155, May 20, 1954, p. 446. 
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FIG, 5.-RELATIONSHIP OF EDDY DIFFUSION COEFFICIENTS TO THE SIZE 
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Laminar Flow.— 


a. Dissolved matter spreads through the combined action of molecular 
diffusion in a radial direction and the variation of velocity over the cross- 
section. 

b. The distribution of concentration of dissolved matter is centered about 
a mean which moved with the mean speedof the flow and is symmetrical about 
it. 

c. The dispersion along the tube is in accord with a vertical coefficient of 


diffusion equal to 5 4) in which r is the radius of the pipe; u denotes the 
a 

mean velocity of the flow; and a2 represents the coefficient of molecular dif- 

fusion of the injected material in the fluid. 

Turbulent Flow.—Taylor predicates that the universal distribution of 
velocity in a pipe obtains (that is, the ratio of the difference between the veloc- 
ity at the center of the pipe and the velocity at any point in the pipe to the 
shear velocity, which is the square root ofthe shear stress at the wall divided 
by the density of the fluid, is a function of the radial distance to the point in 
question divided by the radius of the pipe). Further, it is assumed that Rey- 
nold’s analogy is valid (that is, the transfer of matter, heat, and momentum by 
turbulence are exactly analogous). 


a. Matter is transferred longitudinally by convection (the mean velocity) 
with a virtual coefficient of diffusion equal to 10.16 times the radius of the 
pipe times the shear velocity. The shear velocity is a function only of the 
Reynolds number for smooth pipes (that is, the coefficient of friction). 

b. Matter is transferred longitudinally by turbulent (secondary velocities) 
diffusion, with a coefficient of diffusion equal to 0,052 times the radius times 
the shear velocity. (Taylor assumes that the coefficient of longitudinal diffusion 
is equal to the coefficient of lateral diffusion.) 

c. The combined value for longitudinal diffusion coefficient is, therefore, 
equal to 10.1 times the radius times the shear camper In the American sys- 
tem of notation: 


Som 786 Din yiomek es. 


in which at is the longitudinal diffusion coefficient, D denotes the diameter, 
and f is the Weisbach-Darcy friction coefficient. 


For open channels 


a? = 14.28 oa V 2.8 By ene one aN 


in which R is the hydraulic radius, g represents the gravitational acceleration, - 
and, s is the slope of the channel, { 
Most of the previous theory has held that the transfer due to the superim- ~ 
posed turbulent component of velocity is the predominanttransfer. Therefore, — 
most theoretical formulas for eddy diffusion in pipes have been in error by a 
factor of about 200. This is because they have completely ignored the trans- — 
fer which is dependent on the primary velocity. 
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d. The concentration of dissolved material in the pipe at time t and distance 
x is equal to 


=M-=2. -3/2 » c aha? t 
OSS NOS eA Ae (9) 
in which M is the total mass (or activity) of the dissolved material. At any 
given time after dosing the tracer, the concentration distribution is symmetri- 
cal and gaussian. At a fixed station downstream the curve of concentration 
versus time is not symmetrical due to the diffusion of the material as it passes 
by the measuring station. Ordinarily, skewness of the concentration-time dis- 
tribution is not great when the measuring point is at a distance downstream 
greater than 100 diameters. 

e. Numerical values of the virtual coefficient of diffusion may readily be 
determined by taking advantage of the fact that the Eq. 9 is the error function 
equation. Taylor computed a2 from a plot of the data by scaling off one-half 
the distance at which the concentration exceeds 50% of the peak value. Solving 
the equation the following relation is obtained: 


in which Xo 5 is one-half the distance between sections at which concentration 
equals 50% of the peak value, and is the maximum concentration at the meas- 
uring station. 

Eq 10 may be transposed to the form 


-3 (t i? 
ae 2 _ altos (11) 


iby 4xl1n 2 CaP 0 8. 1o we. oer wine: ee 


in which to 5 is one-half the time that the concentration is above 50% of the 
maximum concentration, and x is the distance from the dosing point to the 
measuring station. 
f. Experimental work indicated that the formlation is validfor both smooth 
and rough pipes, since the universal distribution of velocity holds for both. 
g. For curved pipe the universal distribution of velocity does not hold and, 
hence, the coefficient of diffusion should be expected to be greater than that 


computed from Eq. 6. 


OTHER METHODS OF SOLUTION 


Other treatments of Eq. 9 for the evaluation of a2 are possible. Eq. 9 may 
‘be written in the form ; 
1/2 M (x-u t) 
log. wert =) OD 1084 9° Sides 
a Be a(n ra? 4 at t 


(12) 


In this form a linear relation obtains between log c t1/2 and (x-u t)2/t. 


“On a semi-log plot of c t1/2 versus (x-u t)2/t, the slope will be logyge/ (4 a?). 


GR err pe 
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Another method of linearization is obtained by writing Eq. 9 in the form 


2 -2 
et Pigs bitee) Ve tale oes @ 
log = |~“ —“_s- _—_ - __~~__“ | loge... (13) 
Jt 4a” 5 ey este 
“on “2 i L 
2 
c Dest kua gt 
1 
in which log 1Vi4 is the ordinate; cc is the abscissa; and 2h is the 
Coats 1 ‘2 4ay 


slope of the equation. 

It is also possible to evaluate the standard deviation from moments of the 
observed concentration-time distribution. The first two moments of the con- 
centration-time distribution, m, and mg, (Eq. 9) are 


ee 2 cet 
aioe aE Ye Soe NT So ae (14) 
x 
and 
paomccde es Aaah gt fa hm 
nig Pile Cha TS oat tet a ales a eae (15) 
> x 


in which T is the theoretical detention period, X/u, and X represents the total 
length of pipe. 

The second terms on the right-hand side of Eq. 4 and 5 are ordinaril 
small and may be neglected. If they are neglected the following equation for ar 
may be obtained by solving Eq. 14 and 15 for oe 

2 m 
2 Pe a sa AY me ay eager hi. (16) 


ee Gr 
(™) 
A final method is the approximate formulation used by R. S. Archibald, M. 


ASCE, in which18 the standard deviation of the concentration-time distribution, 
2 oe t, is assumed to be one-sixth the time between the first appearance and 


last trace of the tracer at the sampling station. The imperfection in Archibald’s 
criterion is the determination of the time of occurrence of the last trace, 
since the latter is obscured by fluctuation in background activity. Therefore, 
this method of determining the diffusion coefficient is much less accurate 
than either the method of moments or the Taylor method. 


DIFFUSION COEFFICIENTS COMPUTED FROM EARLY SALT-VELOCITY __ 
TESTS IN PIPELINES 


The results of pipe tests reported in the literature have been analyzed, 
and the diffusion coefficients have been determined by the methods described. 
The results are presented in Table 4. 


18 “Radioactive Tracers in Flow Tests,” by R. S, Archibald, Journal, Boston Soc, of 
Civ. Engrs., Vol. 37, No. 1, 1950, p. 77. 
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Linearization (Eq. 12) has the disadvantage that in forcing the data into a 
normal curve, the experimental points that do not fit into this pattern produce 
some scatter, and the determination of the straight line entails a strong ele- 
ment of subjectivity. Moreover, the physical significance of the data is lost in 
linearization. In the Taylor method (Eq. 11) of computing the diffusion coeffi- 
cient the physical sense is retained; it is much simpler and more straight 
forward. 

Taylor’s method requires a knowledge of the mean time of flow (x/u) and, 
therefore, is not completely satisfactory. None of these methods give satisfac- 
tory résults if instantaneous injection is not made or if an obstruction occurs 
in the line so as to reduce the effective diameter or if the pipeline is not 
straight. Taylor’s method (Eq. 10 and 11) does not indicate this. Therefore, 
the Taylor method should be carried out in conjunction with the method of 
moments. The method of moments entails no bias as the mean time of flow 
remains unknown. Then if the proper conditions have obtained, the results 
of the two methods should be approximately equal. 

A comparison of the values in Table 4 of the diffusion coefficients obtained 
by the various methods makes it evident that Taylor’s formula (Eq. 11) holds 
for both oil and water and for pipes flowing full and partly full. 

Ordinarily the coefficient of diffusion obtained experimentally may be 
expected to exceed the theoretical coefficients, since the Taylor formulation 
pertains to straight pipe. Existing pipes usually have curved sections. A few 
of the values in Col. 13 of Table 4 are less than unity. The discrepancy is 
probably due to errors in the determination of the friction coefficient of the 
pipe, Also, as mentioned previously, the linearization method is dependent on 
the choice of the line of best fit. 


FIELD TESTS OF PIPELINES 


A series of tests using radio-tracers was made to determine the virtual 
coefficient of longitudinal diffusion in pipelines. These results were then com- 
pared with similar coefficients obtained from previous pipe tests and with the 
coefficients computed by the Taylor formula, 

The tests were' run in an old water-supply pipeline made evailehte by the 
Metropolitan District Commission and the Conservation Department of Mass- 
achusetts. The pipe tested was part of the conduit system between Whitehall 
Reservoir and Hopkinton Reservoir, These are out-of-service water-supply 
reservoirs and have been assigned to the Department of Consevation. The 
pipe is 20 in. in diameter and 6,400 ft long. The drop in the hydraulic grade 
between the inlet and outlet of the line is 34 ft. The method of testing was to 
allow the pipe to flow full long enough to establish steady state conditions and 
then to “slug” the pipe with a dose of radioactive tracer and dye. At the outlet 
of the pipe a sampler was arranged so that the outlet was sampled every 1 min. 
until appearance of the dye and every 1/2 min. thereafter until the dye dis- 
appeared. The sampling was continued for five or six detention periods. The 
concentration of the dye served as a rough check for radioactive measure- 
ments. A 50-ml volume of each sample collected was precipitated.The pre- 
cipitate was dried and counted with an end window Geiger tube with suitable 


sealing circuits, The results were compared with a diluted sample of the 
initial dose. 
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The results of two of the tests are shown in Figs 6 and 7. The method of 
moving averages was used to compute the tails of the curves. 

In Table 5, a much greater percentage of iodine passes the downstream 
gauging station than does phosphorus. Only a small amount of phosphorus 
reaches the sampling station in test 2 due to the high absorption rates of 
phosphorus by the bacterial slimes on the pipe walls.19 Further credence to 
this theory is providedby C. C. Coffinet al., who state20 “p32 uptake by plants 
and sooplankton is a matter of minutes or hours and not days or weeks.” It will 
be noted in Table 5 that the coefficient of diffusion by Taylor’s method (Eq. 11) 
and the method of moments (Eq. 16) do not agree too well. 

The apparent reason for the failure of the method of moments to give re- 
sults that agree with other methods is: 


1. The method of moments (Eq. 16) uses “tails” of concentration-time data 
which are obscured by variations in background counting rate. 

2. Taylor’s theory is not good for those portions of the distribution-time 
curve at large distances from mean. This is evident from the shape of the 
curves. 

3. The Taylor method (Eq. 11) is not good for skewed curves. The curves 
for the pipe test were skewed because of the air pockets which remained in 

the pipe. In the last test, however, the pipe was vented before the test began, 
and we can assume that most ofthe air was removed. The method of moments 
then gives approximately the same results as other methods. 

Another series of tests were run through a sewer line at New Bedford, 
Mass., in co-operation with Frank Heaney, F. ASCE, of Fay, Spofford and 
Thorndike, consulting engineers. The sewer tested was the main collecting 
sewer of New Bedford, Mass. Its route is along Rodney French Boulevard and 
through Fort Rodman, and it discharges into Buzzards Bay (Fig. 8). The sewer 
is of the usual horseshoe cross section of combined sewers with a maximum 
width of 92 in. and a height of 84 in. The slope is approximately 0.001. In one 
test the activity was introduced at Hudson Street. Samples were collected with 
diaphragm pumps at downstream stations at the drill ground (1,987.7 ft) and 
at the tide chamber (3,429 ft). At these pump outlets the stream was sampled 
every 1 min until the dye appeared and thereafter at 20-sec intervals until the 
dye disappeared at which time 1-min sampling was resumed; typical curves 
are shown in Fig. 9. 

In another test the tracer was dosedatthe screen house, 1,060 ft upstream 
from Hudson Street. The same sampling procedure was used with an additional 
sampling station at Hudson Street. Due to faulty injection the activity in the 
test bled into the system from the screenhouse well slowly. This is reflected 
in the high values for a2, In these conditions it should be expected and it is 
found that the apparent coefficient of diffusion is much larger than the theo- 
retical value. This emphasizes the necessity for instantaneous injection into 
the turbulent flow. Injection into the laminar sublayer will obviously produce 
laminar diffusion and later turbulent diffusion as it bleeds into the turbulent 
flow. Attempts to measure the eddy diffusion coefficient at the outfall were 


EEE LL 
19 “Adsorption and Assimilation of P32 by Bacterial Slimes,” by G. W. Reid, San, 

Engrg.» Dept., John Hopkins Univ., Contract AT-(10-1)-498, March 20, 1950, p. 7. 

: “Exchange of Material in a Lake as Studies by the Addition of Radioactive Phos- 

phorus,” by C.C. Coffin, F.R, Hayes, L.H. Jodrey, and S, G, Whiteway, Canadian Jour- 

nal of Research D, Vol. 27, 1949, p. 221. coe b he 
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10t successful due to the small amount of radioactive isotopes used and the 
nigh solids in sea water. 

The tests in Massachusetts bear out the conclusions drawn from previous 
‘ests. With proper injection methods and measuring procedures the Taylor 
‘ormulation and the method of moments give satisfactory results for circular 
and horseshoe pipes, full and partially flowing full pipes, and oil, water, and 
sewage carrying pipes. 


CONCLUSIONS 


It is relatively easy to determine diffusion coefficients in reservoirs and 
ipelines by use of radioactive tracers in comparison with most other methods 
aeretofore used. Simple dosing and sampling devices can be used to obtain 
rood results. 
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EARLY HISTORY OF HYDROMETRY IN THE UNITED STATES2 
Closure by Steponas Kolupaila 


STEPONAS KOLUPAILA.®5_The writer feels particularly honored by the 
fact that three of the pioneers presented in the paper took part in the discus- 
sion: J. C. Stevens, F. ASCE, W. G. Hoyt, F. ASCE, and M. Blanchard, F. ASCE. 
It is not enough to thank all contributors for their valuable cooperation; let it 


 § 
i, 
F 


uffice to say that the writer is inspired, by the attention of his colleagues, to 
continue further study of the history of hydrometry. ~ 

é J. W. Johnson deserves credit for recalling the “Early Engineering Center 
in California,” a good addition to the chapter “Water in the West.” Hamilton 
nith is world renowned for his authoritative book86 that is now very rare. 


FIG. 63.—HAMILTON SMITH FIG. 64.—LOAMMI BALDWIN 


e classical experiments carried out by Smith are mentioned in modern text- 
ooks. The portrait of this pioneer in hydrometry (Fig. 63) is reproduced 


-@ January 1960, by Steponas Kolupaila (Proc. Paper 2335). Pn 

85 prof, of Civ. Engrg., Univ. of Notre Dame, Notre Dame, Ind. . 
_ 86 «Hydraulics, the Flow of Water Through Orifices, Over Weirs, and Through Open 
Jonduits and Pipes,” London, 1886, 
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here87 by the courtesy of K. Daugelaand N. H. Bedford. Smith died while sail- 
ing down Oyster River. 

Johnson mentioned in his discussionthat the first of the world’s long-distance 
high-voltage power transmission lines was built in California. In connection 
with this statement, G. H. Ellis, M. ASCE, insisted in a letter to the editor that 
the Missouri River Power Company had built a 70,000-v line extending from 
their Canyon Ferry plant to the town of Buthe, a distance of 63 miles, in the 
year 1901, thus claiming to predate the California project. The writer would 
like to point out, however, that the 60,000-v line running 140 miles from the 
Colgate plant to the city of Oakland, California, that was described by Johnson, 
was built in the year 1899, thus refuting Ellis’ claim. 

Hoyt, one of the leading pioneers in the USGS, from 1907 to 1950, recom- 
mended the inclusion of several hydrographers of that institution. It is the 
feeling of the writer that these persons could be more appropriately included 
in the modern history of hydrometry. Hoyt emphasized the manifold indirect 
duties of hydrometric engineers, as disseminators of public information and 
educators to the necessity of conserving our water resources. He expressed 
his opinion, as aconsultant forthe Senate Committee on Water Resources, con- 
cerning the future of those problems intimately affiliated with the history of 
hydrometry in the United States. The annual expenditure of the USGS for in- 
vestigationof water resources, about 25 million dollars, demonstrates the high 
recognition of this important work. 

E. Shaw Cole, F. ASCE, who skillfully continues the illustrious work of his 
father, E. C. Cole, mentioned some new developments, that, the writer feels, 
belong more appropriately to the modern history of hydrometry. 

Blanchard, one of the distinguished pioneers of the U. S. Lake Survey, added 
more informationon the measurements at Sault Ste Marie, and on hisown work 
in different rivers connecting the Great Lakes. 

A. H. Frazier, who inspired this study, contributed valuable improvements 
with which the writer agrees completely. The pioneer work of Loammi Bald- 
win, (Fig. 64) who has been referred to as the father of civil engineering in 
America, should be included in the early history of hydrometry. 

The writer is grateful to Frazier for disentangling the amazing coincidence 
of the two Nicollets, both investigators in one part of our country, but 200 yr 
apart. 

The first was Jean Nicolet, Canadian explorer, born in Cherbourg, France, 
about 1598, who immigrated to Canada in 1618, travelled across the Indian ter- 
ritory, at that time called New France, and was the first white man to visit 
Wisconsin, He didnot reach the Mississippi River, however, and later returned 
to Quebec. He drowned in the year 1642. A river, a country, and a town in 
Canada bear his name, 

The other was Joseph-Nicolas Nicollett, 1786 to 1843, French mathematician 
and cartographer, who came tothe United States in 1832 andexplored the upper 
Mississippi in the years 1836, 1838, and 1839. His report, written in French, 
was translated to English by J. T. Ducatel and J. H. Alexander, and published 
in Congressional records.8 ; 

The improvements to the material concerning W. H. Hall, State Engineer of 
California, and the portrait furnished by Frazier, are highly appreciated. While 


87 «History of the Town of Durham, New Hampshire,” 1913, 
Report intended to illustrate map of hydrographical basin of upper Mississippi 
River made by J. N, Nicollet. Senate Document 237, 26th Cong., 2nd Session, 1843, 
Serial No. 3PO; House Document 52, Vol. 2, 2nd Session, 1845, Serial No, 464, 
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it is true that several methods were initiated by Hall, who published a large 
volume of hydrometric data collected in California, 89 and several annual re- 
ports on irrigation from 1880 to 1888, his methods were mentioned in papers 
by his assistant C. E. Grunsky only, and therefore were attributed to Grunsky, 
although Grunsky himself did not claim them personally. This information 
concerning Hall and Grunsky was corrected in the recently published bibliog- 
raphy. 90 

J. C. Stevens, the eminent pioneer in hydrometry, was particularly favor- 
able to the writer’s intentions. He recalled the origin of the widely publicized 
theory of hydrologicimportance of our forests, that was deliberately falsified 
in order to further conserve the forests. This was done in connection with a 
comparison of droughts in the Ohio basin of 1824 and from 1925 to 1930, It 
was particularly interesting to know the personal opinion of one of the authors 
of that story. Careful hydrometric investigations carried out in experimental 
forest basins in different countries had thoroughly explained the true hydro- 
logic role of forests and clearings. 

I. J. Moskvitinoff presented some conclusions from his hydrometric work 
in Central Asia from 1914 to 1915. The extensive investigations of the accu- 
racy of flow measurements, the best investigations known to the writer, were 
included in the 5th volume of the Report of the Hydrometric Department for 
1914,91 that, regrettably, was never published. Fragments of this work ap- 
peared in Bulletins of the Hydrometric Department in Turkestan in 1917.92 

Among the many letters received by the author in connection with his pa- 
per was a humorous contribution received from George Henry Ellis, M. ASCE, 
who wrote: 


“Years ago, when I was swinging the current meter, there was oftena 
group of spectators asking foolishquestions. The foolishanswer I some- 
times gave them grew, over a period of about thirty years, into the story 
of the counting of the fish.” 


The enclosed ballad entitled “The counting of the fish” of 1946 contains verses 
like “fishes cannot live in water near a power plant, with kilowatts squeezed 
out.” This excerpt is a good example of hydrometric humor. 


ADDENDUM TO APPENDIX II.-SHORT BIOGRAPHICAL DATA 


Loammi Baldwin, b. May 16, 1780, at North Woburn, Mass., d. June 30, 1838, 
at Charlestown, Mass.; 1800 graduated Harvard Univ., built dams and harbors, 
1821 Union canal in Penn., 1825 in Boston, built Connecticut-Hoosac tunnel; 
called Father of civil engineering in America. 


_ 89 «physical data and statistics of California,” by W.H. Hall, Tables and memoranda 
of data, collected and compiled by the State Engrg. Dept. of Calif., Sacramento, Calif., 
1886. 

_ 90 “Bibliography of Hydrometry,” by S. Kolupaila, Notre Dame, Univ. of Notre Dame 
Press, 1961, See references D17, D157, andH 12. 

_ 91 “Bibliography of Hydrometry,” by S, Kolupaila, Notre Dame Univ., Notre Dame 


Ind,, Press, 1961. See reference Y 67, : 
mAyp “Bibliography of Hydrometry, by S. Kolupaila, Notre Dame Univ. Press, Notre 


Dame, Inc., 1961, See reference H 294, 


i 
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FIG, 65.—JOSEPH SAXTON ' 


FIG, 67.—G. W. RAFTER FIG, 68.—D. C. HUMPHREYS 
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William Hammond Hall, b. Feb, 12, 1846, in Hagerstown, Md., d. Oct. 16, 
1934, (only this information was available). 1879 State Engr. of California. 


David Carlisle Humphreys, b. Oct. 14, 1855, in Chatham Hill, Va., d. Jan. 
10, 1921, in Lexington, Va.; 1878 graduated Washington & Lee Univ., 1879 Mo. 
River investigations, 1885 to 1921 Prof., 1904 Dean, Washington & Lee Univ.; 
1895 to 1908 Geological Survey; 1897, James River expedition. 


Elwood Mead, b. Jan. 16, 1858, near Patriot, Ind., d. Jan. 26, 1936, in Wash- 
ington, D.C.; 1882 graduated Purdue Univ., 1886 Prof. of Irrig. Engrg. Colo- 
rado Agric. College; 1899 Chf., Irrig. Investigations, Dept. Agric.; 1907 to 
1913 river studies in Victoria, Australia, 1913 Prof., Univ. California, 1917 
Calif. State Bd. on Irrig., 1923, 1927, Cons., Palestine, 1924 to 1936 Commis- 


FIG. 69.—ELWOOD MEAD FIG, 70.—F. W. HANNA 


+ 


sioner, Bur. of Reclamation, promoter of Hoover Dam; reservoir on Colorado 
R. was named Lake Mead in his honor. 


_ Edwin S.(*) Nettleton, b. Oct. 22, 1831, at Medina, Ohio, d. April 22, 1901,. 
in Denver, Colo.; ed Medina Academy and Oberlin College, Ohio, 1870 Colo- 
rado town planner, 1883 to 1886 State Engr. of Colorado. (* ) Adopted initial, 


no middle name. = 


f Joseph-Nicolas Nicollet, b. July 24, 1786, in Cluses, Haute-Savoy, France, 


d. Sept. 11, 1843, in Washington, D. C.; cartographer, 1817 Paris astronom- 
ical observatory librarian, Prof. Louis College; 1830 ruined on stock market, 
. immigrated to U. S., 1836 to 1839 Miss. River investigations. 
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FIG, 72,—HYDROMETRIC TRAINING CAMP AT EMBUDO, N. M., 1888 
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Hamilton Smith, b. July 5, 1840, near Louisville, Ky., d. July 4, 1900, near 
Durham, N. H.; 1854 work in mining and hydraulic engrg, 1869 to 1881 in Calif- 
ornia, water supply for mining, 1881 Venezuela, London, 1892, South Africa, 
Alaska, mining, 1895 New York, cons. engr. 


Extending cordial thanks to all contributors, the author has taken the oppor- 
tunity to add several portraits located through the kind cooperationof Frazier. 
They are: J. Sexton, the first manufacturer of a current meter in the United 
States; E. S. Nettleton, State engineer of Colorado; G. W. Rafter, State engi- 
neer of New York; D. C. Humphreys, investigator of the James River (cour- 
tesy of the Washingtonand Lee University Library); E. Mead and F. W. Hanna, 
pioneers of the U. S. Bureau of Reclamation (courtesy of the Seattle Public 
Library). The interesting sketchof the early weir measurement93 is also pre- 
sented. (courtesy of Purdue University, Gorr Library). The historical group 
of the members at the first hydrometric training camp, near Embudo, New 
Mexico, in 1888, must also be included in the early history of hydrometry. 


93 «Construction of Mill Dams,” J, Leffel and Co., 1881, 
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GENERALIZED DISTRIBUTION NETWORK HEAD LOSS CHARACTERISTICS# 


Closure by M. B. McPherson 


M. B. Mc PHERSON, 2° M. ASCE.—The considered opinions, questions, ar- 
guments, clarifications, and original supplementary information presented by 
the thirteen engineers who contributed the nine discussions have placed the 
subject in much clearer perspective. 

Little would be gained by reviewing the discussions point by point or by cit- 
ing the multitude of statements with which the writer is in full agreement. 
Rather, the questions raised will be answered and clarifications offered. 

Specific Details.—The writer stands corrected by Robertsonon the improper 
use of the term “relaxation.” 

Constant mentions that the flows in each pipe can be checked during a net- 
work analysis by means of Eq. 6. For all of the network runs used in the pa- 
per, each junction and loop was checked to insure that the algebraic sums of 
all the corresponding flows and head losses were zero, with due allowance for 
attainable precision of meter readings. For the initial run of each network 
series, each and every pipe was checked via the equivalent of Eq. 6. There- 
after, on every second or third run all of the pipes were completely checked 
again to insure that all “fluistors” were operating properly. This is a stand- 
ing operating procedure of the Philadelphia Water Department; it leaves little 
chance for error. 

The error analysis by Constant is a valuable supplement. Some indication 
of relative error should have been included in the paper. 

Constant contends that the procedures presented are principally applicable 
0 field operations. Radziul and Celenza have clearly demonstrated the value 
of these procedures for operational purposes. However, an equally if not more 
important use is in the system analysis phase of design, as will be explained 
subsequently. 

Bitoun has referred to the number of conditions required to define Fig. 4 
ind for determination of ¢ and n in Eq. 2. The “existing” district in Fig. 3 was 
studied to secure a more exact understanding of its characteristics because 
xtensive, synchronized field tests were not practicable due to its large size. 
This preliminary study was of great value in subsequent studies of modifica- 
ions required forfuture conditions. The field data were included in Fig. 4 for 
erification. More than one run was made for the Fig. 4 presentation to prove 
mequivocally that Eq. 1 was valid under the prescribed conditions. It was felt 
hat a single run would not be adequate for that purpose with such a complex 
etwork, even though only one run is required to define Kj. It must be em- 
hasized that a major consideration in the designof public utilities is the pro- 
ision of adequate facilities for future needs. Seldom will an existing system 


_ 4 January 1960, by M. B. McPherson (Proc. Paper 2339). 
25 Prof, of Hydr, Engrg., Civ. Engrg. Dept., Univ. of Illinois, dir Ill. 
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be capable of handling anticipated higher future demands. Also, modern com- 
putational devices make possible more exact and realistic studies. Indeed, 
the inadequacy of many existing facilities can be attributed to the comparatively 
primitive computational procedures available and in vogue during the era in 
which they were designed. 

Radziul and Celenza have clearly demonstrated the efficacy of using Eqs. 1 
and 2 to determine the characteristics of existing prototype systems. The 
modifications planned for these two districts |[ Figs. 5 and 6(a)], either in pro- 
gress or scheduled, will vastly change their existing (Tables 5 and 6) charac- 
teristics. Because variations in piping and source and storage facilities can- 
not be simulated under prescribed loading conditions with a prototype system, 
they must be simulated by means of some type of model. For this simulation 
one network analysis is required to define Kj in Eq. 1 and two analyses are 
required to define ¢ and n in Eq. 2. Obviously, it matters not what computing 
device is used in analyzing a particular loading run, to paraphrase one of Lo- 
max’s observations. 

Lomax states that it was “not so stated, but it is implied that the source 
magnitudes must be proportional as well as the load magnitudes.” This would 
be true only in the specific instance of a network with a single source and na 
equalizing storage: Eq. 1 would apply for all system demands; each local de- 
mand load would be a fixed fraction of Qg; here, because Qp equals Qgq, the 
several loads would also be an equal fixed fractionof Qp. A network fed by twe 
sources is hydraulically quite similar toa single-storage, single-source com- 
bination (see Figs. 5 and 6) when the equalizing storage is contributing to the 
demand (Qp less than Qq). Two or more true sources, whether they be pumps 
or gravity storage, can be generalized by use of Eq. 2. The writer has made 
effective use of this procedure inthe system analysis for asingle-district large 
city that has two gravity storage reservoirs at dissimilar levels, a pumping 
station at a remote corner with a second station ina different quadrant needec 
for the future. 

The system analysis of the single-source case without equalizing storage is 
obviously a routine matter. With multiple sources, or with single or multiple 
sources combined with equalizing storage, the number of hydraulic parameter: 
that must be satisfied increase. The capacity per foot andrelative water level: 
of storage, together with pump characteristic curves, must be reconciled wit 
the network head loss characteristics under a variety of known or anticipate 
service conditions. The demonstration of storage-network-pump balancin; 
that will be presented is a relatively simple example of a final or equilibriur 
phase of a system analysis. The paper was addressed tothe issue of a gener 
alization of network characteristics. In preparing the paper it was feared tha 
a detailed treatment of system analysis would detract from and overextend th 
subject, but it is evident from some of the discussions that the paper shoul 
have been extended somewhat for a more firm appreciation of the potentiz 
uses of the concepts advanced. : 

Applicability of Proportional Load Principle to Non-Hydraulic Networks.: 
Trunk and Maier and Miller have clearly indicated that the proportional loa 
principle is of little value in the analysis of gas system networks. Trunk re 
ports (correspondence dated March 23, 1960) that prior to acquisition of 
McIlroy Analyzer in 1953, his staff “resorted to the trick of thinking of a net 
work asa single equivalent main,” when a complete analysis was not warrante 
via manual computation. The relation between the equivalent pipe method ar 
Eq. 1 will be discussed in a later section. | 
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Mine ventilation networks have been investigated by means of a special Mc- 
lroy Network Analyzer by the U. S. Bureau of Mines, Health and Safety Re- 
search and Testing Center, Pittsburgh, Pennsylvania. 6 The resistance in 
mine airways is proportional to the Square of the volume rate of flow. Al- 
hough the basic procedure for analysis of mine networks is in many respects 
similar to that for other fluid media, loading conditions are rather unique. E. 
Harris, Acting Chief, Ventilation Group, U. S. Bureau of Mines, Pittsburgh, 
— describes the situation thusly (correspondence dated May 19, 


“We do use a form of Eq. 1; however, the application of ‘proportional 
loading’ to mine ventilation would be very limited. Where minor fan ad- 
justments (speed or blade settings) are made with no other changes in 
the system a form of Eq. 1 and proportional loading would be and is ap- 
plied. Where adjustments in distributionare entailed, use of Eq. 1 would 
be limited to major airways in which the resistance factor remains con- 
stant. Other factors would enter into the application.” 


.§ in all network analyses, accurate results can be achieved only through use 
f accurate pressure and quantity field survey data. 

Alan C. Byers, Head, Electronics Branch, The Franklin Institute, Philadel- 
hia, Pennsylvania (correspondence dated May 26, 1960) has carefully studied 
1€ paper: 


“My conclusion, after such a study, was that the paper does not present 
an analysis useful for load flow studies in power system practice. In 
order to avoid having this negative opinion be unchecked, engineers of the 
Philadelphia Electric Company were asked to study the paper also. Their 
electrical engineers in system planning find that there is nothing useful 
to them in the paper. They also had engineers from their steam and gas 
distribution systems, which might be thought to present closer analogies, 
read the paper. Here the paper was found of general interest, but not 
useful in any specific way. It seemsthat in eachof these fields engineers 
have developed analytical methods welladapted to their problems. Doubt- 
less, by drawing the correct analogies, the methods developed in the pa- 
per might also be used in these other fields. But there seems to be no 
advantage to be gained.” 


Reversing the orientation of these considered comments, it is apparent that 
@ engineer engaged in water distribution network analyses should exercise 
ution intransferring or adapting to his use procedures or analogies that have 
en found satisfactory in the appraisalof electric power systems, gas, steam, 
‘mine ventilation networks. The word “water” deliberately was not incor- 
rated in the title of the paper because it was mistakenly assumed that Eq. 2 
ight prove useful in the analysis of non-hydraulic systems. 

Clennon2” has described a linear resistance analogue, and Mcllroy28 the 
n-linear type, for gas network analysis. 

Fire Flow Analysis via Proportional Load Hydraulics.—The statement: 
lowever, a fire load cannot be regarded as a representative proportional 


26 “A Network Analyzer for Solving Mine-Ventilation-Distribution Problems,” by G. 

McElroy, U.S. Bur. of Mines Information Circular 7704, November, 1954. 

27 “Calculators Solve Flow Problems,” by J. P. Clennon, Amer, Gas Assoc., Octo- 
, 1951, pp. 68-70. 

“0 Direct Reading Electric Analogue Computer Helps Solve Distribution Problems,” 

M, 8. McIlroy, Gas Age, July 17, 1952, pp. 31-35, ; 
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load,” was apparently misconstrued by the writers of two of the discussions. 
It was not intended to imply that fire flows cannot be analyzed via the propor- 
tional loading principle. Linaweaver, Geyer, and Wolff have correctly explained 
that a fire flow load per se is non-proportional, just as is the load imposed by 
elevated storage while being filled. An elementary network is shown in Fig. 
10 to illustrate this point. The fire flow is termed Qf and the head loss from 
Qp to the fire flow point is termed Zhy. The data for the four runs plotted in 
Fig. 10 are given in Table 8. It is apparent that Eq. 2 is appropriate for the 
analysis, so long as each non-fire load is a fixed proportion of the total non- 
fire load demand. 


TABLE 8.—SUMMARY OF SIMULATED NETWORK FIRE FLOW RUNS (SEE FIG, 10) 


McIlroy Analyzer Runs 


(6 = 0.386 and n = 2.85) 


Using the standard fire flow extrapolation formula to anticipate run F-4 fron 
runs F-1 and F-2: 


0.54 
Q = 4.00 (Gar) = 7.0 mgd 
and run F-4 from runs F-1 and F-3: 
0.54 
Q = 6.00 a) = 7.6 mgd 


versus the true F-4 value of 8.0 mgd. 

The elementary network given by Tong29 will be usedas a second example 
In Table 9 are given the appropriate details for solution via Eq. 2. Note the 
the true higher fire flow is 1,000 gom whereas for the standard formula iti 
740 gpm. 

These examples tend to exaggerate, to some extent, discrepancies betwee 
the standard formula and true values. Only the most remote point in the net 
work was considered. It may be noted also that the larger the drop to the ret 
erence flow (39.7 ft, 71.9 ft and 2.3 ft) the better the formula approximatio: 
It is understood that the National Board of Fire Underwriters (NBFU) no! 
mally requires a 10 psi drop before a fire flow sample is considered adequa’ 


pa a a a 
29 Discussion of “Friction Losses in Lines with Servi “Ton 

ce Connections,” by A. L. Ton: 
Proceedings, ASCE, Vol. 86, No. 9, November, 1960, p. 160. oe “4 
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4 rating purposes. Anderson reminds us that for adequate design, fire flows 
just be available under maximum day demand conditions. Whereas the stand- 
-d formula yields extrapolated flows that are conservative (that is, low), the 
BFU seldom extends field test data over a large increment beyond th 
ows to estimate capacities at a 20 psi residual. 
It should be noted that the standard formula is paradoxical, inasmuch as it 
strictly applicable only when the fire flow approximates a proportional sys- 


e test 


ce) +0.1 +0.2 +0,3 
LOG Q,/Q, 


FIG, 10,—SIMULATED NETWORK FIRE FLOW RUNS (SEE TABLE 8), 


load. This approximation is fairly exact in the extreme: when the fire flow 
n insignificant part of the total system demandor when it is practically the 
ire system demand. The first instance is a local phenomenon because the 
ct of the additional system head loss attributable to the transport of the fire 
7 through the network is relatively small. In the second instance the fire 


~ 
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flow occasions the entire system loss and the standard formula is then equiv- 
alent to Eq. 1. 

In the preceding two examples there is no equalizing storage. With equal. 
izing storage there would be flow out of the reservoir(s) and the fire flov 
might be regarded as an additional, filling reservoir without a fixed or con- 
trolled water level. The use of Eq. 2 for fire flows with or without equalizins 
storage in conjunction with machine network analyses would be pointless be- 
cause the design fire flow alone is required and normally can be determine‘ 
directly. However, if it is evident thata field test at agiven site would requir 
extrapolation to the required flow, a standard formula extrapolated flow base 
on an intermediate fire demand would be more realistic in terms of ultimati 
NBFU rating. 

Sweitzer30 has presented a clear description of the procedure for a grow 
hydrant flow test. 

Unique Solution Feature of Networks .—Bitounhas referred to a demonstra 
tion of uniqueness given by @’Auriac.12 The demonstration is an inferentia 


TABLE 9.—FIRE FLOW DATA FOR SMALL NETWORK GIVEN BY A. L, TONG 


@ Computed by writer. P Cemputed by Tong. n = 2.36, and = 9.25 x 10-6 


300 CH =aal ) = 740 gpm, versus Qp = 1000 gpm 


4.0 -1.7 


proof, limited to a topological consideration of boundary conditions too general 
treated. It is summarized: “The unicity of the solution depends uniquely « 
the fact that conduit characteristics are increasing functions.” This is to s: 
that in addition to being a continuous function the head loss characteristic 
pipes is bilateral, inasmuch as the magnitude of the loss is independent of t 
direction of flow. (The continuous, bilateral nature of resistive electrical in 
pedances31 makes possible the useof A.C. electrical analogues inthe analys 
of fluid flow networks.) A rigorous proof for m = 2 appears possible throu 
use of special formulations of matrix algebra. 32 Attainment of a rigorous prc 
of uniqueness for network equations with a non-integer valueof m is unlikel 
an inductive proof of uniqueness follows. 

Consider any balanced distribution network for which the flow in each mé¢ 
branch (both ends terminating at intersection junctions) may be described } 


he ee QMO LY Te ee 
Because the network in questionis balanced, the boundary conditions have be 


30 “Basic Water Works Manual,” by R. J. Sweitzer, Amer. Concrete Pipe Agel 
Chicago, 1958. “_ 
e 31 ipeoe in Electrical Engineering,” by C. R. Vail, Prentice-Hall, Inc., Nev 

ork, 1950. 

32 «The Arithmetic Theory of Quadratic Forms,” by B. W. Jones, Math, Assoc 

Amer., 1950, | 
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satisfied exactly: the algebraic sum of flows at each main branch junction is 
Zero, 


i=f 
LI “Sloat oon ewe ti lic a (14) 
i=1 ie ee 


mn which f is the number of Separate flows at a junction ( necessarily three or 
nore in number) and J is the total number of junctions; the algebraic sum of 
ead losses around each loop is zero, 


i=b 
Fe ee aa eee ee OTe (15) 
i=1 beng ROSES 


1 Which b is the number of main branches in a loop and L is the total number 
f loops. 

The combined number of loops L and junctions J in anetwork is equal to the 
tal number of main branches Xb plus one, Immediately the loads are as- 
igned values the branch flows remain the only unknowns. Because the input 
quals the sum of the loads the input junction equation is superfluous. The 
umber of equations L+J=1 equals the number of unknowns Xb and, therefore, 
‘Solution is always possible. 

For the balanced network under consideration Eqs. 13, 14, and 15 have been 
xactly satisfied. Rewriting Eq. 13 


h=(kQ™M-1) (Qi =k 1Ql 20... eee eee ee (16) 


sduces its original bilateral non-linear form to a bilateral linear equivalent. 
he sign for each h in Eq. 15 must be determined from the sign of the corre- 
ending balancing branch Q in Eq. 14 (note that when m = 2 the sign of h is 
trinsically positive in Eq. 13). Because the magnitude and sign of 1Q! is 
<ed and k is a constant, each ky has a singular constant value without sign, 
(ual to the magnitude only of (k QM-1), It follows that the transformation 
to linear equivalents can have ho bearing on the integrity of Eqs. 14 and 15. 
the equivalent linear form of Eq. 16 in kp together with Eqs. 14 and 15 con- 
itute a system with a unique solution; the solution for the original system of 
uations also must be unique. 

Determinants will be used to investigate the simultaneous solution of the 
nection and loop equations for the balanced network with equivalent linear 
anches. For the jth root of the solution 


which Qj is abranch flow, Dj defines the determinant containing the equation 
nstants and D isthe determinant of the equation coefficients. In Dj the con- 
mt for each loop equation would be a zero (except when a booster pump is 
erposed ina loop). However, loads and inputs are non-zero constants in 
junction equations for an active network and consequently Dj is always un- 
lal to zero. Because Dj is unequal to zero, the system of equations is non- 
nogeneous and it is required that D be unequal to zero for a unique solu- 
1.33 D will equal zero in this application only when any two of its rows or 


33 “The Mathematical Solution of Engineering Problems,” by M, G. Salvadori and K, 
i McGraw-Hill Book Co., New York, 1948, 
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columns have proportional elements. The order of Dis n= L+J-1= 2, 
requiring representation for each and every branch flow in each of its rows. 
For a network with two or more loops a distinct and independent set of zerc 
values characterizes each equation, precluding the possibility of proportiona 
determinant elements, insuring a non-zero value for D (and Dj) and conse- 
quently a unique solution. The unique solution will yield a single, real value 
for the root characterising each branch flow. From preceding arguments the 


TABLE 10,—EXAMPLE OF EQUIVALENT LINEAR NETWORK DETERMINANTS 
(Balanced network of Case I, Fig. 2) 
a ee a ee eee 
Eee ——————————————————— ee 
For the system of seven equations in seven unknowns, where h = ky| Q|(Eq. 16), the 
equations in terms of each branch are: 


; ‘ Con 
b Sia ere aa Pani inna: DURE E SGA: siceene ci Deicide 
+, =@. = Qt = 40 
+Qo -Qor = 30 
+ Qor + Qa + Q, = 110 
-Qq +Qqr = 201 
-Q, +Qer = 451 
+0,0240Q, +0.0146Q,1 -0.0180Qq -0.0259Qqr = ¥ 
+0.0016Q, +0,0180Qq +0.0259Qqr -0.0118Q, -0.0246Q..= 
For branch b: 
400 = 0 0 21 0 0 
300 re a 0 0 0 0 
1100 0 +1 +1 0 +1 0 
200 0 0 oy +1 0 0 
450 0 0 0 0 A +7 
0 +0.0240 +0.0146 -0.0180 -0.0259 0 0 
0 0 0 +0.0180 +0.0259  -0.0118  -0.0246 
=p * 
+o =i 0 0 = 0 0 
0 +1 =1. 0 0 0 0 
0 0 +1 +1 0 +1 0 
0 0 0 ze +1 0 ie 
0 0 0 0 0 = +1 
0 +0.0240 +0.0146 -0.0180 -0.0259 0 | 
+0.0016 0 0 +0.0180 +0.0259  -0,0118  -0.0246 


From which Dp # 0, D # 0. Therefore only one value can exist for Qp (Qp=+1643 gpm 


corresponding solution with Eqs. 13, 14, and 15 also must be unique. Table. 
has been prepared to illustrate the non-zero, nonhomogeneous character of t 
determinants, using the Case I balanced network data of Fig. 2. 4 

Only a single source has been discussed, but the proof applies as we ih 
multiple sources with or without storage. For the special case of m = rt 
preceding proof is tacitly rigorous. 
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It appears that Camp and Hazen384 were the first to adapt linear resistance 
electric computing boards to the analysis of water distribution networks, by 
means of an expression with the same meaning as Eq. 16. Atrial and error 
procedure was used to find the set of ko values that would Satisfy the corre- 
sponding k values, using appropriate conversion factors to reduce electrical 
to hydraulic units. Much later, convergence procedures were developed35, 36 
that vastly reduced the number of trials. Strangely, no proof of solution unique- 
ness was offered. Analysis by electrical analogy is based on satisfaction of 
Kirchhoff’s laws (1845); with A.C. circuits, vector rather than algebraic sums 
must be used. Fich and Potter3? give a valuable discussion on the elements 
of network topology and also cite special types of electronic network compo- 
nents that are associated with homogeneous systems of equations (a solution 
in indeterminate form may exist). Van Valkenburg38 gives brief attention to 
resistive network analysis, but in his book and in a number of others reviewed 
no special effort has been made to prove general solution uniqueness, perhaps 
because of the comparatively small importance of purely resistive cases and 
emphasis on general circuit characteristics. 

It may be noted in passing that electrical resistive analogues, using either 
linear or non-linear resistance elements, can be manipulated to satisfy dif- 
ferent values of m for any or all pipes in the network. Reasons why the use 
of mixed values of m is normally unrealistic have been mentioned previously 
(3). Robertson has presented very convincing, comprehensive arguments for 
the use of a single value of m for analyses via digital computer. Lomax has 
called attention tothe fact that laminarflow may occur in some pipes of a net- 
work, for which m = 1, but this consideration should not be critical in other 
than very smallarterial networks or for unusually low minimum demand rates. 
 Bitoun’s demonstration of the validity of Eq. 1 is a valuable contribution. 
Several decades ago one of the standard design procedures consisted of re- 
ducing a network by the equivalent pipe method to a single pipe. In order to 
Ho this all loads had to be consolidated into a single load. The method was 
particularly useful for small network fire flow analyses. While this equivalent 
pipe approach is repeatedly mentioned in the literature, it appears that J. Do- 
land19 was the first to suggest the simplified version of expressing the flows 
in each pipe as a percentage of the input. This principle has been extended 
without proof to multiple loads inelementary network examples in several text 
pooks. While writing the paper, the writer considered developing Eq. 1 via 
pquivalent pipes but decided against it. In order to emphasize the fact that Eq. 
2 is restricted to a single value of m and proportional loads, the simpler case 
without equalizing storage was offered in three levels of complexity as an in- 
iroduction. Because all of the discussions relating to the equations indicate a 
slear understanding of these limitations, the approach used appears to have 
een successful. The demonstration of the validity of Eq.1 by Bitoun could be 


34 “Hydraulic Analysis of Water Distribution Systems by Means of an Electric Net- 
rk Analyzer,” by T.R. Camp and H.L. Hazen, Journal, N.E.W.W.A., Vol. 48, Decem- 
, 1934, pp. 383-407. 

35 “Network-Flow Analysis Speeded by Modified Electrical Analogy,” by H. A. Perry, 

D, E. Vierling, and R. W. Kohler, E.N.R., September 22, 1949, pp. 19-23, 

36 “Use of Alternating-Current Network Calculator in Distribution System Design,” 
y M. V. Suryaprakasam, G, W. Reid, and J. C. Geyer, Journal, A.W.W.A., Vol. 42, 
ecember, 1950, pp. 1154-1164, , 

37 “Theory of A-C Circuits,” by S. Fich and J. L. Potter, Prentice-Hall, Inc., New 

rk, 1958. 

38 “Network Analysis,” by M. E, VanValkenburg, Prentice-Hall, Inc., New York, 1955, 
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extended tothe storage case of Eq. 2. However, a boundary value approach will 
not lead to a full derivation of the latter. If the writer had abandoned the more 
complex case because his repeated attempts to arrive at a solution by direct 
derivation proved unsuccessful, Eq. 2 would not have been developed, albeit 
empirically via a similarity hypothesis. 

Returning to the preceding proof of solution uniqueness, if each Q and con- 
stant in the system of simultaneous’ linear equations for m = 1 is multiplied by 
any particular non-zero constant, the solutionis unchanged except that the roots 
are equal to the original roots times the constant. Bitoun’s demonstration for 
Eq. 1 is indispensable for non-integer values of m. 

System Analysis via Proportional Load Hydraulics.—The information in 
Table 11 is given in response to the request by Linaweaver, Geyer, and Wolff 
for a clarification of Table 4 flow ratios. The corresponding hydraulic grade 


TABLE 11.—EXPLANATION OF FLOW RATIOS IN TABLE 4, ROX. H. as 


Qpr' |QpR/ Qpw: |@pw/Qa} Rsw; 
=QpR 7 mgd 


9 9.45 -8.0 
2 9.45 +9.7 
4 6.6 4,2 
4D 9 6.6 45.3 
5A 26.9 e 18.9 -8.0 
5B 2 E 18.9 +9.7 
5C 17.4 13.2 -4,2 
5D A!) A 13.2 +5.3 
a) -5, _ -8.0 
2 +6, _ +9.7 
4 — -4,2 
9 _ +5.3 


4 (For Qs, minus sign indicates flow out of storage andplus indicates flow into stor- 
age.) (Non-proportional load: Chestnut Hill; both series A and B, 1.2 mgd; both series 
C and D, 0.8 mgd.) 


data for this Roxborough H.S. test series are given in Table 12. Forthis series, 
realistic water levels andstorage rates were assigned for the Rox. standpipes, 
taken thereafter as a reference or control. The other water levels are from 
the network analyses, using this datum. The pump total dynamic head require- 
ments for uniform daily pumping listed in Table 12 would be approximated 
quite well by the provisional pumps described at the bottom of the table [see 
also Figs. 5 to 7, (1a)]. The W.O.L. standpipe water levels are reasonable, 
realistic and satisfactory for a first phase design except for runs 6A and 6B.. 
However, these two runs are for an outage at W.O.L.P.S. on the maximum day, 
an emergency condition under a rare combination of circumstances. The. sys- 
Be layout in Fig. 7 was the first studied, and not the final arrangement adop- 
8 

{ 

As Linaweaver, Geyer, and Wolff have emphasized, the computations pre- 
sented for the Roxborough H.S. District are merely a first trial or feasibility 
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study, and “a thorough analysis would still be required”----“accounting for 
changes in levels of tanks and standpipes.” The procedure for balancing that 
they describe could not have been applied effectively for this particular study 
because the proposed W.O.L. pumping station and standpipe, new pumps for 
the Rox. pumping station, and new major arterial mains were all unresolved 
design factors at the beginning of the study. It should be noted that Table 4 
could be extendedover a considerable number of other combinations by means 
of Eq. 2, for a balance study or to investigate other combinations of inputs. 
Constant is of the opinion that the developments offered in the paper are of 
little value to the person concerned with improvements for a particular distri- 
bution system, apparently because Eq. 2 is restricted toa given loading distri- 


TABLE 12,—WATER LEVELS AND PUMP REQUIREMENTS FOR TABLE 4, ROX. H. S, 


Rox, Stand- 
pipes (Sp), 
W.L. (Con- 
trol) Feet 


Mesyistonal Pumps for Table 4 Study; Acceptable for Above, Special, and Future Demand 
‘Requirements: 


Rated Capacity, Mgd Rated T.D.H., Shutoff, 


Nol] No.2 Nows Feet Feet 
| Rox. P.S. 7 10 10 110 160 
| W.O.L.P.S. 4 10 10 170 260 
: - 


bution, network configuration, and set of flow coefficients. A given system de- 
sign can be reduced to the hydraulic characteristics of four basic components: 
network, pump performance, pumping station and suction source, and storage. 
These must be integrated and reconciled with the design loadings. If the net- 
work proper is the only basic component subject to modification in the design, 
-onstant’s position is well taken, This would be a rather special design case. 
Nevertheless, Eq. 2 could still be of inestimable value in balancing final and 
Jternative network schemes. Entirely too many designs have been based on 
Tagmentary system studies. A comprehensive system analysis is required to 
asure that the benefits claimed for the design recommendations offered are 
physically attainable. From this standpoint it is difficult to see how “the oper- 
tions end of distribution problems” can be considered independently from a 


ia 
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design alleged to be operable. Even in a perfunctory system analysis at least 
four normal consumption demand rates should be studied. Only the network 
analyses for two demand rates are necessary to obtain the constants in Eq. 2. 
The network hydraulics for other rates can then be computed using the equa- 
tion, precluding the necessity for more than two runs for that particular net- 
work. Use of Eq. 2 is particularly effective when two or more basic compo- 


TABLE 13,—EXAMPLE OF STORAGE-NETWORK-PUMP BALANCING USING 
PROPORTIONAL LOAD PRINCIPLE; BELMONT H. sf 


Time In-|Qq (Ave.,| Pumping} Q, (By | Storage|Storage dh via Z Actual 
terval | 9.3), in | Station /Trial for Eq. (2), Pump 
(Starting) mgd | Suction Sp Level $=0.135, Total 
at 6:30 Water =3.70, | Dynamic | Dynamic 
am) in i in feet°|Head, in| Head, in 
hours Feet d feet 


(3) 


0=1 9.6 114.7 
1-2 11.6 112.8 
2-3 11.7 112.4 
3-4 12.2 111.3 
4-5 11.7 111.3 
5-6 a Siu 111.3 
6-7 11,1 111.0 
7-8 10.6 111.3 
8-9 10.2 111.8 
9-10 10.3 a 

10-11 10.5 111.6 

11-12 11,3 110.5 

12-13 12.3 109.3 

13-14 10.8 109.7 

14-15 9.6 110.2 

15-16 9.0 110.9 

16-17 8.4 111.5 

CT=18 7.4 113.4 

18-19 6.1 116.5 

19-20 5.1 119.8 

20-21 5,1 120.3 © 

21-22 5.1 120.7 

22-23 5.1 121.1 

23-24 7.3 117.3 


Notes: 2 For 2.0 million gal storage in 25-ft, over a 1-hr interval. 
b Storage water level at 0 (6:30 a.m.) hours taken at El, 355.00 (Full). 
C See Table 2, 
d Col. (6) - Col. (3) + Col. (7). 
© At 24 hr, El, 355,09; Q, (ave.) for 24-hr = 9.31 mgd | 
f See Fig. 6(a) and Table 2, . 


nents of the system are subject to modification in the design, permitting ine 


vestigation of several alternatives inthe timethat would otherwise be expended 
“on only one combination. “i 


Proof of the operability of a given design is best demonstrated by means’ of 
some sort of system balance. As noted by Linaweaver, Geyer, and Wolff, a 
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detailed hour-by-hour balance is not always necessary or desired. However, 
the more detailed the balancing the more indispensable Eq. 2 becomes, pro- 
vided proportional loading can be assumed, 

The use of Eq. 2 in the balancing of a system can best be illustrated by 
means of a detailed example. A storage-network-pump balance of the Belmont 
H.S. District (Fig. 6(a), Table 2) has been presented (5), but without benefit 
of Eq. 2. The system balancing outlined in Table 13 uses the same demands, 
storage, network and pump as in a previously noted example (5). To simplify 
the details in Table 13, the pumping station suction loss was taken as a con- 
stant (for Qp = 9.3 mgd, that is, the ave.). In this way the clear well water 
level and suction loss could be combined in terms of a “pumping station suc- 
tion water level.” Other than this abbreviation, details are the same as pre- 
‘viously given by the writer (5), but the balance is more exact. 

Each hourly demand in Table 13 was balanced by trial in sequence. Each 
trial value of Qp leads to a uniquely different value of Qp/' Qg; approximately 
Sixty trial combinations were needed for the preparation of Table 13. Without 
Eq. 2 a considerable number of network analyses would have beenrequired for 
comparable exactitude. Recall that only two network analyses are needed to 
determine ¢ and n. 

A realistic system balance is possible only when the basic system compo- 
nents are compatible. The Philadelphia Water Department is developing a 
system balance digital computer program that will incorporate Eq. 2 and sub- 
routines for pump characteristics andpumping station suction loss. The pump- 
ing station hydraulics will be separated because pump characteristics and sta- 
tion losses are non-linear independent design factors. The comparative in- 
uence of diverse pump characteristics andstorage capacities inthe balancing 

or given networks will be investigated under a variety of demand schedules. 
The principal objective of the study is the determination of basic criteria for 
ptimum utilization of equalizing storage andpumping facilities. The study was 

gely inspired by E. Shaw Cole, F. ASCE (1b). It is rather obvious that this 
study would not be feasible in the absence of the generalized network param- 
eter, Eq. 2. 

Corrections and Additions.—A recent article by Kincaid39 should be added 
to the series given in the list of references (6). 
_ The last equation of the paper should be corrected as stated by Constant. 
The word “now” should be substituted for “not” in the first line of p. 59, HY8, 
1960. Inclusion of the following phrase enclosed by parentheses would clarify 
the sentence on page 127, HY9, 1960: “For systems without storage or with 
torage ator near the pumping station, proportional load hydraulics (via Eq. 1) 
(el s----- .” Reference 19 appeared in Vol. 33, No. 2, pages 234 to 236, 
in February, 1941. In the line preceding Eq. 9, on page 41 of HY7 in 1960, “In 
y loop---” should be changed to read, “At any junction---.” 


39 “Value and Use of System Records in Long-Range Planning,” by R. G. Kincaid, 
ournal, A.W.W.A., Vol. 52, February, 1960, pp. 277-283. 
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DRAG AND LIFT ON SPHERES WITHIN CYLINDRICAL TUBES? 
Closure by Donald F. Young 


DONALD F. YOUNG.19_Laursen’s analysis of the results of this investi- 
gation in terms of the conventional critical tractive force is indeed interesting. 
Although, as noted, his analysis is speculative, it does demonstrate how lift 
and drag data of the type obtained may be utilized in sediment transport prob- 
lems. The writer agrees that for the relatively large values of d/D used in 
this investigation the contraction of the flow is certainly significant, and that 
the study should be extended to smaller values of d/D that would be more char- 
acteristic of sediment transport problems. As pointed out in the paper, the 
same technique as was used inthis investigation can be used to extend the range 
bf values of the parameters. However, as the particle size becomes smaller 
it will become increasingly difficult to determine the conditions of impending 
motion. 

With regard to Laursen’s comment concerning the location of the lift it should 
i noted that as far as the equilibrium equations for the particle are concerned 
he resultant fluid force acting on the particle can be resolved into components 
yt any point along the line of actionof the force. The writer has arbitrarily re- 
solved the resultant fluid force into a lift that passes through the center of the 
phere and a drag that does not. It would, of course, be incorrect to arbitrar- 

y assume that both components pass through the center of the sphere. Re- 
tardless of where the line of actionof the lift is located, the condition still re- 
mains that at impending motion, in a direction normal to the pipe wall, the re- 
sultant fluid force must pass through the center of the sphere in order for ro- 
ation not to occur. 

The angle B referred to by Laursen can be readily obtained from Eqs. 1 and 
! and the values of Cp and Cy, given in Table 1. For impending motion N = F 
- 0 and therefore 


pos tants Erroieaees 5 editor rereted) 
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The remarks of Indri concernthe difference in flow conditions that exist be- 
een the case of solid bodies falling ina still fluid and the caseof bodies held 
equilibrium by a moving fluid. This is a distinction that is frequently over- 
ked. Certainly the velocity profile in the stream preceding the solid has an 
ect on the force exerted on the solid. Also, any turbulence in the approach 
ream is possibly significant. This is a problem that has concerned investi- 
tors working with wind tunnel models for some time. 


8 June 1960, by Donald F. Young (Proc. Paper 2526). 
10 Assoc. Prof., Dept. of Theoretical and Applied Mechanics, Iowa State Univ. of 
ice and Tech., Ames, Iowa. 
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For the investigation under discussion the experimental conditions utilized 
were similar to those normally encountered in studies related to the incipient 
motion of suspensions, that is, the particles were at rest and the fluid moving. 
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MODELS PRIMARILY DEPENDENT ON THE REYNOLDS NUMBER® 


Closure by W. P. Simmons, Jr. 


W. P. SIMMONS, JR. 17_ the clarification of concepts and nomenclature pre- 
‘Sented by Priest, Kolf, and Reitmeyer are important to those working with 
closed conduit flow problems. Priest correctly points out that Reynolds num- 
‘ber applies to viscous motion and not directly to turbulent motion. Kolf and 
Reitmeyer likewise point out a distinction and affirm the useof Euler relation- 
‘Ships for analyzing the model data. Chapple’s contribution of a related and 
‘interesting problem well illustrates a facet of Reynolds number significance. 

The paper emphasized that through an understanding of the Reynolds num- 
|ber, models of large closed-conduit systems, in which viscosity might be thought 
of as a predominant factor, may be made relatively insensitive to viscosity. 
‘This is predicated on the model being operated at Reynolds numbers greater 
jthan some critical value below which viscous effects are significant. Such 
jmodels are indeed analyzed with parameters of the Euler variety that are di- 
jmensionally similar but conceptually different from the Froude ones. Thus, 
the statement that “...the data may be analyzed using the same relations as for 
jmodels based on the Froude law” is misleading. It should be amended to read 
,..the data may be analyzed using Euler relations.” As is pointed out, this 
conceptual difference involves no change in model techniques on data analysis. 
A point deserving further attention is that of determining Reynolds numbers 
high enough so viscous forces can be disregarded. By definition, R involves 
ia characteristic length and flow velocity in the fluid system. These will have 
different significance in different systems. For example, for orifice flow ina 
full pipeline an R, based on orifice diameter and flow velocity of about 5 x 105 
appears sufficient.3 Similarly in a draft tube study an R, based on inlet con- 
ditions, of about 3 x 105 was satisfactory.18 No doubt many other instances of 
Satisfactory test results at R’s less that 1x 10° can be cited. However, if 

inimum or near minimum values are to be used, these critical values must 
be determined through experiment. Unless it is worthwhileto make such tests, 
it seems convenient and safeto work inthe range of the 1 x 106 value suggested 


in the paper. 


a June 1960, by W. P. Simmons, Jr. (Proc. Paper 2531). 
17 Hydr, Engr., Div. of Engrg. Lab., Bur, of Reclamation, Denver, Colo. 

18 «Hydraulic Model Studies of Martin Dam Draft Tubes,” by C. E. Kindsvater and R. 
2, Randolph, Jr., Transactions, ASCE, Vol, 120, 1955. 
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WATER EDDY FORCES ON OSCILLATING CYLINDERS2 
Discussion by Donald VanSickle 


DONALD VANSICKLE, ¢ A.M. ASCE.—This paper is a most welcome addi- 
tion to the authors’ series of papers concerning wave forces on piling. 8 The 
accumulation of such data as they have presented is essential before design 
criteria can be established that correctly account for the many factors involved. 

The writer is particularly interested inthe relationship between drag forces 
on piles with tandem spacings because, in another investigation, 9 he studied 
this effect for H-section piles in uniform flow. It is of interest tocompare the 
authors’ test results on oscillating cylinders, their previous results on cylin- 
ders with uniform motion, 8 the writer’s results on H-piles in uniform flow, 
and the results of an NACA investigation of cylinders in uniform flow.10 
To facilitate this comparison, the various test results are summarized in 

ig. 7, in which the data are presented as ratios of the drag onthe downstream 
pile of a tandem pair tothe drag for a single pile for various center-to-center 
pile Spacings. It was necessary to computethe ratios for the NACA tests using 
the given drag coefficients and a drag coefficient of 1.2 for the single cylinder. 
For the authors’ tests on cylinders withuniform motion, the ratios were deter- 
: from the mean values of the drag coefficientsfor the various velocities. 

he results for the H-piles were already in the formof drag coefficients. The 
authors’ results on oscillating cylinders are presented as the ratios of drag 
forces ; however, if the depths and velocities for the tandem cylinders corre- 
spond to those for the single cylinder, the ratios of drag forces should corre- 
ban to the ratios of drag coefficients. It should be mentioned that the center- 
o-center spacings for the square H-piles are pile widths, whereas those for 
ihe cylinders are cylinder diameters. 
__ An examination of Fig. 7 reveals the following points: 


1, All of the curves have essentially the same shape. The ratioof drag co- 
fficients increases with increasedspacing up to about four diameters spacing. 
or spacings in excess of four diameters, the ratio of drag coefficients re- 
ains practically constant. It would appear, then, that a center-to-center spa- 
ing of approximately four diameters is the critical spacing as far as drag 
rces are concerned. 


a November 1960, by Alan D. K. Laird, Charles A. Johnson, and Robert W. Walker 
Proc, Paper 2652), 

_ 7 Hydr, Engr., Turner & Collie Cons, Engrs., Inc., Heuston! Tex, 

8 “Water Forces on Accelerated Cylinders,” by A. D. K. Laird, C. A. Johnson, and 
P af Walker, Proceedings, ASCE, Vol. 85, No. 1, March, 1959. 
_ 9 an Experimental Investigation of the Drag Characteristics of a Model H-Section 
#ldge Pier,” by Donald R. VanSickle, thesis presented to the Univ. of Texas at Austin, 
oy in 1958 * in partial fulfillment of the requirements for the degree of Master of Science, 

The Interference Between Struts in Various Co teaesenees by David Biermann 

ad W. H. Herrnstein, Jr., NACA Report No, 468, 1933. 
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2. The results of the authors’ previous tests on tandemcylinders with uni- 
form motion are in substantial agreement with the results of the NACA tests. 
D. Biermann and N. H. Herrnstein10 suggest that the negative drag forces at 
close spacings are caused by “the vortices produced by the front cylinder 
(which) probably encircle the rear cylinder, impinging on the back surface 
with sufficient force to produce a forward reaction.” The directions of the 
eddies shown in Fig. 6 indicate that this is certainly a reasonable explanation. 
Both sets of test results further indicate, for spacings greater than four di- 
ameters, a leveling off of the drag coefficient for a cylinder with neighbor at 
from 0.25 to 0.40 the drag coefficient for a single cylinder. Biermann and 
Herrnstein19 advance the hypothesis that “the boundary layer of the second 
cylinder is madeturbulent by the rapid changes of velocity and direction of the 
flow behind the first cylinder.” The well-known fact that the drag coefficient 


Oscillating cylinders, 
by authors 


H-piles in uniform flow, 
by Van Sickle? 


Cylinders with uniform motion, 
by authors® 


Cylinders in uniform flow, 
by NACA?° 


Drag coefficient with neighbor 
Drag coefficient without neighbor 


Ratio of 


-0.2 


4 5 6 7 8 
Center-to-center spacing, pile diameters 


9 10 ll 12 13 


5 
FIG. 7.—VARIATION OF DRAG COEFFICIENT RATIO WITH PILE SPACING 


of a cylinder decreases to approximately one-fourth of its normal value whe! 
the boundary layer becomes turbulent11 lends much support tothis hypothesis 
3. The results of the study of tandem H-piles are shown as a shaded are: 
in Fig. 7. The lower line of the area represents a Froude number of 0.2 an 
the upper line represents a Froude number of 0.6. To avoid complicating Fig 
7, the curves for intermediate values of Froude number are not shown; how: 
ever, as the Froude number increases from 0.2 to 0.6, the curves shift upware 
in fairly uniform steps. The shape of these curves is very similar to thes , 
of the curves for cylinders in uniform flow. The negative drag coefficients f 
close spacings can be accounted for as suggested in point 2. Although, for sp 
cings greater thanfour diameters, a similar leveling off of the drag coefficier 


amiss Rouse, ed, Engineering Hydraulics, John Wiley & Sons, New York, 198 
p. 124, 
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or an H-pile with neighbor occurs at values of 0.60 to 0.80 the coefficient for 
. Single H-pile, no convenient explanation is available. The point of separation 
or the H-pile is essentially fixed at the sharp flange edge, and no shift down- 
tream, with a reduction in drag coefficient, occurs when the flow becomes 
urbulent. 

4. In contrast to these observations for tandem piles in uniform flow, the 
uthors’ data for oscillating cylinders indicate that no negative drag forces oc- 
ur. Furthermore, while the drag coefficient levels off for Spacings greater 
aan four diameters, it does so at the same value as for a Single cylinder. In 
ther words, as the authors point out, for oscillating cylinders the effect of the 
pstream cylinder is negligible for spacings greater than four diameters. 


‘It occurred to the writer that because the authors’ results were presented 
8 a ratio of the force on the cylinder with neighbor to the force ona single 
ylinder, the leveling off of the ratio at a value of 1.0 could be explained if the 
scillatory motion, for some unknown reason, caused the boundary layer for 
ue single cylinder to be turbulent. The occurrence of turbulent flow on the 
pwnstream cylinder of the tandem pair would then produce the turbulent drag 
irce, and a drag ratio of 1.0. The authors point out, however, that “the drag 
vefficients corrected to infinite cylinder values were the usual 1.3 to 1.4 for 
1e vertical cylinders extending through the water surface.” Because this is 
.e drag coefficient for laminar boundary layer conditions, it is evident that 
.€ previously described explanation cannot apply. 

Another alternative is that the downstream cylinder of an oscillating tan- 
m pair does not encounter turbulent flow conditions, whereas the downstream 
linder in uniform flow does. There would seem to be noobvious reason why 

ich a condition should exist. 
Experimental error would not seem to be a factor because the authors’ pre- 
ous work, using essentially the same equipment and procedures, gave results 
Substantial agreement with the NACA tests. 
There remains the further possibility that the hypothesis advanced by Bier- 
and Herrnstein to explain the reduced drag on the downstream cylinder 
not valid. The fact that their hypothesis cannot explain the reduced drag on 
e downstream H-pile adds weight to this argument. Such a conclusion, how- 
er, does not explain why piles in uniform flow have reduced drag on the down- 
ream pile at spacings in excess of four diameters, whereas cylinders in os- 
lating motion do not have reduced drag for these spacings. 


Do the authors know of any reason, or could they present a reasonable hy- 
thesis, why this situation should exist? 
Could they state whether or not they actually encountered negative drag 


ces in any of their tests on oscillating cylinders? 
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FISH PASSAGE THROUGH HYDRAULIC TURBINES2 


Discussion by J. F. Muir 


| 

imo. F. MUIR,” F, ASCE.—The research on the passage of young fish through 
ydraulic turbines, by Von Gunten and his associates, is one of the most im- 

,ortant projects arising from the 3.5 million dollar fisheries-engineering re- 

earch program initiated in 1951 by the Corps of United States Army Engineers, 

_ These recent tests furnish further evidence that the Corps of Engineers’ 

eeram, participated in by all the major fisheries organizations in the Pacific 
orthwestern states, has paid handsome dividends. A wealth of data has been 

jublished giving criteria for the safe passage of fish past dams.3,4,5 The 
ter believes that this program and further research will very likely show 

nat rivers can be developed for multi-purposes, including the preservation of 
dromous fish. 

The preliminary tests on the passage of fingerlings through small Francis 
nd Kaplan runners 12 in. in diameter inthe laboratory showed that cavitating 
ionditions in a turbine will kill fish. The rate of mortality is related to the 
nagnitude of partial vacuum present in the water passages. 

_ The full scaletests at Cushman dam confirmed laboratory findings showing 
at significant mortality does not occur among young fish subjected to an ab- 
pt reductionin water pressure from up to 200 psi to atmospheric. The tests 

howed a significant improvement in survival rates resulting from the reduc- 

on in cavitation in the turbine whenthe tailwater level rose 10 ft. The maxi- 

num mortality occurred when the turbine operated under partial loading, a 

ondition known to result in a worsening of cavitation phenomena. 

_ The writer had the opportunity of witnessing the Cushman tests and was par- 

;cularly interested in observing the fish for evidence of the bends, Experi- 

ents at the University of British Columbia (Vancouver, B. C.) have shown 

the bends are likely to occur among fish subjected to a partial vacuum of 
gh intensity for periods of time considerably greater than the time of pas- — 
ge through the Cushman turbines. The writer saw no symptoms of the bends 
long the Cushman fish after the tests and believes that few fish were killed 
* seriously injured by the effects of apartial vacuum unaccompanied by cavi- 


ca 


& May 1960, by Glenn H. Von Gunten (Proc. Paper 2812). 
2 Head, Dept. of Civ. Engrg., Univ. of British Columbia, Vancouver, Canada. 

3 “Progress Report on Fisheries-Engineering Research Program,” Corps of Engrs., 
§. Army, Pacific Div., November, 1956. 

“Progress Report on Fisheries—Engineering Research Program,” Corps of Engrs., 
3, Army, Pacific Div., July, 1960. 

Annotated Bibliography on Fish Passing Facilities, Committee on Power Project 
g and Design, Power Div., ASCE, September, 1960. 
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One interesting feature of the Cushman tests was the efficiency of the large 
net used totrap fingerlings. The mouth of the net completely covered the draft 
tube exit and trapped approximately 90% of the fish passing through the turbine 
during the tests. Ivan J. Donaldson and his colleagues are to be commended 
for their pioneering work on the design of large nets for this purpose. 

The research of Von Gunten and his associates points to a further program 
designed to solve the problem of successful passage of fish through turbines. 
One important aspect of this problem is the determination of relationships be- 
tween turbine specific speed, the cavitation parameter sigma, and fish mor- 
tality. The Cushman turbines are placed at an elevation about 1 ft above tail- 
water level at high tide. In a future plant similar to Cushman, sigma could be 
increased more than 50% by lowering the turbines about 20 ft. One wonders 
about the extent of improvement in fish survival rates with such a substantial 
increase in the value of sigma. 

In the light of the author’s findings, it would be of interest if he would dis- 
cuss the changes in operating conditions called for at existing plants on the 
Columbia River system. During the migration period, is it nowstandard prac- 
tice to avoid turbine operation under low loadconditions inplants in which young 
fish pass through turbines on their way to the sea? 

It wouldalso be of interest tohave Von Gunten’s comments on the possibility 
of safe passage of migrants through the Grand Coulee and Chief Joseph tur- 
bines. It would now appear that the rehabilitation of the salmonruns in the up- 
per Columbia River is within the realm of future possibilities. 
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Journal of the Hydraulics Division 
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“Uniform Water Conveyance Channels in Alluvial Material” by cares B. 
imons and Maurice L. Albertson, May 1960, Proc. Paper 2484. 


64. Under the heading “Design Procedures: Estimating W Knowing P or 
Versa” in lines 1 and 5 change the word wetter to wetted. _ 
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2725(ST1)°, 


: (WW1), 2727(EM1), 2728(EM1), 2729 (WW1), 2730 (WW1), 2731 (EM1), 2732(6M1), 2733 (Ww1), 
(S) a aaneunt), siete, 2737 (PL1), 2738(PL1), 2739(P11), 2740(PL1), 2741(@M1), 2742672), 
(EM1), 2744(WW1), 2745(WW1), 2746(6Mi), yo ela 2748(BM1), 2749(WW1), 2750(WW1)°, 2751 
(SM1)°, 2753 (PL1)°, 2754(6T2)¢, 2755 : 

scat aes, 2758(AT1), 2769(CO1), 2760(HY¥2), 2761(IR1), 2762(0R1), 2763(HY¥2), 2764 
Bratati) grrsUR), 2176(Hy2), 217TH), BT 19(SAs}, ATOLSTS), Ar —OGIT2), 2TEL(RY), 2182 
¢ Perk 2777 , ; ; 
“ PrastaAsle. saresicon, 2785(CO1)°, 2786(IR1)°, 2787(ST3)¢, 2788(AT1)°, 2789(HW1), 


2795(SM2), 2796(SM2), 2797(SM2), 2798 
), 280 2804(SM2)°, 2805(ST4) 


2), 2791(SM2) 


), 2807(wwa), 2608(HY3), 2609(ww2), 2810(H'¥S), 2011(ww2), 2812(aya), 2619(Ww2), ° 


4a 83 , 2816(H¥3), 2817(HY3), 2818(SA3), 2819(WW2), 2820(SA3), 2821(WWw2), 2822 
223 (HY), 2624(8AS), 2625(HVS), 2026(5A3)°, 2627(H¥S)°. 
grouped by divisions, we 


~ 


MOLINEAUX 
E A. ELSENER — 


> ; 2 
; vad, Nope 


